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A theoretical analysis of the electrothermal ramjet concept of launching space bound pay- 
loads at high acceleration levels is described. A simple system design which could be used to 
implement the concept is identified. The operation of this system is broken down into two 
phases. In the light gas gun phase the payload is accelerated to the velocity at which the 
second or ramjet phase can commence. Preliminary idealized models of operation are examined 
and shown to yield overall energy efficiences for a typical earth escape launch of 60 to 70%. 
When shock losses are incorporated these efficiences are still observed to remain at the 
relatively high values of 40 to 50%. 


17 Key Words (Suggested by Author(s)) 

Electrostatic Thruster 
Hollow Cathode 
Electrothermal Ramjet 

18 Distribution Statement 

Unclassified - Unlimited 

19 Security Classif (of this report) 
Unclassified 

20 Security Classif (of this page) 
Unclassified 

21 No of Pages 
149 

22 Price" 


‘ For sale by the National Technical Information Service, Springfield Virginia 22161 


NASA-C-168 (Rev 10-75) 


1 





TABLE OF CONTENTS 


Topic Page 

Abstract i 

Mercury Hollow Cathodes 1 

Ion Beamlet Vectoring 2 

Controlling Ion Currents to Discharge Chamber Surfaces .... 2 

Ion Extraction using Closely Space Grids 4 

The Influence of Stray Magnetic Fields on Ion Beam 
Neutralization 5 

Enhancement of Ion Beam Currents Through Space 

Charge Compensation 7 

Introduction 7 

Electron Injection Between the Acceleration Grids ... 8 

Electron Injection from Upstream of the Screen Grid . . 16 

Electron Injection from both Sides of the Screen Grid . ]8 

Multiple Energy Level Electron Injection 21 

Conclusions 34 

The Electrothermal Ramjet 36 

Theoretical Model of Diffuser Shock Losses 38 

Results for a Typical Supersonic Launch Cycle 42 

Conclusions 48 

Constant Pressure Acceleration of the Electrothermal Ramjet. . 49 

Introduction 49 

Uniform Acceleration Flows - General Description. ... 50 

Piston Driven Isentropic Accelerator 54 

Pistonless Single Chamber Accelerator 63 

Conclusions 70 

References 72 

Appendix 74 

A Model for Mercury Orificed Hollow Cathodes: 

Theory and Experiment 75 

Ion Beamlet Vectoring by Grid Translation 86 

The Flexible Magnetic Field Thruster 98 

Ion Extraction Capabilities of Very Closely 

Spaced Grids 108 

ii 



Table of Contents (Continued) 


4 


Topic 

The Influence of Stray Magnetic Fields on 
Ion Beam Neutralization 

The Electrothermal Ramjet 

Distribution List 


Page 


118 

128 

139 




T 


111 



LIST OF FIGURES 


Figure No . Ti tl e Page 

1. Configuration without Electron Injection 9 

2. Configuration with Intra-grid Electron 

Injection 9 

3. Property Profiles for Optimal Intra-electrode 

Electron Injection 14 

4. Effect of Injected Electron Current on 

Extracted Ion Current 15 

5. Configuration for Electron Injection from 

Upstream of the Screen Grid (Plane "A") 17 

6. Configuration for Simultaneous Upstream and 

Intra-grid Electron Injection 17 

7. Property Profiles for Optimal Injection of 

Two Electron Groups 22 

8. Configuration for Multi-energy Electron 

Injection 24 

9. Effect of Number of Optimally Injected Electron 

Beams on Extracted Ion Current 29 

10. Injected Electron Current Density Required to 

Yield the Maximum Extracted Ion Current Density . . 31 

11. Effect of First Stage Electron Injection Current 

on Extracted Ion Current 33 

12. Electrothermal Ramjet Schematic 37 

13. Typical Ramjet Launch Profile for a 10 kg Mass 

Accelerated at 30,000 g 44 

14. Uniform Acceleration: Time vs. Distance 52 

15. Accelerator Design 52 

16. Efficiency of Piston Driven Accelerator 61 

17. Rate of Heat Addition: Piston Driven Accelerator. . 61 

18. Pressure and Temperature: 

Pistonless Accelerator 68 

iv 



List of Figures (Continued) 


Figure No . Title Page 

19. Efficiency of Pistonless Accelerator 68 

20. Rate of Heat Addition: 

Pistonless Accelerator 69 


v 



MERCURY HOLLOU CATHODE RESEARCH 


Dan Siegfried 

An extensive experimental investigation of mercury hollow cathodes has 
been completed. A model describing the basic physical processes taking 
place within the cathode has been developed. It can be used to predict, 
to first order, important cathode operating parameters such as emission 
length and insert temperature. The analytical formulation of the model is 
based on the concept of an idealized "ion production region" which is defined 
as the volume circumscribed by the emitting portion of the insert. The 
energy exchange mean free path for primary electrons is used as a criterion 
for determining the length of this region. An ion production region aspect 
ratio (dia.-to-length) of two is suggested as a design criterion for 
minimizing keeper voltage. The model accounts for electrons produced in 
the ion production region both by surface emission and by volume ionization. 
Surface and volume energy balances are used to predict plasma density and 
plasma potential in this region. An empirical relation is presented which 
can be used to estimate cathode internal pressure (a necessary input to the 
model) using the mass flow rate, discharge current and cathode orifice 
diameter. Calculations based on the model are compared with experimental 
results. 

This research is described in some detail in Reference 1, which was 
released recently, and it will therefore not be reproduced here. In 
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addition a summary of this work is described in a recently completed paper 
which is included in Appendix A of this report. 



ION BEAMLET VECTORING 


J. M. Homa 

An experimental study of ion beamlets emerging from two-grid optical 

systems that were vectored by translating the screen/accelerator grid holes 

3 

relative to each other was completed. A paper describing this work is 
included in Appendix A of this report. A more detailed description of 
this research will also be presented in thesis form and present plans call 
for this work to be released as a NASA Contractor report. The results of 
this research can be used to design accelerating systems in which the 
direction and focus of emerging beamlets are important. Deflection and 
divergence angle data were collected for two-grid systems as a function 
of the relative displacement of these grids. At large displacements, 
accelerator grid impingement was observed to become a limiting factor and 
this defined the useful range of beamlet deflection. Beamlet deflection 
was found to vary linearly with grid offset angle over this range. Values 
of deflection-to-offset angle ratio and the useful range of beamlet deflec- 
tion were determined as functions of grid hole geometries, perveance levels, 
and accelerating voltage levels. The divergence of the beamlets was found 
to be unaffected by deflection over the useful range of beamlet deflection. 

CONTROLLING ION CURRENTS TO DISCHARGE CHAMBER SURFACES 

J. R. Brophy 

The study of mechanisms that enable a designer to control the loss 
rates of ions to various discharge chamber surfaces in an ion thruster 
continued to focus on the use of the flexible magnetic field ion source 



during the past year. The resulting research has been described in a paper^ 
included in Appendix A of this report. 

The flexible magnetic field thruster is a unique research tool for 
studying the behavior of direct current electron-bombardment ion thrusters. 

It utilizes a long wire anode shielded by a circumferential magnetic field. 

The magnetic field around the anode, which is produced by passing a direct 
current of up to 150 amperes through the wire, restricts primary electron 
flow from the discharge plasma to the anode. Different magnetic field con- 
figurations can be created by routing the anode wire in various ways through 
the discharge chamber. Divergent, cusped and multipole magnetic field con- 
figurations have for example been created in this manner. The thruster is 
also designed so ion currents to various internal surfaces can be measured 
directly and these measurements facilitate calculation of the distribution 
of ion currents within the discharge chamber. Experiments indicate that the 
distribution of ion currents inside the discharge chamber is strongly dependent 
on the shape and strength of the magnetic field but independent of the dis- 
charge current, discharge voltage, and neutral flow rate. Through proper 
selection of the magnetic field shape and strength, the fraction of the ions 
produced in the discharge chamber that is directed toward the screen grid has 
been increased to over 50$. Measurements of the energy cost per plasma ion 
indicate that this cost decreases with increasing magnetic field strength 
due to increased anode shielding from the primary electrons. Energy costs 
per argon plasma ion as low as 50 eV were measured. The energy cost per beam 
ion was found to be a function of the energy cost per plasma ion, extracted 
ion fraction and discharge voltage. Part of the energy cost per beam ion 
goes into creating many ions in the plasma and then extracting only a fraction 
of them into the beam. The rest of the energy goes into accelerating the 
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remaining plasma ions into the walls of the discharge chamber. Measurement 
of ion fluxes across a virtual anode surface appear to indicate that ions 
cross this surface with velocities approaching their random thermal velocity 
rather than the Bohm velocity. 

ION EXTRACTION USING CLOSELY SPACED GRIDS 
Dean Rovang 

Experimental work has shown that the ion current density that can be 
extracted through two grid systems can be enhanced if small diameter aper- 

5 

tures in closely spaced grids are used. This work is described in a paper 
contained in Appendix A. It has focused on accelerator systems with small 
screen hole diameters (<_ 2.0 mm) operating at net-accelerating voltages of 
100, 300, and 500 V. The apparatus developed for this study is shown to be 
well suited for measuring the impingement-limited perveance, electron back- 
streaming, and electrical breakdown characteristics of two-grid ion accel- 
erator systems. Results suggest that the impingement-limited perveance is 
not dramatically affected by reductions in screen hole diameter to 1.0 mm. 
Impingement-limited performance was found to be dependent on the grid 
separation distance, the discharge-to-total accelerating voltage ratio, and 
the net-to-total accelerating voltage ratio. Results obtained using closely 
spaced grids having small hole diameters and operating at fixed low values 
of net accelerating voltage show that the onset of high impingement current 
occurs at progressively higher current density levels as the net-to-total 
accelerating voltage ratio is increased . This unexpected result suggests a 
new philosophy of high current density grid operation at high net-to-total 
accelerating voltage ratios. When grids are operating at this condition 
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it appears the onset of electron backstreaming ultimately determines the 
extractable ion current density. Beam current densities as high as 25 mA/cm 2 
were obtained using grids with 1.0 mm diameter holes operating at a net 
accelerating voltage of 500 V. 

THE INFLUENCE OF STRAY MAGNETIC FIELDS ON THE ION BEAM NEUTRALIZATION 

Feng Yu-Cai 

Work has continued on the study of the influence of magnetic fields 
in the vicinity of the neutralizer on its capacity to couple to and 
neutralize an ion beam emerging from a thruster. This work is described 
in a paper 6 also included in Appendix A. This paper describes the experi- 
mental investigation in which a comparison of the ion beam neutralization 
characteristics of a local neutralizer (within ^ 5 cm of the beam edge) 
with those associated with a distant one ('v 1 meter away from thruster) 
was made. The influence of magnetic fields in the vicinity of the neutral- 
izer cathode orifice which were either parallel with or normal to the 
neutralizer axis were studied. Plasma property profiles which reflected 
the influence of the magnetic fields were measured. The results indicated 
that magnetic fields at the region of a neutralizer cathode orifice in- 
fluence its ability to couple to the ion beam. They show that there is a 
potential jump from the neutralizer cathode orifice to the plasma which 
exists close to the orifice. The magnitude of this potential jump increases 
as the axial component of magnetic flux density increases. A magnetic field 
perpendicular to the neutralizer axis induces a potential rise a few centi- 
meters further downstream of the neutralizer cathode. Test results indicate 
the path length from the neutralizer to the beam along magnetic field lines 
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also has an effect on the beam-to-neutral izer plasma potential difference 
but this effect was typically smaller than the other two effects. The 
results suggest that in order to minimize the adverse effects of stray 
magnetic fields, a region as free of magnetic fields as possible should be 
selected for placement of the neutralizer. 



ENHANCEMENT OF ION BEAM CURRENTS 


THROUGH SPACE CHARGE COMPENSATION 


Feng Yu-Cai 


Introduction 

Langmuir' 7 has discussed the possibility of increasing the space charge 
limited electron current between two electrodes by injecting positive ions 
into the electron space charge and thereby partially neutral izating its 
detrimental effect on electron extraction. In the following work the 
neutralization of the positive ion space charge induced by injecting negative 
changes into the region between the grids of a positively charged ion accel- 
erator has been studied theoretically. If it is possible to do this the 
ion current density of the extracted beam could be increased. This could 
lead to an increased thrust density for electrostatic thrusters. 

Several options for charged particle injection are considered and the 
results suggest that the extracted ion beam density can be increased greatly 
by negatively charged particle injection. The ratio of the extracted beam 
current density to the conventional space charge limited value can be in- 
creased through the use of charged particle injection by a factor 30 or 
more, depending on the current density and velocity distribution of the 
injected negatively charged particle beams. The ratio of injected particle 
current density to the ion beam current density required to effect this 
change is of the order of the square root of the ratio of the extracted ion 
mass to the injected particle mass. In the following development the injected 
negatively charge particles will be designated as electrons and the extracted 
beam will be designated positive ions. The analysis is carried out in general 
terms however and the injected particles could just as well be negative ions. 
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Further, the results are applicable to the case where positive ion injection 
is used to enhance the extracted current densities of negative ion beams. 


Electron Injection Between the Acceleration Grids 

Consider an infinite plane screen grid "A" at potential and a 
parallel plane cathode (accelerator grid "C" at zero potential located at 
a distance a from "A" as shown in Fig. 1). Assume an infinite supply of 
ions without appreciable initial velocities is available at the screen grid. 
If these ions are extracted at the space charge limit, the potential 
gradient at the screen grid will be zero, and the current density is giveh 

g 

by the Child-Langmuir equation 


!f° fZe , 1 / 2 V /2 

9 ^m.j ' a? 


( 1 ) 


where e is the ion charge, e Q is the permittivity of free space and nn is 
the mass of the ions. 

Let us now consider the effect of introducing negative electrons with- 
out initial velocity, uniformly distributed over a plane "B" which is at a 
distance b from "A" as suggested in Fig. 2. Between "B" and "A" there is 
as a result of this injection an electron current per unit area j . Because 
of the partial neutralization of the ion space charge, the ion current per 
unit area from the screen grid will increase to a new value, j.. We 
assume that the ions and electrons do not collide with gas molecules nor with 
each other and that no appreciable number of ions or electrons is lost by 
re-combination during the passage between the electrodes. Let v.(x) be the 
velocity of the ions at any point p, which is a distance x from the screen 
grid "A". The ion current density at p is then p^(x). The corresponding 
quantities for the electrons are denoted by the subscript e. 

Then in region A-B one has 
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PLANE "A" PLANE "C" 



0 a * 

Fig. 1 Configuration without Electron Injection 



Fig. 2 Configuration with Intra-grid Electron Injection 
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p i (x) = v^TxT 


where 


2(V A -V)e 
A 

1 »' m. 

The electron charge density in this same region is given by 

j. 


( 2 ) 


(3) 


p e (x) = FTxT 


(4) 


where 

v e (x) = 

V is the potential at a point p and Vg is the potential at plane "B" . 
Poisson's equation may be written in one dimension as 




2(V-V B )e 

iii 


(5) 


d 2 V 


•.iBE 

m e 


! °V5 


V)e 


(6) 


Using the definitions, 

, <p = V/V A , A = x/a (7) 

and combining Eqs. (1) and (6) one obtains 




o 


where <f>g = Vg/V^ as suggested by the definition Eq. (7). 

In the space charge limited case the boundary condition at the screen 

9rid is (&) A , 0 ■ 0 • < 9 > 
Applying this same zero gradient condition at plane "B" for the integration 
of Eq. (8) gives 

(t-tg)** + - a(l-<P B ) h . 



( 10 ) 
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The normalized position of plane "B" (X R ) can be determined by integrating 


this expression to obtain: 

ri* J 


x -1 - 

B 4 j 


T~] f + 


(l-4>) % - <*(!-*„ ) h V 


In the region B-C of Fig. 2 only ions are present and their charge 


density is given by 


P«(x> =TTT 


1' ' v i lxj 


V^x) =' 


2 e(V.-V) 


Poisson's equation is given in this region by 


’o 


Integrating this equation one obtains 


/ \ 2 i. 

( djM _ lj^ _i /. , p 

1 dx j - ~ jr (1_ + ) + c • 

\ / 


The constant of integration for Eq. (15) is again selected so that 4f- = 0 


at <p = <f> B and the equation becomes 
ja. n r J.- t h r 


d(j) _ 4 f i 1 T n /i , 

dX ' ' 3 [~ (1 '^ " a(1 "V 


Integrating this equation from X = X^ to 1 one has 


T [j~] 


Adding Eqs. (11) and (17) we obtain the ratio of extracted beam current 
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density with injected electrons to the space charge limited ion current 
without electron injection. It is given by 


ii -_a 
\ ' ’ 6 


| f j^ot( 4>- tf> B ) 2 + 

A 

P B r ^ 

I [(WV 2 -« 


(l-*)* 8 - «(l-* B )^p dcf. 


dV] d(f> 


(18) 


In the space charge limited case the boundary condition at plane "B" is 


( d± 

\ dA 


= 0 


(19) 


B 


and from Eq. (16) we see this means 

a = 1 or 


' e 

Eq. (18) can be integrated to give 

| /*1 t t _j s 

I + (!-+)** •- (Wg) 1 * ] <4 + | (1- VTT b ) 3/2 

J* n 


( 20 ) 


ii -.I 
\ ' 16 


(21) 


+ 40-* B r* (i./w B )^ 


For the case where <j> B = 1 , Eq. (21) gives 




( 22 ) 


This means that if plane "B", from which electron emission occurs, coincides 
with the screen grid, the extracted beam current density cannot be increased 
by electron injection. For the case where <j> B = 0, i.e. locating the electron 
emission plane at the accelerator grid, Eq. (21) can be integrated 
numerically to obtain 
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i- = 1.87 . 


(23) 


This is the same result obtained by Langmuir for the case where ions were 

8 

injected into the electron space charge. Differentiating Eq. (21) with 
respect to (jig, and setting the result equal to zero one obtains the value 
of <t>g that results in the maximum ion current. This occurs when <j> B = 0.4 
at which potential Eq. (21) gives 
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Numerical integration of Eqs. (10) and (16) has been accomplished for 
this optimal case where 4> B = 0.4 and a = 1. This integration yields the 
profile of electrical potential which can be used in turn to obtain the 
potential gradient, and ion velocity profiles. Fig. 3 shows the normalized 
versions of these profiles. For this optimal case the electron emission 
plane is located at b/a = 0.43 as shown in the figure. Either increasesor 
decreases in the optimal value of b/a cause the peak electric field to 
decrease from the value shown in the lower plot and this in turn causes the 
ion current density to decrease from its maximum value. It is noted that 
the electric field has been normalized using the ratio of the intra-grid 
potential difference to the grid separation distance. This is the electric 
field that would exist between the grids if the ion and electron densities 
between the grids were negligible. The electric field for the case of 
electron injection is observed to be over 50% above this reference electric 
field. Figure 4 shows how the maximum positive ion current that can be 
extracted through the grids varies with the electron injection current and 
the potential of the plane from which the electrons are emitted. The ex- 
tracted ion and injected electron currents have been normalized using 



ELECTRIC FIELD ( Ea/V A ) NORMALIZED POTENTIAL OR VELOCITY 


POTENTIAL ( V/V A ) 




b/a POSITION ( x /a ) 


T5 


Fig. 3 Property Profiles for Optimal Intra-electrode 
Electron Injection 



NORMALIZED EXTRACTED ION CURRENT 
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respectively the conventional space charge limited ion current and the mass 
weighted electron current of Eq. (7). The figure suggests the highest ion current 
enhancement is realized at when emission occurs at tj> g= 0.4 as suggested pre- 
viously. It also shows that a reduction in the injected electron current 
from the maximum value (a=l ) is accompanied by a dramatic decrease in the 
extracted ion current. This points out the drawback of electron injection 
because a=l implies a large electron current (Eq. 20) and this in turn 
suggests a substantial input power to the electrons. 

Electron Injection from Upstream of the Screen Grid 

Let us now consider the effect of introducing negative electrons which 
are shot into the region between screen and accelerator grid from the upstream 
side of the screen grid with their velocity directed toward the accelerator grid. 
Because of the retarding field seen by the electrons as they enter the region 
between the grids, their velocity will decrease and at some distance from 
"A" (at plane "B" shown in Fig. 5) their velocity will become zero. They 
will then be accelerated back toward the screen grid. As before, we assume 
that the ions and electrons collide with neither the gas molecules nor with 
each other and that no appreciable number of ions or electrons is lost by 
recombination during the passage between the electrodes. As a consequence 
of electron injection to neutralize the ion space charge, the ion beam 
current density extracted by the grids will reach a value j^ . 

The development of this problem parallels exactly the problem just 
solved. The only difference in the development being that the density 
the injected electrons at any given location is twice that for the case 
where electrons are emitted from a plane between the grids. 

Solution of this problem shows that there is an optimum injection 
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0 b a 

Pig. 5 Configuration for Electron Injection from Upstream of the 
Screen Grid (Plane "A") 


ARBITRARY PLANE "F" 

PLANE (POTENTIAL V F ) 



0 f b a 


Fig. 6 Configuration for Simultaneous Upstream and 
Intra-grid Electron Injection 
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energy corresponding to the location of the plane "B" where the velocity 
goes to zero. As one would expect 


P B opt 


= 0.4 


and 


<Wmax ” 4 ' 6 ' 


Hence the same 4.6 fold increase in the ion current density that was 
obtained using intra-grid electron injection can be effected by optimal 
electron injection from upstream of the grids. Numerical integration of 
the equations in this case can again be used to determine potential, 
potential gradient and ion velocity profiles. These profiles are identical 
to those shown in Fig. 3. The plot of extracted ion current as a function 
of injected electron current is the same as Fig. 4, except for the fact 
that the peak ion currents occur at a normalized injected electron current 
(a) of one half. Hence one can conclude that the effect of electron injection 
is the same regardless of the means of getting the electrons into the ion 
acceleration region. The advantage of upstream electron injection is the 
required electron current is half of that required for intra-grid injection 
(a = 1/2 vs . a = 1 ) . 


Electron Injection from both Sides of the Screen Grid 

It is possible to inject electrons from both upstream and downstream 
of the screen grid simultaneously in order to neutralize the space charge 
within an ion beam. Figure 6 shows a physical model of this situation in 
which electrons are injected from upstream of the screen grid at plane "A" 
with a current density j e ^ and a velocity sufficient to enable them to reach 
plane "F". These electrons pass both directions through the region A-F as 
shown in Fig. 6. Electrons are also released at plane "B" with zero velocity 
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and they pass through region A-B going toward plane "A" as shown. As in 
the cases considered previously collisions are neglected and one computes 
the electron and ion charge density in each region as a function of the 
velocity of the electrons injected from upstream, the location of electron 
injection from between the grids and the current densities of the two in- 
jected electron beams. The charge densities are then used in the Poisson 
Equation describing each region, the resulting equations are integrated 
and potentials and potential gradients are matched at planes "B" and "F". 
The basic procedure is similar to the ones followed for the cases just 
discussed. As a result of this analysis one obtains the following expres- 
sions for the nondimensional ized locations of the planes "B" and "F" and 
the ratio of the ion current to the ion current without electron injection. 



2 
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The potential gradient expressions describing each of the regions of 
Fig. 6 are: 
for region A-F 

N l h i 


+ ot (1 -4> r ) 
2 


1 

hi ' 


( 28 ) 


for region F-B 


Ha = " f [ 7 “] + ~ [ 2 “ 1 0-4 > f )' 2 + ^(Wg)' 2 ] 


h 


(29) 


for region B-C 


^.i[y^^. K0V + O2(lV ] 


(30) 


Application of the space charge limited conditions at planes "B" and "F" 
yields the normalized injected electron currents from Eqs. (28) to (30). 


J i iV m e 


and 


20-4-pr 2 ( < f > F ~ < t , B ) !$ 




m. 


iptj" (2a 1 -l)(l-<j) F )' 2 

° 2 ’ m e ~ Up-^rMl-^r 2 * m e 


(31) 


(32) 


Numerical integration of Eqs. (28) to (30) yields the potential variations 

through the ion acceleration region as a function of a , a , <j) c and <j> c . 

l 2 e ‘ 

The maximum ion current density is observed for 


a x = 0.42, 


a 2 = 0.30, 
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4»p = 0.4 

<j)g = 0.16 = <J>p . 

At the condition corresponding to these parameters the normalized ion current 
is 


£i 

J io 


11.7 


and the normalized total injected electron current (from both sources) is 
required to be 

^■ 0 . 7 ^ - 

J i V m e 

Hence in this case the ion current density can be increased over an order of 
magnitude above the value without injected electrons. As before, however, 
very large injected electron currents are required. Figure 7 shows the 
potential, potential gradient and velocity profiles associated with this 
maximum ion current density operating condition. In this case the electric 
field peak near plane "A" (x/a = 0) is observed to be greater than it was for 
the case of electron injection at a single energy. It is this higher electric 
field that induces the higher ion current densities observed. It is noted 
that optimal electron injection from upstream at two energy levels produces 
profiles that are identical to those shown in Fig. 7. In this case however 


a 2 = 0.15 and hence 


j i fmT * / m. 


Multiple Energy Level Electron Injection 

The logical extension of the preceding section of this paper is the 
analysis of the case where electrons are injected into the ion acceleration 
region at various energies. By doing this the injected electrons can induce 
space charge limited flow boundaries at an arbitrary number of locations 
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Fig. 7 Property Profiles for Optimal Injection of Two Electron Groups 
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within the ion acceleration region. 

The basic one-dimensional model of this process is shown schematically 

in Fig. 8. The screen grid is represented by the plane "A" and is at a 

potential V^. The accelerator grid represented by plane "C" is assumed to 

be at zero potential and is located a distance a from plane "A". Electrons 

are ejected from plane "A" at velocities v p , v , v . . .v all 

% e2 o e 3 0 e n 0 

directed at the accelerator grid (plane "C"). Because of the retarding 
field seen by the electrons as they enter the region between the planes their 
velocities will drop to zero at planes P , P 2 . . . . P n which are at poten- 
tials V , V ... V respectively. The electrons will, after they come to 
1 2 n 

rest, be accelerated back toward plane "A". Again it is assumed that the ions 

and electrons do not collide with each other or with neutral particles between 

the grids. As a result of electron injection at current densities j , j . . 

e l e 2 

j corresponding to injection velocities v , v ... .v the ion 

i 0 \ 

current density will reach a level j.. The ion charge density anywhere 
between the planes "A" and "C" is given by Eqs. (2) and (3) as in the 
preceding cases discussed. The electron charge densities will, however, 
depend on which region is being considered. In the region between planes "A" 
and "P ", for example the total charge density is the sum of charge densities 

p e ’ p e * ' * ’ P e„ associated with each of the incident beams as well as 
l 2 n 

the reflected ones that pass back through this region. 
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In the preceding expressions V is the potential at a location x between 
planes "A" and "P" . Poission's Equation in this region becomes: 



nondimensional ization using the parameters defined in Eqs. (1) and (7) 
enables one to rewrite Eq. (34) in the form 
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(35) 


In this equation j. is the ion current density without electron injection 
'o 

and 


a k = 



(36) 


Integrating Eq. (35) and applying the space charge limited boundary condition 
^ = 0 at A = 0 one obtains 
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Integration of this equation using the boundary condition <j> = 1 at A = 0 
gives the normalized location (A ) of the plane P in terms of the potential <j> 1 


A = 


/J 1 Vs 1 

- n n - 


^ 2a k ((f)-<fi k ) 2 + (1-0) 2 - 2a k (l-0 k ) 2 

1 

\l 

II 


d<f> 


(38) 


In other regions groups of injected electrons that do not reach those regions 
are not included in the summations but the procedure followed above to 
obtain potential gradients and plane locations is similar. In the general 
region between planes P- , and P. for example one obtains 

J * J 


d£ _ if il_ 

dA 3 l J 1b 


2a [< ((j)-(}i| < )' 2 + (l-<j O' 5 - 2«| < (l-<j)| < ) 2 

k=j k=l 
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and 


j, \% 


r + 


A • = A.- t + f( — -° 
J J-l 4 \ j. 


j-1 n 


/ J- i II 

+ (l-<j>) % - J 2cx k (l-4, k ) !s 
k=j k=l 


d<j) (40) 


The preceding expressions are valid for j = 2, 3 . . . . n. For the 
last region to the right in Fig. 8 (Region P n -C) one obtains 
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which can be integrated to give 



(41) 


(42) 


Adding the equations of the form of Eqs. (38), (40) and (42) describing 

the locations of the planes "P ", "P "... ."P " one obtains the normalized 

12 n 

ion current density ratio: 
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Applying the space charge limited current condition at each of the planes 

"P ", "P " "P_" to Eqs. (37), (39) and (41) one obtains the following 

12 n 

set of simultaneous equations 
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Simultaneous solution of Eqs. (44) was accomplished using Cramer's Rule 

to obtain expressions for the constants a . . „ . a in terms of the 

l n 

normalized potentials <f> . „ . . <j> n . The previous cases considered in this 

paper suggested that the maximum ion current density was realized when the 

normalized potentials of the planes "P " and "P " were <j> = 0.4 and 

1 2 1 
2 9 

<J) = $ = 0.16. Using the simplex method for cases of two and three 

2 1 

injected electron groups this result was verified for the two group 

problem and the potentials for the three group system were found to be 
2 3 

<j>1 = 0.4, <j> 2 = (j>i = 0.16 and ^3 = = 0.064. The logical extension of 

this is to establish injected electron energies such that the generalized 
potential is 

<f> k = (0.4) k . (45) 

Numerical integration of Eq„ (43) carried out using values of 
determined using Eq. (45) together with the solutions of Eq. (40) suggested 
that this criterion does indeed result in the maximum normalized ion current 
density for the general case. Normalized ion current densities have been 
computed for one through ten stages of electron injection and the results 
of this work are shown in Fig. 9. This figure shows that the ion current 
that can be extracted from a plasma increases monotonically with the number 
of stages of electron injection to a value over thirty times the uninjected 
Child's Law limit. This value of extracted ion current is realized when 
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Fig. 9 Effect of Number of Optimally Injected Electron Beams 

¥ 


on Extracted Ion Current 
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ten stages of injection are used. The normalized electron current density 
required for each injection stage at this optimized condition as a function 
of the number of injection stages is shown in Fig. 10. The greatest 
electron current densities are seen to be required for the first (lowest 
energy) stage regardless of the number injection stages used. The total 
electron current density required would be the sum of all of the values 
given by the curves that apply for the number of stages being used. For 
the cases of upstream electron injection investigated this value lies between 
0.5 and 0.6 times the ion current and the square root of the extracted ion 
to injected electron mass ratio. 

As an example consider the case of electron injection from the upstream 
side of the screen grid at two energies. In this case the potential, 
potential gradient and ion velocity profiles at the maximum extracted ion 
current density operating condition are exactly the same as those shown in 
Fig. 7 for the case of two injected electron groups, one injected from 
upstream of the grids and the other from between the grids. For the example 
where upstream injection of both beams is used, the maximum extracted ion 
current density is observed for 


“1 

a 2 

*1 

*2 


0.42 

0.15 

0.4 
2 




1 


from Fig. 10 

from Eq. 45 

0.16 ' 


At this condition the normalized ion current (j./j,- ) from Fig. 9 is 11.7 

1 ‘o 

and the total injected electron current has the same value as previously determined. 
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Fig. 10 Injected Electron Current Density Required to Yield the 
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Comparison of these results with those obtained previously for upstream 
injection of the primary electron beam (i.e. at <f> ) and intra-grid injection 
of the secondary one (i.e. at <}> 2 ) shows that they are the same except for 
the fact that is 0.15 rather than 0.30 in the present case. This is 
simply a reflection of the fact that only half of the injected electron 
current density is required when upstream injection is used and the electrons 
pass through the space charge region twice rather than a single time. 

It is appropriate to ask whether electron injection at a fraction of 
the optimum values given in Fig. 10 might not give substantial increases 
in extracted ion currents (i.e. a substantial fraction of the values shown in 
Fig. 9). The results shown in Fig. 11 suggest the answer to this question. 

In this figure the normalized extracted ion current is plotted against the normalized 
first stage of injected electron current for a typical electron injection 
case (three stages of injection with the second and third stages at optimal 
current densities). Unfortunately the figure shows that the ion current does 
not begin to increase substantially until the normalized electron current 
approaches its optimum value (0.4 in this case). The figure also suggests 
one should just keep increasing the electron current above 0.4 to obtain 
still higher ion currents. It is however physically impossible to increase 
the electron current above the optimum value without inducing changes that 
alter the potential distribution between the grids in such a way that the 
normalized ion current decreases. Mathematically this is manifested in 
the equations by solutions for injected electrons current density that 
become imagainary. 

Itshouldbe pointed out that whil e the analysis presented here has assumed 
the injected negative particle currents would be produced using electrons, 
the analysis is valid for negative ion injection. In this case the m g 
appearing in the development would be the injected negative ions. 
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Concl usions 

The current density of a beam of charged particles being accelerated 
between two parallel plate electrodes can be increased substantially by 
injecting oppositely charged particles into the region between the electrodes. 
Multiple beams of oppositely charged particles can be used to produce the 
most substantial increases in the current density of the primary beam. The 
current density of this beam increases monotonically with the number of 
optimally selected, injected particle beams being used. There is an optimum 
injected particle current density and energy associated with each beam that 
yields the maximum primary beam current density for a given number of in- 
jection stages. It is perferable to inject the neutralizing particles from 
outside of the acceleration region rather than from within this region to 
minimize the required injected particle current densities. However either 
release of these particles in the acceleration region or injection with 
sufficient energy to carry them against an adverse electric field to the 
proper intra-grid potential produces the same solutions. The total required 
injected particle current density, regardless of the number of stages used, 
is of the order of the ion current density times the square root of the ex- 
tracted ion-to-injected particle mass ratio. The primary beam current 
density that can be achieved decreases rapidly below its maximum value with 
small decreases in the injected particle current density below the optimum 
value. The fact that the extracted ion-to-injected particle mass ratio is 
high if electrons are the injected particle means that the injected particle 
electron current density would have to be very high for typical ion current 
densities. It is because they would have to be so high that the technique 
would probably be impractical for use in space propulsion applications. 
Negative ions could be injected and the required injected particle current 
density would decrease substantially. The required negative ion current 
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densities would however be comparable to the positive ion current densities 
produced in the beam for this case. Efficient production of such a beam 
however poses the same problem as generation of the positive ion beam that 
would be the objective of the final device. Recent experiments ^ 0 suggest 
that these concepts are of interest however in ground based applications 
where power demand is a lesser concern. 



THE ELECTROTHERMAL RAMJECT 


Ben D. Shaw 

In order to extend the practical uses of space for mankind it is de- 
sirable (1) to reduce the cost for delivery of a payload into space from 
earth, (2) to simplify the systems needed to make these launches and (3) 
to improve their reliability,, One system that holds some promise of 
achieving these objectives involves the acceleration of a payload to a 
sufficiently high velocity at the earth’s surface so it would be able to 
undergo passage through the earth's atmosphere and still have adequate 
velocity to enable it to meet its mission objectives in space. One candi- 
date propulsion system that has been analyzed and appears to hold promise 
for use in such a launching scheme is the electrothermal ramjet.^ It 
appears to be capable of operation at the high thrust levels needed to 
achieve the high launch velocities desired with a reasonably short launch 
tube. 

The concept of this launch scheme is illustrated in Fig. 12. To effect 
the launch, propellant, preferrably a low molecular weight gas, would be dis- 
tributed in the launch tube upstream of the projectile shown. This projec- 
tile, which is also the payload, would be brought to an initial velocity 
where efficient ramjet operation could begin by driving it with a light gas 
gun driver.^ At this point the projectile shown proceeding down the launch 
tube in Fig. 12 would begin to process the propellant in the launch tube 
through its annular intake diffuser. The diffused propellant would then 
pass through a heat addition zone in the ramjet where energy would be added to 
raise its stagnation temperature. Heat addition would be effected using either 
electromagnetic energy radiated from pulsed lasers or directly deposited 
energy from pulsed electrical discharges located along the tube wall. 
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These devices would be triggered sequentially as the projectile passed down 
the launch tube. The heated propellant would then pass through the nozzle 
and would be expelled at a sufficiently high pressure and velocity to 
produce the desired thrust. The thrust level desired would actually be 
attained by adjusting the initial temperature and pressure of the propellant 
in the tube, the rate of heat addition to the propellant, and the ramjet 
geometry. 

A preliminary analysis of this concept was conduced during the grant 
period and the results of this work were presented in the paper 11 indued in 
Appendix A. In this work it was assumed that the diffusion process illustrated 
in Fig. 12 could be model led as an isentropic one. It was recognized that 
this diffuser model was idealized and that the flow through the actual dif- 
fuser would encounter a system of shocks. The preliminary analysis based on 
the isentropic assumption suggested that the concept was workable and that 
reasonable energy utilization efficiences could be expected with it. 

Theoretical Model of Diffuser Shock Losses 

It was recognized that the diffusion process within an electrothermal 
plug ramjet (Fig. 12) would be non-isentropic because of shock losses. For 
this reason the first extension of the preliminary thermodynamic analysis 
of the electrothermal ramjet has been the incorporation of these losses. 

For supersonic operation a conical shock will either be attached to the nose 
of the projectile or a bow shock will be driven ahead of it. Which will 
occur depends on the velocity of the projectile and the geometric and boundary 
conditions associated with the flow. Of the two options the attached conical 
shock is preferrable because the stagnation pressure loss across a conical 
shock system is less than that across a bow shock seeing the same free stream 
flow. Consequently, flow conditions that enhance attachment of a conical 
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shock are desirable and a conical shock model has therefore been used in 
the analysis. 

In order to incorporate shock losses into the complete analysis of the 
device, the flow conditions immediately downstream of the conical shock 
shown in Fig. 12 were first modelled. In an actual diffuser of the type 
shown, the attached conical shock would be followed by complex series of 
reflected shocks terminated by a normal shock. The location of this termi- 
nating normal shock would depend on the geometry of the device and the flow 
boundary conditions. Ideally the stagnation pressure loss associated with 
the initial conical shock and the series of reflected shocks will be less 
than that of a single conical shock terminated by a normal shock. However, 
a review of the literature on hypersonic wind tunnel diffusers (Ref. 12 for 
example) has shown that the actual stagnation pressure recovery that can be 
expected is very close to the stagnation pressure recovery prediced theo- 
retically by assuming that all of the flow goes through the single conical 
shock followed by a normal shock. Because of this, the conical /normal shock 
model suggested in Fig. 12 was used in the analysis. 

The thermodynamic properties of the working fluid will vary over the 
range of temperatures and pressures encountered during ramject operation 
primarily because of the effects of dissociation. For this analysis, however, 
it was desirable to avoid the complications of temperature and pressure de- 
pendent specific heats and gas constants. In addition, increases in propel- 
lant temperature predicted during ramjet operation are generally accompanied 
by corresponding increases in pressure. This being the case, it is not too 
unrealistic to use fixed values for the specific heats and the gas constant. 
The results contained in this analysis were obtained assuming ideal gas 
behavior and diatomic hydrogen was the working fluid. A gas constant (r) 
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of 4160 joules/kg°K and a ratio of specific heats (y) of 1.4 were used. It 
is also noteworthy that the exhaust pressure for the device turns out to be 
quite high for the cases examined to date. This suggests that recombination 
rates in the exhaust will be high, that frozen flow losses will be small and 
this should also tend to make efficiencies computed using the ideal gas 
assumption reasonably accurate. 

The model suggested in Fig. 12 implies an axi-symmetric (two-dimensional) 
conical flow field immediately downstream of the conical shock that must 
somehow be matched to the one-dimensional model used to describe the flow 
downstream of the normal shock. The marrying of these flow fields necessitated 
the development of an averaging procedure that provided single values of flow 
properties at the conical flow field exit for input into the one-dimensional 
normal shock relations. This was accomplished by determining a mass weighted 
average of each property using the mass, energy, and momentum conservation 
equations. Downstream of the postulated normal shock in the diffuser, all 
flow throughout the ramject was assumed to be well mixed and one-dimensional 
in nature. The flow exiting the diffuser passes through a constant area, 
frictionless passage where heat addition at a subsonic Mach number is assumed 
to occur. Next the flow passes through the nozzle where isentropic expansion 
to the maximum attainable supersonic Mach number is achieved under the con- 
straint of the available exhaust area. 

It is also assumed that the ramjet exhaust pressure can be maintained 
at that value dictated by the isentropic expansion through the nozzle. This 
static exhaust pressure required at the projectile base will vary with pro- 
jectile position and control of its magnitude would have to be achievedinan 
actual device by proper design of the launch tube. Such a tube would consist 
of two sections, a cavity section into which propellant processed by the 
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ramjet would be collected and the launch section through which the pro- 
jectile would pass. It has been assumed in this analysis that the cavity 
section would be designed to have the length and cross-sectional area needed 
to accommodate the propellant and to facilitate maintenance of the necessary 
projectile base pressure. While it seems possible to do this intuitively, 
the unsteady fluid mechanical analysis of this process has not yet been 
attemped. 

The basic equations of fluid mechanics (conservation of mass and energy 
coupled with a thermodynamic path equation) were applied to all processes 
occurring within the ramjet. Solutions were obtained numerically. Propellant 
velocities were computed in the frame of reference of the moving projectile. 
The fact that the tube wall, which is the outer boundary to the flow, is 
moving in this frame of reference was assumed to have no effect on the flow 
behavior. This analysis assumes that the projectile has been accelerated to 
a high enough velocity before the initiation of ramject operation so that 
attachment of the conical shock to the diffuser nose is assured. Thrust was 
computed by applying the momentum equation in the frame of reference of the 
control volume moving with the projectile. 

It should be noted that the boundary conditions associated with the 
annular flow ramjet are different than those associated with conventional 
ramjet operation. In the annular flow ramjet the inlet and outlet flow areas 
must be equal and the flow area in the heat addition zone must be less than 
this area. Further the ramjet inlet and outlet pressures will be different. 
For the conventional ramjet on the other hand the inlet, outlet and heat 
addition region flow areas are not constrained, but the inlet and outlet 
pressures must be equal . 
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Results for a Typical Supersonic Launch Cycle 

In the analysis procedure that has been described it is assumed that 
the projectile will be acted upon by either the ramjet thrust or some ex- 
ternally applied force needed to produce the desired acceleration. As a 
result of this assumption the projectile trajectory is known. This in turn 
implies the projectile velocity is known and hence the propellant flow velocity 
into the ramjet diffuser is also known at each position along the launch 
track. In order to minimize the mechanical complexity of the projectile 
it has also been assumed that the projectile geometry remains fixed through- 
out the launch cycle. The diameter of the launch tube is allowed to vary 
as a function of axial position however so efficient ramjet operation can 
be achieved. Further, the present analysis requires that the flow through 
the ramjet is choked at the nozzle throat. While choking could occur at 
the diffuser throat rather than the nozzle throat, these two throat areas 
were selected so that choking would occur in the nozzle because this gave 
the best performance. Using the analysis one determines how the tube diam- 
eter and heat addition rate should vary so that a prescribed thrust can be 
maintained throughout the period of ramjet operation. 

A typical case of electrothermal ramjet operation has been examined. 

In this case hydrogen propellant is placed in the launch tube at a pressure 
of 30 atmospheres and a temperature of 300°K. A 10 kg projectile having a 
16 cm diffuser throat diameter (Fig. 12), an 8 cm center body diameter and 
a 16 cm nozzle throat diameter is used. These values were selected because 
they gave reasonable performance and were at the same time sufficient to 
ensure structural integrity during launch. It is also required that the 
thrust be sufficient to sustain the projectile at a constant acceleration 
level of 30,000 g until it reaches a final velocity of 15,000 m/sec. 
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The results of the analysis in which diffuser shock losses are in- 
cluded are compared to results obtained using an isentropic diffuser model 
in Fig. 13 for this particular case. Figure 13a shows how the projectile 
and exhaust velocities vary during the launch cycle. The projectile velocity 
(relative to the stationary launch tube) is observed to increase linearly 
with time as a consequence of the constant thrust and constant projectile 
mass that have been specified. The exhaust velocities (relative to the 
projectile) computed for the isentropic and non-isentropic cases are shown 
in Fig. 13a for times greater than 10 to 16 msec. Prior to these times, 
the conical shock would not attach, adequate thrust could not be produced 
regardless of the level of heat addition to the propellant, or the tube 
diameter became so large that the assumption of a one-dimensional flow model 
was seriously violated. During these initial phases a gun mode of operation 
would be required to accelerate the projectile to the velocities at which 
proper ramjet operation could be initiated. The parameter appearing in 
this figure (6) is the half angle of the diffuser cone. A half angle of 
about 22.5° is observed to produce exhaust velocities that lie closest to 
the isentropic ones. At half angles greater or less than this value (5° and 
45° are shown) the exhaust velocities are less at all times. The data of 
Fig. 13a show that •v 50 msec would be required to accelerate the projectile 
to a final velocity of 15 km/sec. This would require a launch tube about 
400 m long. 

Figure 13b indicates how the launch tube diameter would have to vary 
with time in order to maintain unity Mach number at the nozzle throat for 
these same isentropic and non-isentropic cases. It is important to hold the 
nozzle throat Mach number at unity to get the best performance from the 
device. If the tube diameter is less than the value indicated in the figure 
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Fig. 13 cont. Typical Ramjet Launch Profile for a 
10 kg Mass Accelerated at 30,000 g 
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at a given time and all other input parameters are fixed, then the conical 
shock will detach and losses will increase. If the tube diameter is greater 
than this value the heat addition will be insufficient to maintain unity 
Mach number at the nozzle throat, the nozzle flow will go subsonic and poor 
ramjet efficiency will result. In the case of Fig. 13b one again observes 
that a diffuser cone half angle of 22.5° produces results that are closest 
to the isentropic ones. 

Figure 13c illustrates the pressures required at the projectile base 
for the examples of ramjet operation being considered. The base pressures 
for the non-isentropic cases rise as the projectile travels down the launch 
tube. The non-isentropic cases are observed to exhibit higher exhaust pres- 
sure requirements over most of the launch sequence. Once again a diffuser 
cone angle of 22.5° yields conditions that approach the isentropic ones most 
nearly. It is of interest to note that at the beginning of ramjet operation, 
for each of the non-isentropic cases considered here the high projectile 
base pressure is the source of all of the required thrust, and the momentum 
change of the flowing propellant across the ramjet actually detracts from the 
thrust. As ramjet operation continues, the momentum change of the propellant 
begins to dominate the thrust, although not completely. For example, at the 
end of the operating period for the 22.5° cone angle device, = 35% of the 
thrust is being provided by the pressure differential across the projectile, 
and z 65% of the thrust is being provided by the momentum change of the 
propel lant. 

Figure 13d shows the temporal variation in propellant exhaust temperature 
that is induced by the heat addition profile shown in Fig. 13e. The tem- 
peratures and thermal input powers shown are high, but not excessive in 
light of the short launch times involved. The power going into the pro- 
jectile is also shown in Fig. 13e and the ratio of this power to the thermal 
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input power at a given time is termed the instantaneous efficiency. This 
ratio is plotted as a function of time for the cases considered in 
Fig. 1 3f . 

This figure shows that for the case under consideration the non- 
isentropic diffuser having a cone half angle of 22.5° again produces results 
closest to the isentropic ones. This suggests then that there is an optimum 
diffuser shape and that the optimum conical diffuser has a cone half angle 
of 'v 22.5°. The isentropic efficiency is observed in Fig. 13f to be 
substantially greater than that for the best non-isentropic case (6 = 22.5°) 
over most of duration of the launch. The efficiency of the electrothermal 
ramjet equipped with this optimized realistic diffuser is however still very 
attractive. It exhibits an efficiency near the end of the launch of about 
50% and the mean efficiency over the complete cycle beginning at 10 msec is 
observed to be greater than 40%. This efficiency is much greater than the 
efficiencies of chemical and light gas gun launchers and is competitive with 
efficiencies predicted for rail guns and mass drivers. 

The reason why an optimum diffuser cone half angle is observed in the 
data of Fig. 13f can be understood by considering the fluid flow conditions 
associated with the limiting cone half angles (i.e. values of 0° and 90°). 
For the present model we are requiring subsonic heat addition so the flow in 
the diffuser section must pass through a normal shock wave at some location. 
For a specified incoming (free-stream) Mach number, each of these cones 
coupled with the diffuser, will produce identical stagnation pressure re- 
covery effects. The very shallow angle 0°) cone will produce a Mach wave 
at the diffuser entrance across which the Mach number of the flow remains 
constant. The flow will then enter the supersonic diffuser, where it must 
encounter a normal shock in order to reach the subsonic flows required in 
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the heat addition region. This means the stagnation pressure loss across 
the shock will correspond to the free-stream Mach number. The blunt cone 
will propagate a bow shock, which is essentially a normal shock, ahead of it. 
The bow shock will cause the flow entering the diffuser to become subsonic, 
and the stagnation pressure loss will again correspond to the free stream 
Mach number. This analysis has shown that a moderate cone half-angle 
(between 0° and 90°) and the corresponding conical shock coupled with a 
normal shock will produce more efficient stagnation pressure recovery than 
will either of these single normal shocks. 

Conclusions 

There is an optimum cone angle for the diffuser in an electrothermal 
ramjet which yields maximum energy efficiency. For a typical case con- 
sidered here it is ^ 22.5°. The efficiency of the electrothermal ramjet 
is degraded from the isentropic efficiency when account is taken for shock 
losses occurring in the diffuser. Overall energy efficiencies for a typical 
mission, where acceleration to 15,000 m/sec is effected, still lie in the 
very attractive range of 40 to 50% when these losses are considered. 



CONSTANT PRESSURE ACCELERATION OF THE ELECTROTHERMAL RAMJET 


Charles E. Mitchell 


Introduction 

No ramjet engine can operate efficiently until it has been accelerated 
to some minimum operational velocity. For the electrothermal ramjet system 
discussed in Reference 11 and elsewhere in this report this minimum opera- 
tional speed is quite high. Preliminary calculations performed for a 
vehicle mass (m^) of 10 kg, an acceleration rate (a) of 30,000 g, and 
employing diatomic hydrogen as the propellant indicated that this minimum 
speed could be as high as 5000 m/sec (Ref. 11). Although certain designs 
facilitate ramjet operation at considerably lower speeds it is assumed con- 
servatively in the following analysis that the starting system would have 
to bring the ramjet to 5000 m/sec. 

One method of obtaining speeds of this order in short time periods at 
acceleration levels of the order of those assumed for the electrothermal 
ramjet is acceleration in a light gas gun system. Devices of this type 
have been both analysed and constructed and are discussed in the literature 
(see for example Ref. 13). In what follows two acceleration systems based 
on the light gas gun concept will be described and analysed. Both systems 
are designed to produce a high level of constant acceleration (and therefore 
constant base pressure) for the ramjet. The mode of energy addition re- 
quired to produce the acceleration is taken to be similar to that assumed 
for the ramjet mode of operation; that is, electrothermal in nature and 
controllable in time. 

A general description of the acceleration process will be presented 
first. This will be followed by the analysis of an isentropic flow, piston 
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dr iven dual chamber device. This type of system is typical of previous 
light gas gun designs. Next a non-isentropic flow, single chamber design 
will be described and analysed. This latter device differs from conventional 
designs but appears attractive in the current application. Results of 
calculations for a specific example will be presented and the relative 
performance of the two designs compared. 

Uniform Acceleration Flows - General Description 

If a constant pressure is maintained at the rear surface of the ramjet 
then, in the absence of frictional effects, a constant acceleration of the 
ramjet vehicle will occur. (It can be pointed out in passing that available 
data (see Reference 14) indicate that frictional effects between the 
projectile (ramjet) and the launch tube should cause a decrease in velocity 
obtained of about 5%.) The basic problem is to devise a method by which a 
constant ramjet base pressure can be maintained. In Reference 11 a quasi- 
steady expansion was assumed so that all the gas in back of the ramjet 
(in the launch tube) was at a constant pressure during acceleration. The 
high pressure level necessary was to be maintained by heat addition. Since 
sound speeds of the order of 3000 m/sec and vehicle speeds up to 5000 m/sec 
are desired, it is clear that the quasi-steady approximation cannot be correct, 
and that the flow behind the ramjet is inherently unsteady. 

One unsteady flow which does provide constant pressure at the ramjet 
is the flow which results if all fluid particles in back of the ramjet are 
accelerated uniformly. That is, the gas in back of the ramjet and the 
ramjet itself are accelerated as a solid body. Conceptually, this is the 
flow situation which would result if the ramjet and all the gas behind it 
began accelerating at the constant value a starting at some arbitrary time, 
say t = 0. The determination of the variation of the dynamic and thermo- 
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dynamic variables for this flow as a function of axial location x and 
time t is relatively straightforward. 

Consider first the pressure-distance condition which must exist at 
t = 0. Since solid body acceleration is occurring, then for the gas in back 
of the ramjet it must be true that 


d7 Q =-P« (46) 

where P is local pressure, p local density and x Q the axial distance along 
the tube at t = 0. This is of course, recognized as the usual gravitational 
equation with g replaced by a. If a property relationship between P and p 
is given (say isentropic or isothermal for example) then P as a function of 
x Q along the tube axis could be found at t = 0. For later times the way 
in which the pressure must increase at any x location is determined by the 
uniform acceleration condition itself. For any fluid particle (and for the 
ramjet) 


x = 


x 


o 


+ 


at 2 

2 


(47) 


where x is the instantaneous location of the particle and x q is its location 
at t = 0. Particle paths on an x, t diagram for the case where x q of the 
ramjet is zero, are shown in Fig. 14. Also, since the acceleration is 
uniform, the state of any fluid particle cannot change in time. That is 


DP 

Dt 


= 0 , 





(48) 


where is the material derviative and T and s are temperature and entropy, 
respectively. From Eq. (47) the identity (x q value) of the fluid particle 
which passes location x at time t is immediately found as 


x 


o 


x 


at 2 

2 


(49) 
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a = 30,000 g 
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Fig. 14 Uniform Acceleration: Time vs. Distance 
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Fig. 15 Accelerator Design 
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Then, also, 

dx 

st~ = - at (50) 

for any constant x location. 

Combining Eqs„ (50) and (46) gives 

a 2 pt (51) 

at any x station along the tube. Also at any x station 

u = ext . (52) 

In order to integrate Eq. (51) a statement which relates P to p along 
the line t = 0, and permits the integration of Eq. (46) must be given, 
that is 

P 0 = P 0 (p 0 > • (53) 

Since the pressure and density of every fluid particle remain constant for 
t >_ 0, then Eq. (53) is valid for all times 

P = P(p) . (54) 

Integration of Eq. (51) using Eq. (54), and application of Eq. (52) permits 
the specification of all thermodynamic and dynamic variables at any x 
location in the uniformly accelerating flow as a function of t. 

The practical importance of such a specification to the ramjet accel- 
eration problem can be seen as follows. If the ramjet is located at x = 0 
at t = 0 and if P, u, p etc., take the values required by Eqs. (51) - (54) 
along x = 0 then the flow field for x >_ 0, t >_ 0 will be a uniform accelera- 
tion flow field. This means that if a flow with the required physical param- 
eters at x = 0 can be suppl ied to the tube at x = 0 starting at t = 0 then a 


constant pressure acceleration of the ramjet wil 1 occur. In practice this could be 
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achieved through the use of a large diameter driving chamber attached to 
the launch tube at x = 0. Appropriate variation of pressure and density in 
the driving chamber (by piston motion or heat addition for example) would 
then give the required values of P, p, and u at x = 0. Two specific designs 
for providing the required parameters at x = 0 have been considered,, The 
first of these, a dual chamber piston driven accelerator, will be discussed 
next. 


Piston Driven Isentropic Accelerator 

A sketch of this device is shown in Fig. 15. The general acceleration 
sequence would occur as follows. First, with the ramjet held at x = 0 the 
pressure in the compression chamber would be raised to the value required 
for the desired acceleration level. For an acceleration of a, this pressure 
would be 


m a 

P = P = _i 

C 00 A-r 


(55) 


where P Q0 is the pressure at x = 0 and t = 0, is the ramjet mass, Ay is 
the launch tube cross sectional area and P c is the compression chamber pres- 
sure. The ramjet would then be released and the solution of Eqs. (51) to 
(54) would specify the required conditions at x = 0. These would be met 
through appropriate motion of the piston which would in turn be effected by 
heat addition in the heating chamber. This would continue until the ramjet 
reached its operational speed. A complicating factor is the possibility of 
choking the flow at x = 0. This will be discussed below. 

Since the heat addition occurs in gas which is physically separated 
from the gas which will enter the launch tube, it is reasonable to assume 
that all gas entering the tube has the same entropy value and that the 
flow field in both the compression chamber and launch tube is homentropic 
(constant entropy everywhere). If in addition the gas is taken to be 
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perfect (ideal gas with constant specific heats) then the conditions at 
x = 0 can be expressed in closed form as functions of time, as can the 
heat addition rate in the heating chamber. It will also be assumed that 
the diameter of the compression chamber is large compared with that of the 
launch tube so that the properties in the compression chamber are constant 
in space except for a negligibly thin acceleration region which connects 
the chamber to the tube and which can be taken to the quasi-steady. Finally, 
the mass of the piston and the friction associated with its motion are 
neglected. 

In solving the problem thus posed it is convenient to define the 
following dimensionless variables 


P = P/P 

oo 


X = 



u = 



P = p/p 00 , T = T/Tqq » 


Q H - _ aL 

H ~ Pa A ’ L 32 » 

00 00 C 00 


Qu = 


a = 


oo 



where the subscript oo indicates a quantity evaluated at x 
a„ = /yRT „ is the reference state sound speed and y and 

oo 00 r i 

of specific heats and gas constant, respectively. 

The isentropic condition then gives 

P = p Y 


= 0 and t = 0, 

R are the ratio 


(56) 


which is a special case of Eq. (54). Using this relationship and non- 
dimensional izing Eq. (51) becomes 

1 

= y ( "F ] Y t . (x = 0) (57) 

dt \ / 

Integrating, and using the ideal gas law, P = p T, gives at x = 0 
(superscript o) 
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T° = [l + izjjU t 2 ] 

5° = [i + 



1 

Y-l 


?> = 1 + P 


|0 _ U_ _ t 

a [(y-DF/ 2 + l]* 5 


(Mach Number). 


(58a) 

(58b) 

(58c) 

(58d) 


Equation (58d) 
The result is 


can be solved for t*, the dimensionless time at which M° = 1. 



(59) 


Clearly, the solution as outlined above is only valid to this point. 

However, it must take some time for the information that the flow has choked 
to reach the ramjet. That is, there must be a period of time after t = t* 
during which the ramjet will continue to accelerate at a with a base 
pressure P Q0 . The exact value of the time involved is determined by the 
intersection of the right running characteristic emanating from x = 0, 

¥ = [3“]^ > with the ramjet path on the >T, t plane. Without giving the 
details (see Ref. 14 for example) the resulting intersection time is 




max 


(y-1) 


(60) 


This is the maximum time that the ramjet can experience uniform acceleration 
and constant base pressure, under the assumptions made here. 

If the gas in the tube and compression chamber is assumed to be hydrogen, 
some numerical calculations can be made to assess the importance of the 
maximum acceleration time condition. For hydrogen, y = 1.4 and thus t* = 1.118 


and Uv = 1-708. Since u m3V = t . and u m:w = u m=v a„. an estimate 
max max max max max 00 
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of the T qo necessary to reach a certain u max in time t max can be made„ 

T u m ax 1 2 

For hydrogen this relationship is T Q0 = I yjg - 3I • If u max is the conser- 
vatively high value of 5000 m/sec then T Q0 = 1472°K. This temperature does 
not seem unreasonably high. Moreover, the maximum temperature which would 
occur would be at t = t* when T max = 1840°K. If it is desirable or necessary 
to operate at lower temperatures and still reach 5000 m/sec, it is possible 
to accelerate the ramjet in a non-uniform fashion after t = t m:w . The nature 
of the flow field in the tube after tf = t* is, of course, non-uniform and 
must be determined using the method of characteristics. Solutions indicate 
nearly uniform acceleration at the ramjet can still be achieved by appropriate 
tailoring at x = 0 (see Ref. 15). For the purpose of the present analysis, 
the temperature condition will be taken as not overly restrictive and the 

t of interest will be taken to be t < ¥ „ 

— max. 

The values of P", p" and T as given by Eqs. (58a, b, c) lead through the 
quasi-steady assumption and the isentropic perfect gas relationships, to 
the following expressions for T , p , P , the non-dimensional temperature, 

V L 

density and pressure in the compression chamber 

1 + (y- 1) t 1 2 (61a) 

1 

[l + (y-D t 2 ] Y_1 (61b) 

1 _L_ 

[l + (y-1 ) P] H . (61c) 

The ideal gas law written for the compression chamber is 

P„ ? A. - R T„ (62) 

c c c c 

where m c is the mass of hydrogen in the compression chamber at any instant 
and 5 is the instantaneous length of the compression chamber (see Fig. 15). 
Also the quasi -steady outflow at x = 0 implies that 
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dm 


1 


c _ 


dt - -W - [l +-=lt2] Y - 1 t . (63) 

Combination of Eqs. (62) and (63) leads after some manipulation to 
the following expression for c , the non-dimensional compression chamber 
length 


!_ = 


1 - 0 


[l + (y - 1 ) t 2 ] 


1 

Y~ 


(64) 


Z Q is the compression chamber length at t = 0. It is related to the final 

clearance length in the compression chamber which exists at t = t*. For 

Ay i _ _ _ _ 

example, if y = 1.4, = yg- , and ? -> 0 at t = t* , then = 0„0528. This 

corresponds to 5 = 1 .540 m if u is 5000 m/sec. The quantity is the 

lllaA ■ ■ 


dt 


non-dimensional piston speed. For the example just discussed 

|~~| takes a maximum value of 0.0622 which corresponds to a piston speed of 
dt 

about 183 m/sec for u max = 5000 m/sec. The piston speed is thus quite small 
compared with the sound speed in the chamber 3000 m/sec and the maximum 
ramjet velocity 'v 5000 m/sec). This suggests that the assumption of uniform 
pressure, temperature and density in the compression chamber is probably not 
unreasonable. 

Several methods for producing piston motion of the type required by 
Eq. (64) have been suggested in the literature (Ref. 15) „ In the present 
analysis it will be assumed that energy can be added to the heating chamber 
behind the piston by electric discharge or laser addition methods similar 
to those envisioned for the electrothermal ramjet itself. In particular it 
will be assumed that the rate of energy addition can be controlled closely 
enough to meet the requirements of the piston motion presented above. The 
necessary heat addition profile can be determined from application of the 
first law of thermodynamics to the heating chamber. This is 
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dU, 




W„ = 


H 


H dt 


(65) 


where Q H is the required heat addition rate, is the power necessary to 

du H 

move the piston as required by Eq. (64), and is the rate of energy 
storage in the gas in the heating chamber . 

If the piston is massless and frictionless (these effects can be 
shown to be relatively small) then 


u = _p dJL a = _ p !■!£. a 
W H k H dt H C K C dt M C 


( 66 ) 


with determined from Eq. (64), while 


m r 5i 

dt m H L v dt 


(67) 


where m^ is the (constant) mass of gas in the heating chamber and C v is the 
specific heat at constant volume for the same gas. Application of the 
ideal gas law and substitution of Eqs. (66) and (67) into Eq. (65) gives 
after some rearrangement 


OO 00 C 


_ _Y_ 

" Y-l 


, fl J_ _L , 

|^[l + 2 L[h(y- 1) F] Y '' jt (68) 


where L = L -f- and L is the combined length of the compression and heating 
oo 

chambers (see Fig. 15). The value of L is somewhat arbitrary, though it does 
have an impact on the temperature rise in the heating chamber. This is 
readily seen upon application of the ideal gas law to the gas in the heating 
chamber. This gives 



( 69 ) 


If 5 = 0 at t = t* = 1 .118 then 
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T * 

= 4.133 . (70) 

'H l " 5 o 

It is clear from either Eqs. (69) or (70) that small values of L lead to 
large temperature rises. The maximum value of T^, T^* still depends on 
T^ q . Since T Hq need not be equal to T a separate control on the maxi- 
mum temperature is possible beyond that implied by the choice of L. In 
the numerical results to be presented shortly, L was taken quite arbitrarily 
to be 25 q . This means that at t = 0 the driving chamber is separated into 
two equal volumes by the piston. Other choices of L would lead to higher 
or lower values of T^ for the same T^ Q . As will be pointed out later 
the choice of L also affects the performance of the device as far as energy 
conversion efficiency is concerned. 

The efficiency of the accelerating process at any instant of time can 
be measured by the ratio of the rate of change of the kinetic energy of the 
ramjet vehicle to the rate of heat addition to the heating chamber Q^. The 
rate of change of the kinetic energy of the ramjet is 


d(m u 2 /2) 


dt 


= u4r- = ni a 2 t = m. aa t = PA. at . 
r dt r t oo oo t oo 


Thus, using Equation (68) 
d(m u 2 /2) 


n = 


dt 


<H 


111 

Y 


[l +*j^P] 7-1 + 2L ^ [H( Y .l )P]tH 


-1 


(71) 


The dependence of n on t for y = 1.4, = 16, L = 2 e, q = 0.1056 (L = 3.08 m) 

is shown in Fig. 16. The maximum efficiency is seen to be about 0.065 at 
t = 0. The efficiency decreases with increasing t and drops to 0.026 at 
t = t* = 1.118, the end of the heat addition period. These values can be 


compared with the equilibrium, constant pressure expansion value of 


111 


0.286, 
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Fig. 16 Efficiency of Piston Driven Accelerator 



Fig. 17 


Rate of Heat Addition: Piston Driven Accelerator 


-62- 


which could be realized if the ramjet were accelerating slowly (this is the 

r-1 


value used in Reference 11). The factor multiplying in Equation (71) 
can thus be considered to be the loss incurred because of the unsteady 
acceleration process occurring in the piston driven device. Physically 
this loss can be traced to the fact that the pressures in the compression 
chamber and heating chamber must be maintained at values increasingly larger 
than P Qo , the pressure at the base of the ramjet, as time increases. 

The effect of L on efficiency (n) can be seen in Eq. (71). Increasing 
values of L decrease efficiency through the second term in brackets. There 
is a limit on how small U can be. If U = If (i.e. no heating chamber volume) 
then n at t = 0 would increase to 0.106 and to 0.045 at "t = 1.118. Of course 
this represents an impossible situation physically, since, as can be seen from 


Eq. (70) that T H -*■ °° 


— — __ 1 U 

if L -*■ 5 . The choice L = 2 e gives j — = 4.13, 

0 0 1 Ho 


still a substantial temperature rise. In general, relatively modest in- 
creases in efficiency can be obtained at the cost of large increases in 


Ho 

The heat addition rate (Q^) as a function of time (t) is shown in 
Fig. 17. From this figure it can be seen that the heat addition rate is 
nonlinear in time and must increase more and more rapidly as the acceleration 
process progresses. The maximum power required for accelerating a 10 kg 
ramjet to 5000 m/sec would be 370 gigawatts (at t = 1.118). 

An overall efficiency for the acceleration process can be defined as 
the ratio of the kinetic energy of the ramjet to the total amount of heat 
added to the heating chamber. Thus 


m r u r /2 

n (y-1 ) 

0 t . 2 


/ Q H dt 


— 2 
t 


(UtltV ' 1 -(!*£) 

c m T rt T j 


(72) 



Equation (72) is valid for any time, t". However it must be remembered 
that for t <_ t* = 1.118, Q^= 0. That is, once the flow chokes at x = 0, the 
power to the heating chamber can be cut off while the ramjet continues at 
constant acceleration until t = 1.708. The period from "t = 1.118 to 
t = 1.708 thus represents a "free" acceleration period during which the over- 
all efficiency continually increases (and n is ro ). At t = t* = 1.118 the 
overall efficiency (n rt ) is 0.0386. At t = t „ this efficiency reaches 
0.0901. This last value represents the overall efficiency for the complete 
acceleration process. During ramjet mode operation, the efficiency of energy 
conversion remains above 40%, even considering shock losses. Consequently, 
it seems that once the ramjet can begin operation it is to be much preferred 
over the piston driven device just described as a mass accelerator. 

Pistonless Single Chamber Accelerator 

In the accelerator design discussed above a separation of the heating 
chamber from the isentropic compression chamber by a piston was assumed. 

This is consistent with the practice common in high velocity light gas gun 
design. The question arises as to what would happen if the piston were 
simply removed and both heating and compression occurred in a single chamber. 
Certainly a measure of design simplification would result. What would be 
the effect on efficiency? 

In attempting to answer this question, the first thing which becomes 
clear is that the analysis of the problem is considerably more difficult 
mathematically. This is because the unsteady gasdynamic flow in the system 
can no longer be modelled as a constant entropy (homentropic) flow. In 
general each fluid particle in both the chamber and the tube will have a 
different entropy value associated with it since it is generated in the 
heating chamber at a different entropy. However, there is no reason to 



expect that the entropy of each particle, once generated, might not remain 
constant. That is, it seems reasonable to assume that = 0, as before, 
for x > 0. The solution for the flow field downstream of x = 0 still 
poses a difficult problem in the method of characteristics, one that cannot 
apparently be solved in closed form, but that is amenable to numerical 
techniques (see Ref. 16). This solution has not yet been attempted. However, 
a reasonable estimate of the performance of the single chamber design, based 
on the uniform acceleration model discussed earlier, can be made. 

As before, it is assumed that a solid body type acceleration of the gas 
and ramjet occurs starting from t = 0 with the ramjet at x = 0. That is, it 
is assumed Eqs. (46) to (54) apply. The basic physical design of the accel- 
erator is also similar. Referring to Fig. 15, removal of the piston leads 
simply and immediately to the design we wish to consider. The parameter L 
now represents the heating chamber length, ahile £ has no significance. The 
ratio of A c to Ay will be taken to be large as in the dual chamber design 
so that a quasi-steady isentropic (not homentropic! ) flow field connects the 
heating chamber with the launch tube. 

Because of the homentropic assumption appropriate in the piston driven 
accelerator, it was possible to replace Eq. (54) by the relationship 
F = (Eq. (56)) and finally develop closed form expressions for P° , p° 
and T°, the pressure, density and temperature at x = 0, as functions of time t. 
in the single chamber design no relationship of this type can be written 
down a priori. In fact, the relationship between pressure and density is 
imp! icit in the complete solution. Though any choice of P(p) will lead under 
the uniform acceleration assumption to the P°, p° and T° required to support 
the solution at x = 0; this flow must in general be generated in the upstream 
region by a combination of piston motion and heat addition in the chamber 



-65- 


immediately before x = 0. For the choice P* = p^, only piston motion was 
required. For the present case where only heat addition is allowed the 
relationship between P - and p at x = 0 is unknown. If TP = "p (isothermal) 
then both piston motion and heat transfer out of the chamber would be 
requi red. 

The solution is pursued by writing down the necessary conservation 
and state equations at 3T = 0, in the heating chamber, and for the region 
connecting them. 

At 7 = 0 ; 

= yp° t (uniform acceleration condition) (73a) 

dt 

P° = p° T° (ideal gas equation of state) (73b) 

In the heating chamber : 

= p" H T h (ideal gas equation of state) (74a) 


dT. . A t -o 


Tu dP. 


1 r , H m, H 

— A — 1 H — 1 + — _ 

dt P H L P H dt 


(conservation of mass) 


• ^ j d^ ^ dPn 

Qh = p « l v-i 


(first law of thermodynamics) 


Across the connecting region ; 


= T° + t (first law of thermodynamics) 


dP, 

dt 


H 


1. lH 

2 T° 


r Li 


r 

“fti 

_ 


A- “0 . ■ 
2 Y t_X^_L£_l 

A Cp H L 


(74b) 


(74c) 


(75a) 

(75b) 


(isentropic condition) 
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The derivation of this system of equations requires extensive algebraic 
manipulations, which have not been reproduced here. The seven equations 
in the seven variables P°, p°, T°, P^, p^» 7^, ^ can be solved numerically 
in a straightforward way by the method of finite differences using a forward 
marching technique in t, recognizing that at t = 0 all of the variables are 
unity except Q^, which is zero. 

The next question which arises concerns how far the integration in t 

should proceed. Certainly if the desired final speed for the ramjet is 

— — u max 

u„,„„, then one answer is t , = u„,„ = — — . For this solution it is 
max max max a„^ 

oo 

only necessary to specify a QQ (or T Q0 ) in order to find the maximum time the 
ramjet can experience uniform acceleration tL, w . Unfortunately, two aspects 
of the problem would be left out with this conclusion. First, the question 
of choking at x = 0, and second the "free acceleration" period which must 
occur after heat addition is terminated. The question of choking does not 
turn out to be of major importance in this case because of the high tempera- 
tures which are present at x = 0. The Mach numer at x = 0 is 



Numerical integration of the governing system of equations (Eqs. (73) through 
(75)) indicates that T° increases rapidly enough so that M° = 1 is not 
approached in reasonable values of t and is never approached for some parameter 
selections (especially small U values). Thus, in this case there does not 
seem to be a "natural" termination of the uniform acceleration solution, as 
there was in the homentropic accelerator. For a given T Q0 , then, higher 
final velocities can probably be achieved in the pistonless design. 

The problem of the "free acceleration" period is more difficult to 
evaluate. Certainly it must take some finite time for information to travel 
along the right running characteristic which intersects the ramjet path in 
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x, t space. In order to find the exact time involved it is necessary to 
apply the method of characteristics to the non-homentropic flow field down- 
stream of x = 0. As was mentioned earlier this has not been done yet. 
Consequently, we are left with making an estimate of what this period might 

be. It is reasonable to expect that it might be of the same order as that 
for the homentropic flow field. 

For the sake of preliminary comparisons between the piston driven and 
pistonless accelerator it will be assumed that heat driven acceleration 
occurs until t = 1.2 and that from t = 1.2 to 1.708 a "free acceleration" 
period occurs. These numbers are arbitrary and will be corrected when the 
characteristic solution is completed. The parameter L", the non-dimensional 
length of the heat addition region, also appears in the equations. As in the 
piston driven case, small values of L give better performance and higher 
maximum temperatures, large values give reduced performance and maximum 
temperatures. The exact dependencies are more obscure in the pistonless 
case. For the preliminary calculations that were performed for the sake 
of comparison with the piston driven accelerator, an L value of 0.05 was 
chosen. This is close to the dimensionless length of the compression 
chamber of the dual chambered accelerator. 

Results of calculations are shown in Figs. 18, 19 and 20. The variations 

of P H , T h and P° with t are shown in Fig. 18. At t = 1.2, the end of the 

heating period, = 2.32. = 3.80, P° = 1.78. Dimensionally these 

correspond to P H = 344 MPa (49,900 psia), T H = 5597°K, P° = 264 MPa (38,300 psia) 

for the exampl e where a 16 cm dia, 10 kg ramjet is accel erated to 5000 m/sec. 

The maximum pressures attained are about half of those attained in the piston 
driven case. Acomparison between maximum system temperatures for the two dif- 
ferent designs is not really possible, because of the arbitrary nature of T^ q 
in the piston driven case. If T Hq were taken to be equal to T^ , then the maximum 
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Fig. 18 Pressure and Temperature: Pistonless Accelerator 


0 . 1 5 1 


0. 10 


u- 0.05 

LU 


7 = 1.4 
D c /D t =4 
L = 0.05 

umax = 5000 m/sec 


0.25 0.50 0.75 I.C 

NON-DIMENSIONAL TIME (T) 


Fig. 19 Efficiency of Pistonless Accelerator 
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Fig. 20 Rate of Heat Addition: Pistonless Accelerator 




system temperature in the pistonless accelerator would be about 40% of that 

of the piston driven accelerator, for the example chosen here. 

In Fig. 19 o the instantaneous efficiency defined earlier is shown 

as a function of time (t). In contrast with the piston driven case (Fig. 16), 

the efficiency increases with time, reaching a maximum value of about 11.9% 

at t = 1.2. The total change in efficiency during the acceleration period 

is 30%. For the piston driven device it decreases by a factor of 2.5 during 

acceleration, starting from a maximum of 6.5%. 

As shown in Fig. 20 the heat addition rate Q u increases nearly linearly 

h 

with time (i.e. power input would have to increase linearly with time to 

provide the uniform acceleration desired). The maximum of Q occurs at t = 1„2 

H 

and corresponds to a power input of 86.8 gigawatts. This compares favorably 
with the maximum power requirement of 370 gigawatts for the piston driven 
accelerator. The overall efficiency of the acceleration, n Q , including the 
"free acceleration" period from t = 1.2 to T = 1.708 is 23.4%, more than 
double that of the piston driven device. 

Conclusions 

Calculations indicate that gun-type, constant base pressure acceleration 
of an electrothermal ramjet vehicle to the speed necessary for effective 
ramjet mode operation is feasible. The efficiency of the acceleration 
process is not high compared with the ramjet itself. Overall efficiences 
(conversion of thermal to kinetic energy) of less than 25% are predicted 
for idealized accelerator operation. These compare unfavorably with effi- 
ciencies greater than 40% for ramjet mode operation, even when shock losses 
are considered. Two types of accelerator designs were evaluated: a dual 

chambered piston driven accelerator, and a single chambered accelerator with 
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no piston. Of the two, preliminary comparisons strongly favor the pistonless 
accelerator. It is mechanically simpler, provides higher efficiencies, has 
the potential for reaching higher maximum velocities, requires a much 
simpler power addition profile (linear vs nonlinear), and reaches lower 
maximum pressures during the acceleration process. 
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A MODEL FOR MERCURY ORIFICED HOLLOW CATHODES. 
THEORY AND EXPERIMENT* 

Daniel E. Siegfried** and Paul J. Wilburf 
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Fort Collins, Colorado 


Abs tr act 

A model is presented which provides a nsefnl 
qualitative description of the basic physical 
processes taking place within a mercury orificed 
hollow cathode and can predict, to first order, 
important cathode operating parameters such as 
emission length and insert temperature. The 
analytical formulation of the model is based on the 
concept of an idealized "ion production region" 
which is defined as the volume circumscribed by the 
emitting portionof the insert. The energy exchange 
mean free path for primary electrons is used as a 
criterion for determining the length (L )of this 
region. An ion production region aspect ratio 
(D/L ) of two is suggested as a design criterion 
for minimizing keeper voltage. The model accounts 
for electrons produced in the ion production region 
both by surface emission and by volume ionization. 
Surface and volume energy balances are used to 
predict plasma density and plasma potential in this 
region. An empirical relation is presented which 
can be used to estimate cathode internal pressure 
(a necessary input to the model) from the discharge 
current and cathode orifice diameter. Calculations 
based on the model are compared with experimental 
results. 

Nomenclature 

a Q - theoretical constant (1.2xl0 6 A/m 2 °K 2 ) 

A c - area of end boundary of ion production 
region (m 2 ) 

A e - insert emission area (m 2 ) 

A s - total surface area of ion production 
region (m 2 ) 

D - insert inner diameter (m) 
d Q - orifice diameter (mm) 
e - electronic charge (Coulombs) 

E — electric field at insert surface (V/m) 

Ip - total discharge current (A) 

I e - insert electron emission current (A) 

- total ion current to cathode surfaces (A) 

- Bohm current density (A/m ) 


* Work performed under NASA Grant NGR-06-002-112 

** Research Assistant, Department of Mechanical 
Engineering, Associate Member AIAA 

t Professor, Department of Mechanical 
Engineering, Member AIAA 


j - field-enhanced, thermionic emission 
current density (A/m 2 ) 

I - normalized emission current ( Eq.l ) 

k - Boltzman’s constant (1.38x10 23 J/°K) 

L - insert emission length (m) 
e 

m - propellant mass flow rate (mA equivalent) 

mj - ionic mass (kg/ion) 

n - electron density in ion production 
6 region (m 3 ) 

—a 

n^ - ion density in ion production region (ra ) 

n Q - total neutjal atom density in ion production 
region (m ) 

P - internal cathode pressure 

q^ - energy flux associated with de-excitation 

of excited states on the emission surface (W/m 2 ) 

q. — energy flux due to plasma radiation (W/m 2 ) 
ph 

6.. - insert thermal power loss (W) 
th 

T - electron temperature in ion production 
6 region (°K) 

T g - insert emission temperature (°K) 

V - plasma potential in ion production 
p region (V) 

Sj - ionization potential (10.4V for Hg) 

e - permittivitity of free space 
0 (8.85x10 F/m) 

c - primary electron energy, equivalent to 
pr plasma potential V p (eV) 

a ^ ^ - electron-electron elastic mean free path (m) 

Agj - effective, energy exchange, elastic 

mean free path for primary electrons (m) 

A In - effective inelastic mean free path 
for primary electrons (m) 

X - primary electron, energy exchange 
mean free path (m) 

- average effective work function (V) 
j» s - average surface work function (V) 

Introduction 
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The results oi previous experimental investi- 
gations of mercury orificed hollow cathodes [1,2] 



showed that —70% of the total cathode current is 
due to surface electron emission from a narrow band 
(~2mm long) on the downstream end of the low work 
function insert and that the dominant surface emis- 
sion process is probably that of field-enhanced, 
thermionic emission. The results also indicated 
that the remaining 30% of the current is due to 
volume ionization which takes place in the very 
intense plasma discharge which is established adja- 
cent to the emitting portion of the insert. A 
phenomenological model describing the emission and 
plasma production processes within the cathode was 
developed based on the results of those experiments 
[2]. The important cathode processes are shown 
schematically m Fig. 1. Surface electron produc- 
tion is from a narrow band on the downstream end of 
the insert, primarily from field-enhanced ther- 
mionic emission. These electrons are accelerated 
across the plasma sheath and, thereby, pick up suf- 
ficient energy to ionize mercury atoms (via mul- 
tistep excitations) in the region adjacent to the 
insert. The total discharge current is made up of 
the sum of the currents due to both surface and 
volume electron production. The ions produced m 
the process leave the production volume at the Bohm 
velocity and are neutralized at internal cathode 
potential surfaces. Ions striking the insert sur- 
face provide the energy input which is required to 
maintain the insert at the emission temperature. 

By defining an idealized ion production volume 
as the region adjacent to the emitting portion of 
the insert (indicated by dashed lines in Fig. 1) 
and assuming that the plasma properties are uniform 
in that region, the physical description above can 
be represented analytically in a rather simple 
form. The analytical model, as initially proposed, 
was presented in Ref. [2] and further developed in 
Ref. [3]. The development of the model is dis- 
cussed in detail in those two references, so only 
the important assumptions and a summary of the 
final results will be presented here. 

Summary of Mod e 1 

~ I 

Table I gives a summary of the equations 
developed in the model along with the physical 
basis for each equation and a brief comment indi- 
cating the important assumptions on which it is 
based. It should be noted that all of the equa- 
tions in Table I are in MRS units EXCEPT for the 
empirical expression for the pressure (Eq. I. 11) 
which has the units indicated in parentheses. Each 
of the equations in Table I will be discussed 
briefly below. 

The total cathode discharge current 1^ is the 
sum of the surface electron emission current I 
plus the current of volume produced ions returning 
to the surface Ij. This is represented by Eq. 1.1 

of the table, where j.. is the field-enhanced 
til 

thermionic current density coming from emitting 
area of the insert (emission length L^) and 
is the Bohm current density of ions crossing the 
boundary of the ion production region which has a 
surface area A . This equation assumes that the 
current due to surface emission processes such as 
cold-field emission, photoemission, and secondary 
emission due to ions and excited state atoms is 
negligible. These emission processes are discussed 
m Ref. [3] where it is shown that the electric 
fields at the insert surface are too low to cause 
appreciable field emission and that photoemission 
is negligible because photons produced within the 


^INSERT 


Fig. 1 Schematic of Ion Production Region 



plasma by resonance de-excitations are effectively 
trapped due to the high neutral atom density, so 
cannot reach the surface m appreciable numbers. 
Secondary emission currents due to ions and excited 
state atoms are estimated to be less than 5% of the 
total current, so they are also neglected. Equation 
1.1 also neglects electron production within the 
orifice region and downstream of the orifice. 
Experimental results indicate that this is a rea- 
sonable assumption since these sources contribute 
only -7% of the total discharge current. 


Equation 1.2 of Table I gives the current den- 
sity J t jj for field-enhanced thermionic emission 
from a surface at temperature T with an effective 
work function of The effictive work function 

j» e is the surface work function f (a material pro- 
perty) minus the reduction in the work function due 
to the electric field E at the emission surface and 
can be calculated using Eq. 1.3. The electric 
field in Eq. 1.3 is the field at the emission sur- 
face associated with the potential drop across the 
plasma sheath. It is estimated based on the 
theoretical analysis presented by Prewett and Allen 
[4] for the double sheath which forms at a hot 
cathode surface emitting electrons into a plasma. 
For the range of operating conditions typical of 
mercnry hollow cathodes, Prewett and Allen's gen- 
eral analysis yields the approximation for the 
electric field given in Table I by Eq. 1.4. This 
approximation holds for plasmas where the ratio of 
plasma potential (V ) to electron temperature (T ) 
is -10 and the normalized emission current 



kT 

<--> 


1/2 


3/2 


-■-1 


( 2 ^) 

e 


( 1 ) 


—3 

is less than ~10 . These criteria are satisfied 

for typical cathode conditions. 

The Bohm current density is estimated from 
an energy balance on the emitting surface. In such 
a balance, the power due to heating from ion neu- 
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Table I 


Summary of Equations Used in Model 


Physical Basis 

Comments 

Eqnation 

No. 

Current balance 

Thermionic electron emission 
from insert only; ion flux 
based on Bohm criterion. 

b ml .*h m J th A e + Ji A s 

1.1 

Field-enhanced 
thermionic emission 

Model neglects all surface 
emission mechanisms, except 
this one. 

2 e ?e 

Jth = ‘o 1 ; kF> 

s 

1.2 

Effective work 
function 

Based on electric field 
at emission surface 

1/2 

r e r s 4ne^ 

0 

1.3 

Double sheath analysis 

Approx, based on theoretical 
analysis in Ref. [41. Bolds 
for J < 10 3 (J from Eq.l) 

rn kT 1 1/2 r eV 1/2 

b “K, ] r 1 * 2 ^’ 

r 

1.4 

Insert energy balance 

Neglects energy input due 
to excited states and plasma 
radiation. 

J * ‘ * I »][ 2 vV't>] 

where a = (V +e .-0 ) 

P i r & 

1.5 

Energy balance on 
ion production region 

Neglects energy loss due 
to excited states and 
plasma radiation. 

t r 5kT I 1 

T P- w. n JlV *• \ 

1.6 

Bohm criterion 

Assumes uniform or 
average plasma properties. 

n ® elk^/mj ] 172 

1.7 

Existence length for 
primary electron 

Criterion for emission/ion 
production region length. 

L = 2k 
e pr 

1.8 

Semi-empirical 

Energy exchange mean free 
path based on results 
of computer model. 

X = 
pr 

6.5x10 17 n 10 3 n V 

« + e_p 

v 2 2.83xl0 23 -1.5n 

P 0 

-i 

1 

Ideal gas lav 

Assumes heavy particle 
temperature equal to 
insert temperature. 

P - n k(T + T ) 

n = s 5- 

O k T 

s 

1. 10 

Empirical 

Neglects orifice plate 
thickness effect. 

V = ■— (13.7+7.82 I )xl0 ~ 3 (Torr) 

d 2 

0 

where ra(mA) and d^tmm) 

1. 11 
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tralization and de-excitation of excited states is 
equated to the power conducted and radiated from 
the surface plus the power removed by emitted elec- 
trons. The equation describing this is 

j A (V + e -<# )+q. A +q . A = 6..+I f* (2) 
l e p l H dx e M ph e th e r e ' ' 

where e is the effective work function of the sur- 
r e 

face, j» is the surface work function (a material 

property), is the thermal power transferred 

away from the surface, is the potential drop 

across the plasma sheath, and e, is the ionization 

potential. The terras q . A and q, A are the 
pn c qz c 

energy input due, respectively, to photons leaving 
the plasma and to surface de-excitation of excited 
atomic states. The photon flux is negligible 
because the resonance radiation is effectively 
trapped within the plasma. The term 4 ( j I A e due to 
excited states is believed to be relatively 
small though not necessarily negligible. Calcula- 
tion of this term and various factors affecting it 
will be discussed in more detail m the last sec- 
tion of this paper. The term <5^ in Eq. 2 is the 
thermal power lost from the insert. In general, it 
is a function of the insert temperature T and must 
be estimated for the specific cathode thermal con- 
figuration on the basis of conduction and radiation 
from the insert. Neglecting the terms 

q^A^ and q^^ in Eq. 2 and solving for the ion 
current density j ^ gives Eq. 1.5 of Table I. In 
expressing Eq. 1.5 in terms of the total discharge 
current Ip, Eq. 1.1 was used along with the fact 
that the total surface area of the ion production 

region is A -A +2A where A is the area of the end 
sec c 

boundary of the region. 

The plasma potential is estimated by a similar 
energy balance on the ion production region. 

Energy is convected into and out of this volume by 
the motion of the various particle species crossing 
the boundaries of the ion production region. This 
particle motion transports not only the random and 
directed kinetic energy of the particles but also 
their excitation potential energy. Again neglect- 
ing the energy flux due to excited state atoms q, 

• dx 

and due to radiation q « leaving the volume, the 
ph 

energy balance can be written simply as 


Ve = e i J x + 2 h 


(3) 


Here V^I^ is the rate at which energy is brought 
into the volume by the primary (surface emitted) 
electrons, e^I is the rate at which energy leaves 
the volume in the form of ionized atoms, and 
(5/2 kT e Ip) is the rate at which energy leaves the 
volume due to the convection of Maxwellian elec- 
trons through the orifice ( 5/2 kT is the enthalpy 
of the Maxwellian electron gas). Other energy 
terms such as the directed kinetic energy and 
enthalpy of the ions and neutrals can be shown to 
be small [3] and are neglected. Using the current 
balance of Eq. 1.1 and solving Eq. 3 for V results 
in Eq. 1.6 of Table I. p 

Calculation of the electric field E using 
Equation 1.4 requires values for the plasma density 
n and the electron temperature T . Assuming that 
tie plasma density (n^ ~ n^) is uniform throughout 


the ion production region and that the ions leave 
the region at the Bohm velocity (kT^/m^) , an 
equation for the ion current density from the pro- 
duction region can be written. Rearrangement of 
this equation gives Eq. 1.7 where j i is the ion 
current density from Eq. 1.5. The electron tem- 
perature is also needed here. However, neither Eq. 
1.4 nor 1.7 is very sensitive to the electron tem- 
perature because it appears as a square root. 
Probe measurements in the ion production region 
show that over the normal range of mercury hollow 
cathode operation the electron temperature is rea- 
sonably constant at 0.71 + 0.1 eV. Therefore, it 
is suggested that this measured value be used m 
the equations. 

To complete the model one requires a means of 
estimating the length of the ion production region. 
This length is determined by the collisions! 
processes taking place within the ion production 
region. Energy input to the ion production region 
occurs mainly through the current of primary elec- 
trons emitted from the insert. The dominant energy 
exchange paths for a primary electron entering this 
region are inelastic collisions with a mercury atom 
(typically 10® collisions per second) and elastic 
collisions with Maxwellian electrons (typically 10° 
collisions per second). The elastic collision 
cross-section for an electron-electron collision is 
proportional to the inverse square of the projec- 
tile electron's energy. Therefore, once a primary 
electron gives up energy, either by small angle 
elastic or by inelastic collisions, it’s collision 
rate increases and it is rapidly thermalized, los- 
ing it's identity as a primary. The fact that the 
primary electrons are required for the excitation 
reactions which sustain the plasma discharge sug- 
gests the length of the ion production region is 
probably on the order of the primary electron 
energy exchange mean free path. Primaries can lose 
energy both by elastic collisions over a mean free 
path of X and by inelastic collisions over a mean 
free path of Therefore, the effective mean 
free path for the loss of a primary electron would 
be 


X pr " (1/X ln + 1/X E1> 


-1 


(4) 


From where they are created, primary electrons will 
be scattered upstream to some extent by elastic 
collisions. Therefore, a reasonable criterion for 
the length of the ion production region is probably 
one to two primary electron energy exchange mean 
free paths. It will be seen in the next section 
that a length of two mean free paths gives the best 
agreement with experimental results. This cri- 
terion is indicated in Table I by Eq.1.8. 

The primary electron mean free path criterion 
as well as the concept of an ion production region 
are clearly idealizations. The extent of the ion 
production region in reality will not be sharply 
defined; properties will not be uniform throughout 
the region and ionization will occur upstream of 
the boundary. In addition, ions produced in the 
region will diffuse upstream, heating the insert in 
the region upstream of the boundary. So far the 
diffusion of ions upstream has not been discussed, 
although to be consistent with the model, the 
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plasma density oust fall off rapidly upstream of 
the boundary. If this were not the case, signifi- 
cant ion heating of the insert and related electron 
emission would be expected to occur upstream of the 
boundary defined by the mean free path criterion. 
Experimental plasma density profiles show that, in 
fact, the axial plasma density drops off exponen- 
tially in the upstream direction. This indicates 
that significant ion heating probably does not 
extend far upstream of the ion production region 
and is consistent with the observed rapid fall off 
in surface temperature upstream of the emission 
region. This observed behavior is also in agree- 
ment with calculations based on ambipolar diffusion 
of the ions and electrons. Such calculations show 
that, by itself, volume recombination in the plasma 
can result in a density decrease upstream of the 
production region more rapid than that observed 
experimentally. In the real situation, of course, 
both volume ionization and wall recombination are 
also taking place in conjunction with volume recom- 
bination to determine the actual plasma density 
profile in this upstream region. 



Fig. 2 Effect of Total Neutral Atom Density 
on Primary Electron Inelastic Mean 
Free Path 


In order to use Eq. 4 to estimate the emis- 
sion length it is necessary to be able to calculate 
the elastic (Ajjj) and inelastic (Aj^) energy 
exchange mean free paths. Elastic electron-ion and 
electron-atom collision frequencies are relatively 
low. Therefore, the effective elastic mean free 
path for energy exchange is approximately equal to 
the electron-electron mean free path, or 


A =A = 
El ee 


2 

e 

_E£_ 


6.5x10 


.-17 


(5) 


siderable range, the mean free path for inelastic 
collisions was found to be sensitive only to total 
neutral density and to primary electron energy. 
The Maxwellian electron density n was also found 
to have a slight effect on ^j n * However, this 
effect is rather small amounting to an increase m 
Xj n of less than 20% when electron density is 
increased by a factor of five. For our purposes 
here, this effect can be neglected. 

Recognizing that the mean free path varies as 
the inverse ot the neutral density n , the results 
of Fig. 2 are fitted with the following expression 


Here the expression for X is based on a Coulomb 
ee 

collision between a primary electron with energy 

e and a low energy Maxwellian electron [5], and 
p r 

n is the Maxwellian electron density. The inelas- 

tic 

mean free path Is given by the following expression 


Ay = 

In 


(6) 


Z n 


a “ip 


where is the collision cross-section for pro- 
duction p of excited state p from a target particle 
of type a having a density n . The summation in 
Eq. 6 is made over all of tiie important excitation 
reactions in the ion production region. Evaluation 
of using Eq. 6 requires the densities of the 
various excited states of mercury in the ion pro- 
duction region as well as the cross-section for all 
of the important reactions. Peters [6] has com- 
piled the necessary collision cross-section data 
and developed a computer model which calculates 
excited state densities for a mercury discharge. 
The application of Peters' computer model for this 
purpose is described in detail in Ref. [3], so only 
the important results will be presented here. 


The computer model was used to calculate 
excited state densities and inelastic mean free 
paths over a wide range of input parameters 
typical of hollow cathode conditions. The results 
of the computations are shown in Fig. 2 where Aj. 
is plotted as a function of total neutral density 
n using primary energy as a parameter. Although 
a?l of the input parameters were varied over a con- 


In 1 n 

o 


■1.5) 


10 3 e 


pr 


(7) 


This expression fits the results shown in Fig. 2 
within + 5% over the full range of the parameters 

except for values ot equal to 6.0V where A^^ is 

overestimated by ~20% in the intermediate density 
range. Using Eq.5 for the elastic mean free path 
and Eq. 7 for the inelastic mean free path, the 
effective primary electron mean free path for 
energy loss due to both elastic and inelastic col- 
lisions A is given by Eq.1.9 of Table I. It has 
also been P assumed in this equation that the primary 
energy e pr is equal to the plasma potential V p . 
Equation 1.9 provides an easy means of estimating 
the insert emission length. Of particular impor- 
tance is the fact that over the normal range of 
cathode conditions, the results of Eq. 1.9 are 
dependent mainly on tho neutral density and are not 
very sensitive to the plasma potential (V is nor- 
mally 8— 12V ) or electron density. This P enables 

one to make a reasonable estimate of A based only 

pr 

on typical plasma conditions and the neutral den- 
sity, It is worth pointing out that the curves in 
Fig. 2 are relatively flat at densities 
greater than ~2 x 10 16 cm 3 resulting in typical 

values ot A near one millimeter. This is in good 
pr 

agreement with the experiments, which show that for 

internal pressures over a few Torr (n - 2xl0 3 ** 
3 o 

cm ), the emission is confined to a region ~2 mm 

long or less on the downstream end of the insert 
but that the region extends significantly upstream 


80 



MICA SPACERS 


» 


at pressures below a few Torr. Comparison with 
additional experimental results will be presented 
in the next section. 

Finally, it is necessary to estimate the total 
neutral density n Q in the region upstream of the 
orifice. The presence of the cathode orifice plate 
simplifies this problem somewhat because most of 
the pressure drop is across the orifice and the 
pressure within the cathode cavity upstream of the 
orifice is essentially constant. The total pres- 
sure at any point m the cathode is the sum of the 
partial pressures of each species. Using the ideal 
gas law this can be expressed as 


P “ n kT +n. kT, +n„ kT 
e e i i o o 


( 8 ) 



where k is Boltzmann's constant, n is the density, 
T is the temperature, and the subscripts e, i, o 
refer to electrons, ions, and neutral atoms respec- 
tively. Given the local pressure, each of the tem- 
peratures, and the plasma density, Eq. 8 can be 
solved for the neutral gas density. Invoking 
quasi-neutrality of the plasma (i.e., n^ ~ n.) and 
equilibrium of the heavy particles (ions and neu- 
trals) with each other at the insert temperature, 
Eq. 8 takes the form given in Eq. I. 10. The local 
pressure P used in Eq, I. 10 is determined from a 
correlation of experimental results based on actual 
pressure measurements made over a wide range of 
mass flow rates, discharge currents, and orifice 
diameters [3], This correlation is shown in Table 
I as Eq. I. 11 where the units are P (Torr), m (mA 
equivalent), d (mm) and I (A) . It is believed to 
be accurate to within +30%. 

Given the mass flow rate m, the orifice diame- 
ter d Q , the insert diameter D, the total discharge 
current L, the surface work function , the 
insert thermal power loss 6 and the physical 

constants (e, k, m,, e , e,, a ), the equations in 
X O 1 o 

Table I can be solved for all of the other parame- 
ters except the electron temperature. As discussed 
above, a value of 0.71 eV is a reasonable assump- 
tion for the electron temperature. Because Eq. 1.2 
can not be solved explicitly for the surface tem- 
perature T , the solution of these equations is 
necessarily iterative. However, the equations con- 
verge to a solution very rapidly (five place accu- 
racy in six or seven iterations). The results of 
calculations based on this set of equations will be 
compared with experimental results in the next sec- 
tion. 

Com p arison with Experim ent 

An earlier paper describes an experiment in 
which a special quartz tube test cathode was used 
to determine the magnitude and origin of the vari- 
ous hollow cathode currents [21. The predictions of 
the model will be compared here with the results of 
that experiment. 

A schematic of the test cathode used in the 
experiment is shown in Fig. 3. The cathode was 
fabricated from a quartz tube such that all of the 
internal surfaces were covered with tantalum foil 
and each surface of interest was isolated as a 
separate electrode. The insert was 2.2 mm long and 
had a diameter of 3.9 mm. The cathode had an ori- 
fice diameter of 0.96 mm. Both the insert and the 
cathode tube were coated with the low work function 


chemical R-500*. The cathode discharge was coupled 
to a cylindrical anode which was completely 
enclosed within a stainless steel ground tube 
covered on its downstream end by a fine mesh 
screen. The current from each of the surfaces 
shown in Fig. 3 was measured separately for total 
discharge currents I Q of 1.3, 2.3, 3.3 and 4.3A at 
a mercury flow rate m of -100mA. Data were also 
collected for an emission current of 3.3 A at 
internal cathode pressures P ranging from 1.3 to 
3.5 Torr. The internal pressure was measured using 
a U-tube manometer and the insert temperature was 
measured using a micro-optical pyrometer. 

Emission Length 

A critical test of the model is its ability to 
predict the length of the insert emission region 
since this parameter strongly affects the predicted 
insert surface temperature for a given discharge 
current and surface work function. The model 
assumes that the insert emission length is between 
one and two primary electron mean free path 
lengths. The primary electron mean free path 1 
was calculated with Eq. 1.9 for the experimental 
conditions. In making these calculations, the neu- 
tral density was estimated from the measured 
cathode pressure using Eq. I. 10 and the primary 
electron energy was taken to be 8.7 eV. The pri- 
mary electron energy could have been estimated 
using Eq. 1.6, but it will be seen shortly that 
this underestimates the plasma potential by -2V. 
The value of 8.7 eV was chosen on the basis of 
plasma property measurements [3] which showed that 
over the range of normal operating conditions for 
this size cathode, the plasma potential in the ion 
production region is nearly constant at 8.7 +0.5 V. 

The results of the mean free path calculation 

discussed above are plotted as the solid curve 

(2k ) in Fig. 4. Insert emission lengths deter- 

pr 

mined from the experimental results are indicated 
in the figure as circles. The agreement between 
the experimental data points and the calculated 
curve is quite good and supports the assumption of 
the model that is approximately two. The 

data points above four Torr would have been in even 
closer agreement with the curve except that there 
was no way in this experiment to discriminate emis- 
sion lengths shorter than the 2.2 mm long insert 
segment. Insert emission lengths, for conditions 

* A double carbonate mixture - (Ba/Sr) COj - 
manufactured by the J.R. Baker Chemical Co., 
Phillipsburg, New Jersey. 
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INTERNAL PRESSURE (Torr) 

Fig. 4 Comparison of Emission Length with 
Primary Electron Energy Exchange 
Mean Free Path 

where emission was taking place not only from the 
insert (surface 2) but also from the cathode tnbe 
(surface 1), were estimated by assuming that the 
current density for emission from the tnbe was the 
same as the measured average current density of the 
2.2 mm long insert at that same condition. 

It is worth noting hero that, for operation at 
a typical internal cathode pressure of ~3 Torr, the 
2mm emission length is the same as the cathode 
insert radius. These conditions correspond to a 
pressure-diameter product of ~1 Torr-cm. It is 
interesting to compare this with experimental 
results reported by Lidsky et al. [7] for an open 
tube cathode. Working with a number of different 
gases (Hj, H , A, N^,), they determined experimen- 
tally that the emission region locates itself where 
the local pressure-diameter product is -1 Torr-cm. 
A pressure-diameter product of this magnitude 
corresponds to a primary electron mean free path 
which is on the order of the cathode radius. In 
the open tube cathode, the local pressure varies 
significantly along the tube. It is reasonable in 
this situation that the ion production region 
locates itself where conditions are such that pri- 
mary electrons can reach the centerline of the 
cathode. In the orificed hollow cathode, on the 
other hand, most of the pressure drop occurs across 
the orifice plate so that the pressure within the 
tube is rather constant along the axis of the 
cathode and can be adjusted independently without 
significantly affecting the emission location. 
However, a primary electron mean free path on the 
order of the insert radius is probably a desireable 
operating condition here as well. This is so even 
though the orificed cathode can be forced to 
operate under conditions both with shorter and 
longer primary electron mean free paths. That a 
primary electron mean free path on the order ot an 
insert radius is a desireable operating condition 
is also suggested by experimental results [1,8] 
which show that a keeper/discharge voltage minimum 
will occur at a pressure-diameter product of a few 
Torr-cm. 

The above discussion suggests that, for design 
purposes, the insert should be chosen to have a 
radius on the order of a few primary electron mean 
free paths. This would correspond to an insert 
radius which is about the same as the emission 
length and an ion production region which has an 


aspect ratio (D/L ) of approximately two. Such an 
aspect ratio is e also convenient with regards to 
application of the model since it should result m 
an ion production region with fairly uniform plasma 
properties (a basic assumption of the model). 

Other Predictions 

The results of the experiment can also be com- 
pared with the model's prediction of insert tem- 
perature, plasma density, and plasma potential. In 
order to do this, calculations using the model were 
performed based on the following considerations 

1) Rather than the mass flow rate, the measured 
pressures were used as input to the model. 
Equation I. 11 could have been used to determine 
the pressure from the flow rate and orifice 
diameter, but in this case it underestimates 
the pressure by about 20%. A probable explana- 
tion for this discrepancy is that the effective 
diameter of the tantalum, foil-lined orifice 
was slightly smaller during operation than its 
measured diameter when cold. A small error in 
the orifice diameter results m a large error 
in pressure because it appears in Eq. IV. 1 as 
an inverse square. 

2) The surface work function for the insert was 
not known exactly, although it was believed to 
be between 1.8 and 2.0 eV. The value of 1.94 eV 
used in the calculation was arbitrarily chosen 
from this range to give agreement between the 
measured and the calculated temperature at the 
operating conditions where 1^ was 3.3 A and m 
was -100 mA. The comparison U be tween calculated 
and measured temperatures will, therefore, be a 
relative one valid mainly for checking the 
functional dependence on discharge current and 
pressure predicted by the model. 

3) The thermal power loss used for the calcul- 
ations was taken from the dashed curve of Fig. 
A.l in Appendix A. 

4) The electron temperature was assumed to be the 
average measured value of 0.71 eV. 

5) The ratio of emission length to mean free path 

L /X was assumed to be two. 

e pr 

The results of calculations using the model 
based on the assumptions discussed above are shown 
in Figs. 5, 6, and 7. Figures 5a and 5b show the 
effect of discharge current and pressure on the 
emission temperatures predicted by the model for 
the assumed average surface work function of 1.94 
eV. It should be recalled here that agreement with 
the experimental temperature at a current of 3.3A 
was assured by selecting a value of 1.94 for p . 

The significant feature of Fig. 5a, therefore? is 
that the model accurately predicts the effect of 
discharge current on the emission temperature. 

Such good agreement is not obtained for the 
effect of pressure on the emission temperature. 
Fig. 5b shows emission temperatures plotted as a 
function of internal pressure for a discharge 
current of 3.3 A. Unfortunately, at the pressure 
conditions indicated by the solid circles, insert 
temperatures were not measured m this experiment. 
The data indicated by the solid symbols were, 
therefore, estimated using the results of another 
experiment [1] which showed that for a cathode of 
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DISCHARGE CURRENT (A) 

Fig. 5a Effect of Discharge Current on 

Calculated and Measured Emission 
Surface Temperatures 



0 12 3 4 5 6 


INTERNAL PRESSURE (Torr) 

Fig. 5b Effect of Internal Pressure on 

Calculated and Measured Emission 
Surface Temperatures 

* similar construction operating at a discharge 

current of 3.3A the maximum insert temperature 
decreased by ~ 7°C per Torr when internal pressure 
vas increased. However, it is believed that this 
s estimate is reasonably reliable, so it would not 

account for the differences between the two curves 
in Fig. 5b. The main reason for the differences 
between the two curves is that the temperatures 
indicated by the circles are the maximum emission 
temperature while the calculated temperatures (cir- 
cles) are average values based on the assumption of 
uniform emission over the entire emitting length. 
While for short emission lengths, the uniform emis- 
sion assumption is probably valid, the assumption 
apparently breaks down as the length of the emis- 
sion region increases. This suggests that, for low 
pressures, emission temperature gradients become 
important and the simple, lumped parameter model 
does not provide an accurate description of the 
emission temperature. 

The ion production region plasma density 
predicted by the model^j^ plottfd in Fig. 6a as a 
function of discharge durrentV rl, £or comparison, the 
plasma density n^ adjacent to the upstream side of 
' the orifice plate is plotted as the circles. These 

data points (n^) were estimated from the measured 
current 1^ based on the Bohm criterion (Eq. 1.7) and 
represent the average plasma density adjacent to 
the downstream boundary of the ion production 
v region. Plasma densities are shown in Fig. 6b 

plotted as a function of internal pressure for a 
discharge current of 3.3A. In both cases (Figs. 6a 
and 6b), the agreement between experimental values 
and those predicted by the model is reasonably good. 



DISCHARGE CURRENT (A) 


Fig. 6a Effect of Discharge Current on 
Calculated and Measured 
Plasma Densities 



INTERNAL PRESSURE (Torr) 


Fig. 6b Effect of Internal Pressure on 
Calculated and Measured 
Plasma Densities 


although in Fig. 6b the curve shapes for the meas- 
ured and calculated values are different. This 
difference may be attributable, at least partially, 
to the assumption in the model of a uniform ion 
production region. This assumption is expected to 
hold best for operation at conditions where the 
primary electron mean free path is on the order of 
the insert radius. As the internal pressure is 
increased beyond this point (a few Torr), the pri- 
mary electrons have a lower probability of reaching 
the centerline of the cathode, so that ion produc- 
tion is increasingly confined to an annular region 
bounded by the emission surface (as opposed to the 
assumed cylindrical volume). Such a reduction in 
volume could account for the fact that the average 
measured values of the plasma density adjacent to 
the orifice plate are nearly constant for pressures 
above a few Torr (Fig. 6b). 

The fraction ot the total discharge current that 
is due to volume electron production processes 
depends mainly on the plasma density. For the 
plasma densities predicted by the model, this frac- 
tion ranged from 0.30 to 0.34. This is in good 
agreement with results based on measured currents 
which showed current fractions due to volume ioni- 
zation ranging from 0.27 to 0.31. 


83 



Finally, the plasma potential predictions of 
the model are shown in Figs. 7a and 7b where they 
are plotted, respectively, as a function of 
discharge current at a constant mass flow rate 
(~100 mA) , and as a function of internal pressure 
at a constant discharge current (3.3A). Although 
the plasma potential was not measured in this 
experiment, results from Langmuir probe measure- 
ments made during another experiment on a similar 
cathode showed that the plasma potential in the ion 
production region was ~8.7V [3]. The average meas- 
ured value of 8.7V from that experiment is indi- 
cated in the Figs. 7a and 7b as the dashed line. 
The value is shown as a constant (horizontal line) 
because the probe measurements indicated no clear 
correlation of the plasma potential in the ion 
production region with either current or pressure. 
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Fig. 7a Effect of Discharge Current on 
Calculated and Measured 
Plasma Potentials 
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Fig. 7b Effect of Internal Pressure on 
Calculated and Measured 
Plasma Potentials 


It should be recalled that both the plasma 
density and plasma potential predictions are based 
on energy balances for the ion production region 
and neglect the energy term due to the the flux of 
excited atoms. A rough estimate of the energy flux 
due to excited states was made in Ref. [3]. The 
estimate was based on excited state densities cal- 
culated using Peters' computer model and assumed 
that the excited atoms leave the ion production 
region at their thermal velocity carrying with them 
their excitation potential energy. Furthermore, it 


was assumed that all of the excitation energy is 
deposited at surfaces outside of the boundary and 
that the de-excited atoms return with only their 
random thermal energy. The excited state energy 
flux based on those assumptions was found to be 
rather large [3] and gave predictms for both 
plasma density and plasma potential which were not 
realistic. On the other hand, the results in Figs. 
6 and 7 indicate that the model predicts plasma 
densities and plasma potentials of the right magni- 
tude when the excited state power term is 
neglected. (In both cases the model would provide 
better agreement if some fraction of the calculated 
excited state power were included - 0 to 301), The 
over estimation of the excited state energy flux 
may be due to a number of factors. For example, it 
may be that the excited state densities near the 
boundary and, therefore, the particle fluxes are 
actually over estimated. Another possibility is 
that the particle fluxes are estimated correctly 
but that there is some mechanism for returning 
energy to the volume, such as reflection of excited 
particles or of resonance radiation from cathode 
surfaces. However, even in light of these uncer- 
tainties associated with calculation of the excited 
state energy flux, the model provides a useful and 
reasonable qualitative description of the basic 
physical processes taking place within the cathode. 

Based on the above experimental results, it is 
suggested that for design purposes the average 
measured values oi the electron temperature (0.71 
eV) and plasma potential (8.7V) should be used when 
making calculations with the model. The results 
also indicate that the best agreement between model 
and experiment is obtained when the ion production 
region aspect ratio (D/L ) is on the order of two. 
This is also a condition which is expected to pro- 
vide good cathode performance. 


Cone lusions 

A phenomenological model has been presented 
which describes the physical processes underlying 
the operation of mercury orificed hollow cathodes 
of the type used in ion thrusters. This descriptive 
model is in good qualitative agreement with the 
experimental results indicating that our basic 
understanding of the important physical processes 
for these devices is essentially correct. By assum- 
ing an idealized ion production region within which 
most ot the plasma processes are concentrated, this 
phenomenological model has been expressed analyti- 
cal as a simple set of equations which relate 
cathode dimensions and specifiable operating condi- 
tions, such as mass flow rate and discharge 
current, to such important parameters as insert 
temperature and plasma properties. Comparison with 
experimental results shows that, if the excited 
state energy flux is neglected, the model provides 
reasonably accurate predictions of emission length, 
emission surface temperature, plasma density and 
fraction of discharge current due to volume ioniza- 
tion. The plasma potential prediction of the model 
is rather low; and it is suggested that the average 
measured value of 8.7V be used in place of Eq. 1.6 
when making calculations using the model. An ion 
production region aspect ratio (D/L ) of two is 
suggested as a design criterion for minimizing 
keeper voltage. 
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APPENDIX A 

Thermal Power Loss fr om Emissi on Surfa ce 

The rate of heat transfer (Q ) away from the 
emitting portion of the insert due to conduction, 
convection and radiation is required for the energy 
balance calculations discussed above. An accurate 
value for the power loss <5.. for the experimental 

til 

cathode is not easily determined, although bounds 
can be placed on its value and a reasonable esti- 
mate can be made. A maximum value can be deter- 
mined by assuming that end of the quartz tube which 
holds the insert segment is at the same temperature 
as the segment (perfect thermal contact). The 
power loss is then calculated based on conduction 
down the quartz tube and radiation from its outer 
surface to ambient (T & ~ 300°K) . The quartz tube 
would then be analogous to a cylindrical fin one 
end of which is at the insert temperature. The 
results of this analysis are plotted in Fig. A.l as 
a function of insert temperature and are indicated 
by the upper (maximum) curve. The minimum value of 
is calculated by assuming radiation from both 
outer and inner surfaces of the insert to surfaces 
at 700 C and by neglecting all losses due to con- 
duction and convection from the insert. A tempera- 
ture of 700 C was selected because the adjacent 
surfaces - non-emitting portion of insert, start-up 
heater, etc. - were at a temperature less than the 
minimum temperature detectable with the optical 
pyrometer (~700°C) used in making temperature meas- 
urements. The minimum value for 6 . calculated in 

tn 

this way is plotted as the lower curve in Fig. A.l. 



INSERT TEMPERATURE (°C) 


Fig. A.l Thermal Power Loss for 2 mm Long 
Tantalum Foil Insert Segment. 

Finally, a more probable value for 6.. is calcu- 

tn 

lated by assuming direct radiation from the insert 
in addition to some contact between the insert and 
the quartz tube. The heat transfer between the 
insert and tube is estimated to be due to a combi- 
nation of direct contact and transfer via the 
intervening mercury vapor. The heat transferred to 
the tube is then assumed to be conducted down its 
length and radiated from its median diameter to 

surfaces at T ~ 300°K. The results of this 
a 

analysis are plotted in Fig. A.l as the dashed 
(best estimate) curve. 
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Abstract 

An experimental study of ion beamlet vectoring 
Is described. The results can be used to design 
accelerating systems in which the direction and 
focus of emerging beamlets are Important. Deflec- 
tion and divergence angle data were collected for 
two-grid systems as a function of the relative dis- 
placement of these grids. At large displacements, 
accelerator grid impingement becomes a limiting 
factor and this defines the useful range of beamlet 
deflection. Beamlet deflection was found to vary 
linearly with grid offset angle over this range. 
Values of deflection-to-offset angle ratio and use- 
ful range of deflection are given as functions of 
grid-hole geometries, perveance levels, and accel- 
erating voltages. The divergence of the beamlets 
was found to be unaffected by deflection over the 
useful range of beamlet deflection. The grids of a 
typical dished-grid ion thruster are examined to 
determine where over the grid surface the grid off- 
sets exceed the useful range. This defines the 
regions on this surface where high accelerator grid 
impingement Is probably occurring. 


Nomenclature 

d a accelerator grid hole diameter 

d s screen grid hole diameter 

E! electric field, acceleration region 

E 2 electric field, deceleration region 

f thin lens focal length 

h dish depth 

j current density 

J beam current per hole 

£e effective acceleration length 

l acceleration length 

£g grid separation distance 

m-| ion mass 

NP/H normalized perveance per hole 

q charge 

R net-to-total acceleration voltage ratio 

R c radius of curvature 

t a accelerator grid thickness 

t s screen grid thickness 

Vf) net accelerating voltage 

V 7 total accelerating voltage 

a beamlet divergence half-angle 

B beamlet deflection angle 

B u useful beamlet deflection range 

6 grid offset angl e 

6 U useful grid offset angle range 

e relative grid displacement 

e a fractional change in accelerator grid 
hole spacing 

e 0 free space permittivity 
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* Research Assistant, Department of Mechanical 
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Introduction 

The study of grid translation effects on ion 
beam optics has two important applications in the 
design of electrostatic thrusters. First, for a 
dished grid pair, hole misalignments are used to 
redirect the emerging beamlets in order to minimize 
off-axis thrust. A study that would yield correla- 
tions or a model of this phenomenon would facili- 
tate design of these grids. Second, for some 
applications a vectorable beam is desired^ and grid 
translation Is one method of producing a moderate 
thrust vectoring capability. In addition, beam 
focusing may be desirable to increase the current 
density in such ground-based applications as 
sputtering . 2 Again a model or correlations descri- 
bing the effects of hole misalignment on grid 
optics would be useful in designing the hardware 
needed for these tasks. 

Work done by Aston"* has shown the effects of 
such parameters as grid separation, grid hole dia- 
meter and net-to-total accelerating voltage ratio 
on the optical performance of a two-grid system with 
aligned apertures. Optical performance in this 
context was defined as the beamlet divergence vari- 
ation with beam current as well as the maximum 
perveance that could be obtained before accelerator 
grid impingement became excessive. The present 
study of grid translation is Intended to examine 
the effects of grid geometry and operating para- 
meters on grid performance for vectored beams. Of 
particular importance in this study are the degra- 
dation of optical performance as a function of re- 
lative grid displacement, the functional relation- 
ship between beamlet deflection angle and grid 
displacement, and the maximum usable beamlet de- 
fl ection. 

Definitions and Linear Optics Theory 

Figure la shows a two-grid configuration in 
which the accelerator grid hole is displaced a dis- 
tance e relative to the screen grid hole. As a 
consequence the beamlet undergoes a deflection 
through an angle B (degrees) with respect to the 
screen grid centerline, while diverging with a 
half-angle a (degrees). Figure la also illustrates 
pertinent grid dimensions including screen and ac- 
celerator grid aperture diameters (d and d ), grid 
thickness [t and t a ) and the grid separation dis- 
tance [in). S The grid offset angle S also shown in 
Fig. la is defined in terms of the relative grid 
displacement e by the equation: 



Figure lb shows typical potential variations 
with position for the two-grid aperture configura- 
tion of Fig. la. As the figure suggests the ions 
are first accelerated through a total potential, 
V 7 , and then decelerated to a final ion velocity 
corresponding to V^. The ratio of net-to-total 
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where is the maximum current per hole, and d s 
is the screen hole diameter. For the argon propel- 
lant used in this study, the maximum normalized 
perveance per hole is 6.8 x 10' 9 A/V 3 / 2 . It is 
noted that beamlet deflection can only be achieved 
at values of NP/H significantly below this limit. 


Beamlet deflection can be estimated using 
linear optics theory,®*® and modeling the accelera- 
tor aperture of a two-grid system as a thin lens 
with a focal length' f. The focal length is given 
by 


f 



(6) 



Fig. 1. Beamlet Nomenclature and Potential 
Variation. 

voltage is designated as the R-value. 


where Vj is the potential of the accelerator grid 
with respect to the screen grid and Ej, E 2 are the 
potential gradients on the acceleration and decel- 
eration sides of the accelerator grid, respectively. 
The beamlet deflection angle, B is then given by: 


B = j , (7) 


where c is the relative grid displacement. When Ej 
Is assumed to be equal to the total accelerating 
voltage divided by the grid separation distance, 
this becomes 

• (8) 
g 



( 2 ) 


In the case where E 2 «Ej this simplifies to: 


The results of this study will be presented as 
a function of the normalized perveance per hole at 
which the grids are operating. This quanity desig- 
nated herein as "NP/H" is based on the theoretical 
Child's law current density limit for ion extrac- 
tion between parallel plates due to the accumulation 
of space charge in the intervening gap. This law is 
expressed as: 


J 


4 

' 9 e o 



3/2 


(3) 


or 



(9) 


B = - 14.3 (e in degrees). (10) 

g 


The minus sign indicates that beamlet deflection 
occurs in the direction opposite to that of accel- 
erator grid displacement. It should be noted that 
in Eq. 8, E 2 and Ei are of opposite sign so when E 2 
is significant, the ratio (e/e) should be greater 
than the value given by Eq. 10. 


where j Is the maximum current density, e 0 is the 
permittivity of free space, q is the charge of the 
ion, m, is the ion mass, Vj is the total accelerat- 
ing voltage and z is the acceleration length. If 
one considers the effective acceleration length 
appropriate to Fig. la to be the one from the 
plasma sheath to the intersection of the screen 
hole centerline and the upstream surface of the 
accelerator grid then the length in Eq. 3 is given 
by 


*e = + V 2 + d s W 

the theoretical limiting value of normalized per- 
veance per hole is given by 


Apparatus 

The beamlet vectoring study has been conducted 
on an 8 cm dia. divergent field Ion source in which 
the grids have been masked down to a 1 cm 2 active 
area. This masking of the grids results in very 
uniform plasma properties at each aperture and 
hence a uniform beam. The thruster uses a tungsten 
filament for both the main and neutralizer cathodes 
and it operates with argon as the propellant. 

The grids are fabricated from graphite sheets 
with a nineteen hole array arranged in a hexagonal 
pattern on them. A screen grid hole diameter of 
2.06 mm was selected for the study and mica sheets 
were used to maintain the spacing between the grids 
While the accelerator is held fixed, the screen 
grid can be translated along one axis during thrus- 
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ter operation by a mechanical system In the manner 
suggested by Fig. 2. The screen grid was selected 
as the translating grid, because such an arrange- 
ment facilitated grid spacing and accelerator grid 
alterations. 


The beam emerging from the grids is probed by 
a collection of molybdenum, flat plate probes each 
1 cm 2 in area, arranged In the manner suggested in 
Fig. 2. Twenty of these probes are situated along 
a horizontal line 14 cm downstream from the grids. 

In addition, a line of probes extends vertically 
downward from each end of the horizontal probe rake. 
By using these three probe rakes It is possible to 
intercept the full beam profile over the complete 
range of grid translation. A 70% open area screen 
is placed In front of the probes, and is biased 
24 v negative relative to ground in order to re- 
flect neutralizer electrons. The probes themselves 
are biased 18 v positive relative to ground to re- 
flect low energy, charge-exchange ions produced 
downstream of the grids. The rake assembly is de- 
signed so it can be rotated in the manner suggested 
in Fig. 2 to facilitate rake axis alignment with 
the beam axis In the plane perpendicular to the 
direction of grid translation. 


ELECTRON SHIELDING 



RAKE ASSEMBLY ROTATION 
TO FACILITATE ALICNMENT 
WITH BEAM 



■A 


Fig. 2. Deflecting Beamlet Test Apparatus. 


Procedure 


A complete set of beam profiles is collected 
by first achieving stable Ion source operation with 
the screen grid displaced to produce a substantial 
negative ion beam deflection. Once the ion beam 
has been probed, the screen grid is displaced in 
0.05 mm increments and the beam is probed after each 
displacement. This process continues until the beam 
Is displaced sufficiently positive so the accelera- 
tor grid impingement current reaches n. 25% of the 
beam current. Since a grid displacement of 0.05 mm 
has been shown to produce measurable beam deflection, 
careful alignment of the grids In the direction 
perpendicular to the direction of grid translation 
and reduction of play in the movable grid are cru- 
cial to the generation of accurate data. 

Beam current density profiles are analyzed 
using a computer routine. This routine first uses 
the current density data obtained from the three 
segment probe rake to compute the current density 
on a circular arc located 17 cm from the grid 
center. From this profile the point of maximum 
current density is identified. An initial estimate 
of the deflection angle is then computed as the 


angle between the thruster axis and the line from 
the screen grid center through the point of maximum 
current density. Next the two segments of the cur- 
rent density profile on either side of the maximum 
are averaged and the total beam current is esti- 
mated by integrating this average profile for an 
assumed axi-symmetric distribution about the line 
of maximum current density. The deflection angle 
value which best fits the data is then obtained by 
iterating assumed deflection angle values and 
solving for the minimum Integrated beam current. 

The final beam current obtained in the analysis Is 
then compared to the measured beam current to in- 
sure that they agree to within 5 to 10%. The di- 
vergence angle (a in Fig. 1) is computed as the 
half angle of the cone enclosing 95% of the inte- 
grated beam current determined In the aforementioned 
analysis. This definition of half angle is the same 
as the one used by Aston3.4 in the study of di- 
vergence characteristics of aligned grids. 

Initially, current density data were collected 
manually by reading a meter that could be switched 
to read each Faraday probe in the rake. These re- 
sults showed the linear relationship between the 
grid offset angle (e) and the deflection angle (b) 
suggested by Eq. 10 but, inconsistencies were 
observed when data from various runs were compared. 
It was determined that these inconsistencies were 
due in part at least to variations in Ion source 
operating conditions that occurred while the probe 
currents were being read. In order to minimize 
this drift the manual data collection was eliminated 
in favor of automated collection using a Hewlett 
Packard 3054 Data Logger. Results obtained using 
the data logger have been very reproducible and 
self-consistent. 

In order to Isolate the effects of various 
parametric changes to be introducted In the experi- 
ment a standard grid geometry was selected and then 
one grid parameter was varied while other parameters 
were held at their standard value. The standard 
case selected is defined by the following geometri- 
cal and operating conditions: 


Screen hole diameter (d $ } 

2,06 mm 

Total accelerating voltage (Vj) 

600 vol ts 

Net-to-total accelerating 
voltage ratio (R) 

0.70 

Screen grid thickness ratio 
(t s /d s ) 

0.18 

Accelerator grid thickness 
ratio (t a /d s ) 

0.37 

Accelerator grid hole diameter 
ratio (d a /d s ) 

0.64 

Grid separation ratio U g /d s ) 

0.49 

As suggested by the above list. 

the geometri- 


cal grid parameters of Fig. 1 are normalized using 
the screen hole diameter, which was held constant 
at 2.06 mm for all tests. The range of variation 
of the parameters investigated in this study were: 

Total accelerating voltage (Vt) 600 - 1100 

vol ts 

Net-to-total accelerating 

voltage ratio (R) 0.35 - 0.90 
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Screen grid thickness 
ratio (t 5 /d s ) 

Accelerator grid thickness 
ratio (t a /d $ ) 

Accelerator grid hole dia. 
ratio (d a /d s ) 

Grid separation ratio Ug/d s ) 


0.18 - 0.37 
0.12 - 0.74 


0.64 - 1 .00 


0.49 - 1.00 


In addition, the discharge voltage was main- 
tained at 40 v and the bell jar pressures were held 
between 6 - 8 x 10 -6 Torr. An argon flow rate of 
10 mA equivalent was used throughout the experiments 
because it was sufficiently high to produce stable 
discharge chamber operation while being low enough 
to keep the charge exchange ion production less 
than one percent of the beam current. 


Results 


CURRENT DENSITY 
PEAK 

I 
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A typical beam ion current density profile 
collected in the manner outlined in the previous 
section is shown in Fig. 3a. As this figure sug- 
gests, the profiles are relatively symmetrical. 

This symmetry is retained within the useful range of 
beam deflections. Using the information on this 
profile one can determine the beamlet divergence and 
deflection associated with this particular grid geo- 
metry and relative grid offset. A typical plot of 
deflection angle vs. grid offset angle is presented 
in Fig. 3b; this one for the standard grid set 
operating at a normalized perveance per hole (NP/H) 
of 1.02 x 10" 9 A/V 3 / 2 where 



Figure 3b shows the expected linear relationship 
between deflection angle (s) and grid offset ( 6 ) up 
to a grid offset angle of n, 15° where the beamlets 
begin to intercept the accelerator grid. The de- 
flection angle at the onset of this non-linear 
behavior is designated Bu (for B useful), as sug- 
gested in Fig. 3b. The data shown is for the 
standard case accelerator grid hole diameter ratio 
of 0.64. At accelerator hole ratios near 1.00 the 
beam deflection beyond B u levels off rather than 
increasing as shown in Fig. 3b. This change in the 
slopes above the point where impingement begins to 
increase is probably related to the observation that 
the beam profiles become increasingly asymmetrical 
at high grid offset angles. It should be noted, 
however, that the non-linear beam deflection always 
occurs at relatively high impingement levels which 
are outside of the normal thruster operating regime. 


o BEAM CURRENT DENSITY PROFILE 


DEFLECTION 
ANGLE [£](deg) 



Fig. 3. Typical Standard Grid Data. 

The coefficient in this equation is reasonably 
close to the theoretical constant of 14.3° 

(1/4 rad) and is also in reasonable agreement with 
previously measured values. 8 It is also noteworthy 
that Eq. (8) predicts an increase in this constant 
if the deceleration field E ? is significant as it 
should be for the data of Fig. 3 where the net-to- 
total accelerating voltage ratio (R) is less than 
unity. 


The experimental deflection characteristics 
illustrated in Fig. 3b can be compared with the 
theoretical relation expressed in Eq. 10. From 
Fig. 3b, a deflection angle (e) of 6° corresponds 
to a grid offset angle (a) of 15.5° within the 
linear region. Using Eq. 1 to determine e, one 
arrives at an experimental correlation of the same 
form as the theoretical value given in Eq. 10, 
namely 

B = - 15.8 y ~ (B in degrees) , (12) 

g 


Typical divergence and impingement character- 
istics for the grid set of Fig. 3 are given in 
Fig. 4. It can be seen in Fig. 4a that the di- 
vergence angle (a) does not vary significantly over 
the useful range of grid offset angles where the 
impingement currents are acceptably low and opera- 
tion would normally occur. With a larger accelera- 
tor grid hole diameter ratio (d a /d s ) the divergence 
angle drops slightly from the value at a significant 
offset to a minimum at zero offset. Even in this 
case, however, the total variation in divergence 
angle over the useful range of grid offset angle 
was not observed to exceed two degrees. The rapid 
decrease in divergence angle apparent in Fig. 4a at 
large grid offset angles is probably due to mter- 
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a DIVERGENCE ANGLE DATA 



Fig. 5. Effect of Grid Offset on Impingement 
Current. 


IMPINGEMENT CURRENT (fi A) 



GRID OFFSET ANGLe[s] (deg) 
b ACCELERATOR GRID IMPINGEMENT DATA 


Fig. 4. Typical Standard Grid Data. 

ceptlon of the most divergent ions from the beam- 
lets by the accelerator grid. In Fig. 4b this 
rapid increase in accelerator grid Impingement cur- 
rent is apparent and it defines the range of useful 
grid offset angles (s u ). In order to quantify this 
range it is defined here rather arbitrarily as the 
grid offset angle at which the Impingement current 
reaches 10% of the beam current. This sets the 
limit on the useful range of beam deflection (6 U ) 
which is indicated on Fig. 3b. 

The useful range of the grid offset angle (6 U ) 
is a function of both the geometrical parameters of 
the grids and the perveance at which they are opera- 
ting. Figure 5 illustrates the effect of perveance 
on typical impingement current vs. grid offset angle 
plots. The two values of useful range of grid off- 
set angle (6 U ) appropriate to the two levels of 
perveance are identified and they show that the 
limiting offset angle drops as the perveance Is In- 
creased. This presumably occurs because the beam- 
let more nearly fills the accelerator grid aperture 
at the higher perveance level. Using a collection 
of data like those shown in Figures 3 and 4, one 
can define a group of parameters that characterize 
the optical behavior of deflected Ion beamlets. The 
parameters that have been selected for use here are: 


1. The useful ion beamlet deflection 
range (p u ) 

2. The deflection-to-offset angle ratio 
in the useful range { B/6 ) 

3. Divergence angle In the useful range (a) 


Effect of Grid Offset on Beamlet Divergence 

The assumption that a single divergence angle 
characterizes the beamlets over the range of useful 
deflection (3 above) Is supported by the results of 
Fig. 6, which show that the divergence angle Is In- 
dependent of grid offset angle over the full range 
of perveance. The divergence angle results also 
agree to within <2° with those obtained by Aston 3 
using aligned grids. In comparing these results It 
should be noted that a slightly different definition 
of the effective Ion acceleration length Is used to 
compute the perveances. 9 In Aston's results the 
effective acceleration length Is defined as: 




Fig. 6. Effect of Grid Offset on Beamlet 
Divergence. 
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V areas in the present work the screen grid thick- 
ness is also incorporated in the definition of ac- 
celeration length as* 

l e = [f*g + V 2 + ^2.] 


Figure 7 shows the effects of net-to-total voltage 
ratio (R) on beamlet divergence (a). These typical 
curves again agree reasonably well with those ob- 
tained by Aston. 3 Because of this agreement di- 
vergence curves will not be presented here for each 
parameter investigated, rather the reader is re- 
ferred for divergence angle data to Aston's work 3 
which utilized fixed grids and Is therefore con- 
sidered to give more accurate divergence angles 
than the present study. 



that this results in a weaker electric field be- 
tween the grids. Linear optics theory suggests 
this will result In less beamlet deflection for a 
given grid displacement (the grid separation appears 
in the denominator of Eq. 8). Figure 9 shows that 
the ratio of the beamlet deflection angle to the 
grid offset angle is independent of the grid sepa- 
ration ratio over the full range of perveance in- 
vestigated. This is not surprising because one can 
combine a first order expansion of Eq. (1) andEq. (9) 
together and obtain the expression for this defiec- 
tion-to-grid offset ratio. It shows this ratio 
should indeed be independent of the grid separation 
ratio to first order for this case of a thin screen 
grid. The independence of the deflection angle-to- 
grid offset ratio also confirms its preference as 
a correlating factor over other factors such as the 
ratio of deflection angle to the linear offset dis- 
placement (B/c). 



Fig. 7. Effect of Net-to-Total Accelerating 
Voltage Ratio on Divergence. 


Fig. 8. Effect of Grid Separation Distance on 
Deflection Angle Range. 


Effect of Normalized Perveance Per Hole 


The effect of increasing the normalized per- 
veance at which the grid set is operating always has 
the effect of decreasing the usable deflection angle 
range. As perveance is increased an unacceptable 
impingement level is eventually observed even at the 
zero grid offset condition. If a moderate beamlet 
steering capability is desired, this suggests the 
grid-hole pair must be operating well below the 
maximum obtainable, al igned-grid perveance. Figure 6 
shows the shape of a typical divergence angle vs. 
perveance curve. One would generally want to 
operate at the minimum of this curve, but this cor- 
responds to a rather high perveance condition where 
the beamlet cannot be deflected substantially with- 
out inducing a high impingement. Consequently, 
although the divergence angle shows little change 
with grid offset, the beamlets would generally have 
to be operating at a perveance below the one for 
minimum divergence in order to facilitate operation 
at a moderate offset without high impingement. All 
of the experiments suggested that perveance had no 
significant effect on the deflection-to-offset 
angle ratio ( 6/5 ) . 

Effect of Grid Separation Distance 

The deflection angle range is observed in 
Fig. 8 to decrease as the separation distance be- 
tween the grids is increased. This can be explain- 
ed by noting that the total voltage remains con- 
stant as the separation distance is increased and 
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Fig. 9. Effect of Grid Separation Distance on 
Deflection-to-Offset Angle Ratio. 


Effect of Accelerator Hole Diameter 


Accelerator hole diameters larger than the 
standard value enable grid sets to operate over 
much larger deflection angle ranges, as shown in 
Fig. 10. This is presumably due to the fact that 
the beamlets have more room to deflect before being 
intercepted by the accelerator grid. Figure 10 al- 
so shows that larger accelerator holes will permit 
a given beamlet deflection at a much higher per- 
veance level; power supply limitations prevented 
operation at impingement limited perveance levels 
with the larger accelerator grid hole diameters 
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used In this study. As Fig. 11 shows, the accel- 
erator hole diameter had no significant effect on 
the deflectlon-to-offset angle ratio. 


I2 r 
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Fig. 10. Effect of Accelerator Grid Hole Diameter 
on Deflection Angle Range. 
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Fig. 12. Effect of Net-to-Total Accelerating 

Voltage Ratio on Deflection-to-Offset 
Angle Ratio. 



Fig. 13. Effect of Net-to-Total Accelerating 

Voltage Ratio on Deflection Angle Range. 


Effect of Accelerator Grid Thickness 


Fig. 11. Effect of Accelerator Grid Hole Diameter 
on Deflection-to-Offset Angle Ratio. 


Effect of Net-to-Total Accelerating Voltage Ratio 

The net-to-total voltage ratio (R) has a 
moderate effect on the deflection-to-offset angle 
ratio (6/6). As the R-value decreased from 0.90 
to 0.35, 6/6 increased from 0.4 to 0.6, as shown 
in Fig. 12. This Is explained by considering that 
at a constant total voltage, a decrease in the 
R-value is accompanied by a decrease in the net 
acceleration voltage and hence in the Ion velocity. 
The ion is then subjected to the lens effect of the 
accelerator aperture for a longer time and this re- 
sults In a larger deflection angle for a given grid 
displacement. The deflection angle ranges as de- 
picted in Fig. 13, show only a slight dependence on 
the R-value as long as It Is greater than about 0.5. 
When the R-value is reduced to about 0.35 however 
the deflection range becomes more restricted at the 
higher perveances. This effect is presumably due to 
the fact that beamlets which become more divergent 
at low R-values begin to intercept the accelerator 
grid at a lesser grid offset angle. 


As shown In Fig. 14, the accelerator grid 
thickness had no significant effect on the deflec- 
tion-to-offset angle ratio. However, In Fig. 15 the 
thick accelerator grid is seen to exhibit a pro- 
nounced reduction in deflection angle range. This 
is considered to be due to the fact that Inter- 
ception of the beamlet on the downstream edge of 
the accelerator grid occurs at a much smaller grid 
displacement with the thickest grid. It is note- 
worthy that a very thin accelerator grid showed no 
improvement In the deflection angle range over that 
observed for the standard case thickness. 
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Fig. 14. Effect of Accelerator Grid Thickness on 
Deflection-to-Offset Angle Ratio. 
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Fig. 15. Effect of Accelerator Grid Thickness on 
Deflection Angle Range. 

Effect of Screen Grid Thickness 

The effect of doubling the standard screen 
grid thickness on the deflection characteristics of 
beamlets is examined in Figs. 16 and 17. This 
change had only a slight effect on the beamlet de- 
flection characteristics. A slight reduction in 
deflection range was observed for the thicker screen 
grid as illustrated in Fig. 16. This is supported 
by previous work by Aston’ which shows that thinner 
screen grid geometries are capable of achieving 
slightly higher perveances. This is consistent with 
the fact that the thinner screen grid exhibits a 
slightly higher deflection angle range at a given 
perveance level. There appears to be a slight dif- 
ference in the deflection-to-offset ratio for dif- 
ferent screen grid thicknesses as shown in Fig. 17. 



NORMALIZED PERVEANCE PER HOLE (A/V V * > 


Fig. 16. Effect of Screen Grid Thickness on 
Deflection Angle Range. 



NORMALIZED PERVEANCE PER HOLE ( A / V '* ) 

Fig. 17. Effect of Screen Grid Thickness on 
Defl ection-to-Offset Angle Ratio. 



An inspection of Eq. 8, which pertains to 
linear optics theory, reveals that the total ac- 
celeration voltage does not appear so this para- 
meter would not be expected to affect the beamlet 
deflection parameters to first order. Figures 18 
and 19 do show however that these parameters were 
affected somewhat by a change in total accelerating 
voltage from 600 v to 1100 v. Figure 18 reveals a 
larger deflection angle range over the complete 
perveance range at the higher voltage. Figure 19 
shows that while the magnitudes of the deflection- 
to-offset angle ratio are comparable at the two 
voltages, changes in perveance seem to induce a 
greater change in this ratio at the higher voltage. 
It should be noted that some of the data on 
Figs. 18 and 19 Involved operation of the ion 
source at a slightly higher discharge voltage and 
magnet current. This did not appear to alter the 
data trends. It was done so the ion source dis- 
charge could be sustained at the higher beam cur- 
rents associated with operation at the 1100 v 
condition. 



NORMALIZED PERVEANCE PER HOLE (A/V**) 

Fig. 18. Effect of Total Accelerating Voltage on 
Deflection Angle Range. 
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Fig. 19. Effect of Total Accelerating Voltage on 
Deflection-to-Offset Angle Ratio. 


Dished-Grid Example 


In this section the beamlet vectoring data 
presented in the previous section of this paper 
will be used to examine the beamlet deflections that 
would be expected in a sample di shed-grid geometry. 
The objective is to determine the effects of such 
design factors as dish depth and beamlet steering on 
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the impingement characteristics of a dished-grid 
set. In order to conduct this analysis the current 
density profile and thermal distortion characteris- 
tics of a typical dished-grid system will be 
examined. 

The parameters used to define the particular 
spherically dished grid set considered here are 
shown in Fig. 20. The beamlet pattern shown sug- 
gests aligned screen and accel grid holes that 
would cause the beamlets to emerge in this divergent 
pattern. Using current density profile data for a 
typical ion source together with reasonable values 
for grid geometric parameters, the radial screen/ 
accel grid hole displacement profile needed to cor- 
rect this condition will be determined. At this 
operating condition the beamlets would be aligned 
with the thruster axis. In addition, the impinge- 
ment limited perveance will be calculated as a 
function of grid radius using the deflection angle 
range data present in this paper. This will enable 
one to determine where high impingement might be 
expected over the surface of the grids. 



Fig. 20. Dished Grid Geometry. 

The thruster selected was the 30 cm dia. 800 
series Engineering Model Thruster. The current per 
hole data required as input for analysis of the 
grids when this thruster is operating at a 2.0 A 
beam current are shown in Fig. 21. 10 These data 
were obtained using mercury propellant, so the cur- 
rent per hole cited in Fig. 21 had to be multiplied 
by the square root of the mercury to argon atomic 
mass ratio to obtain the allowable perveance data 
needed for comparison. In order to determine the 
grid separation distance and relative hole positions 
at operating temperature, the linear expansion model 
of Rawlin, Banks, and Byers^ 1*12 wa s applied to the 
molybdenum grids and support rings. This model may 
not be applicable to all grid assemblies, but it is 1 
useful here to illustrate the application of beam- i 
let vectoring data to the analysis of dished-grid 
accelerator systems. 

When grids are dished and then provided with 
the desired separation distance, the screen-accel- 
erator hole pairs become mis-aligned, and the 
beamlets are directed severely off-axis. This re- 
quires hole-pattern compensation, which for grids 



Fig. 21. Current per Hole Profile for Sample 
30 cm Thruster. 

dished convex downstream (Fig. 20), Implies an ac- 
celerator grid hole pattern with center-to-center 
hole spacings increased or a screen grid with hole 
spacings reduced. By using some compensation a 
grid-hole pair can be re-aligned to the zero offset 
angle condition, however, an additional amount will 
be required in order to cause each beamlet to be 
directed parallel to the thruster axis. In order 
to accomplish this the accelerator grid hole 
spacings are designed to be (1 + e a ) times the 
screen hole spacings. Figure 22 shows how the grid 
hole offset angles vary for a given grid dish depth 
as a function of radius under the influence of the 
indicated compensations (e a 's) and the thermal 
gradients that the grid thermal modell0»ll predicts. 
Also shown in Fig. 22 as a solid line is the grid 
offset required as a function of radius to produce 
para-axial beamlets. It is determined from a de- 
flection-to-offset plot (Fig. 12 for R = 0.70 in 
this case) using as input the angular separation 
between the thruster centerline and the normal to 



Fig. 22. Dished Grid Offset Angles for Varying 

Accelerator Grid Compensation, h = 2.25. 
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the accel grid surface at each radius for the 30 cm 
thruster (0's in Fig. 20). Figure 23 presents the 
same type of data as Fig. 22 but In this case the 
dish depth is less. The magnitude of the dish 
depth affects the value of grid offset needed to 
re-all gn beamlets with the thruster centerline. A 
comparison of Figs. 22 and 23 shows that the 
smaller dish depth requires much less grid offset 
for the outer radius holes, and therefore It would 
be expected that the outer holes could then be 
operated at higher perveance levels. The drawback 
to using a smaller dish depth is that this leads to 
greater variations in the screen hole to accelerator 
hole separation distance as a function of operating 
temperature. It is apparent from Figs. 22 and 23 
that a single value of accel grid compensation (e a ) 
does not produce the required grid offset over the 
entire grid surface for the thermal deformations 
assumed here. This situation worsens as the dish 
depth is reduced. 



optimum dish depth seems to fall approximately mid- 
way between h = 2.25 cm and h = 1.12 cm. At zero 
radius, the estimated perveance profile drops with 
decreasing dish depth. This is due to the decrease 
In separation gap which is a maximum at the center- 
line and is more pronounced at small dish depths 
(the screen grid is assumed to be hotter than the 
accelerator grid and therefore expands more). These 
problems could be alleviated somewhat by varying 
the initial grid separation and accelerator grid 
hole size as a function of thruster radius. 



Fig. 23. Dished Grid Offset Angles for Varying 

Accelerator Grid Compensation, h = 1.12. 

Figure 24 shows perveance limit design plots 
for the sample case being considered. The allowed 
perveance profiles in this figure were determined 
by computing the angular displacement between the 
thruster centerline and normal to the grid surface 
at each radial location (0's in Fig. 20). This in- 
formation is then used to enter the useful deflec- 
tion range plots (Fig. 10 for d a /d s = 1.00 In this 
case) in this paper to determine the maximum per- 
veance allowed for a given deflection angle at 
each radial location before direct accelerator grid 
impingement becomes excessive. Also shown in Fig. 24 
are the estimated perveance profiles that would be 
observed in the 30 cm thruster operating at the cur- 
rent per hole values of Fig. 21 with the dish depths 
indicated on Fig. 24 and the computed grid thermal 
distortions. It can be seen from Fig. 24 that the 
allowed profile for the shallowest dish depth 
(h = 1.12 cm) Is above the corresponding estimated 
profile at all radii. This indicates the impinge- 
ment levels will be low over the entire accelerator 
grid surface in this case. For the greater dish 
depth curves however the allowed profiles fall be- 
low the estimated ones at the greater radii . In 
this case impingement would be expected to be more 
substantial. For the case under consideration the 


Fig. 24. Typical Comparison of Allowed Perveance 
and Estimated Perveance for 30 cm Dished 
Grids of Varying Dish Depths 


A potential problem also exists for large 
diameter multipole thrusters and is worthy of men- 
tion. As seen in Figs. 22 and 23 realtively large 
offset angles may be required for the apertures at 
the outer edges of the grids on large diameter 
thrusters. Because of the flatter plasma density 
profiles in multipole thrusters, these outer radius 
holes may be expected to operate at higher perveance 
values and this would be expected to create even 
more severe impingement problems when the beamlets 
are re-directed parallel to the thruster axis than 
was observed in Fig. 24 for the 30 cm divergent 
field thruster. 

Conclusion 

Ion beamlet deflections of several degrees can 
be accomplished by displacing the screen and accel- 
erator grid holes relative to each other. This 
deflection is limited to a useful range by the on- 
set of high accelerator grid impingement currents. 
This range is primarily dependent on the perveance 
at which the grid is operating and the diameter of 
the accelerator grid aperture. The useful range is 
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also dependent to a lesser extent on the net-to- 12. Rawlin, V. K., Banks, B. A., and Byers, D. C., 

total voltage ratio, grid separation distance, grid “Dished Accelerator Grids on a 30-cm Ion 

thicknesses, and the total accelerating voltage. Thruster," Journal of Spacecraft and Rockets, 

Over the useful range of beam deflections the beam- Vol. 10, No. 1, Jan., 1973, pp. 29-35. 

let divergence angle is independent of deflection. 

In this range the deflection angle is also a linear 
function of the grid offset angle. The constant of 
proportionality between these angles exhibits a 
dependence on the net-to-total accelerating ratio 
but is relatively independent of the other para- 
meters studied. Good agreement between the theo- 
retical approximation and experimental results for 
this constant has been demonstrated. The beamlet 
deflection data developed in this study can be used 
to design large diameter dished-grid accelerator 
systems in which the ion beamlets are to be directed 
along the thruster axis. An important design para- 
meter that is crucial to the maintenance of low 
impingement currents in these designs is the grid 
dish depth. 
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Abstract 


The flexible magnetic field thruster is a 
unique research tool for studying the behavior of 
direct current electron-bombardment thrusters. It 
utilizes a long wire anode shielded by a circum- 
ferential magnetic field. The magnetic field 
around the anode, which is produced by passing a 
direct current of up to 150 amperes through the 
wire, restricts primary electron flow from the dis- 
charge plasma to the anode. Different magnetic 
field configurations can be created by routing the 
anode wire in various ways through the discharge 
chamber. Divergent, cusped and multipole magnetic 
field configurations have been created in this man- 
ner. The thruster is also designed so Ion currents 
to various internal surfaces can be measured directly 
and these measurements facilitate calculation of 
the distribution of ion currents within the dis- 
charge chamber. Experiments indicate that the 
distribution of ion currents inside the discharge 
chamber is strongly dependent on the shape and 
strength of the magnetic field but independent of 
the discharge current, discharge voltage, and neu- 
tral flow rate. Measurements of the energy cost 
per plasma ion Indicate that this cost decreases 
with increasing magnetic field strength due to in- 
creased anode shielding from the primary electrons. 
Energy costs per argon plasma ion as low as 50 eV 
were measured. The energy cost per beam ion was 
found to be a function of the energy cost per 
plasma ion, extracted ion fraction and discharge 
voltage. Part of the energy cost per beam ion goes 
into creating many ions in the plasma and then ex- 
tracting only a fraction of them into the beam. The 
rest of the energy goes into accelerating the re- 
maining plasma ions into the walls of the discharge 
chamber. Measurement of ion fluxes across a virtual 
anode surface appear to indicate that ions cross 
this surface with velocities approaching their 
random thermal velocity rather than the Bohra 
velocity. 

Introduction 

A great deal of effort has been expended over 
the last two decades on research to improve the 
efficiency of electron-bombardment ion thrusters. 
Much of this effort has centered on the containment 
of the high energy or primary electrons in the main 
discharge region of the thrusters, and has resulted 
in a rather interesting evolution of the thruster 
magnetic field configurations. ^ >2 In spite of this 
effort the energy cost per beam Ion (in electron 
volts) for most thrusters is still on the order of 
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15 to 30 times the first ionization potential en- 
ergy of the propellent atoms.* For example, the 
J-series thruster^ using mercury, operates with a 
discharge power of 190 eV/beam ion at a beam cur- 
rent of 2 amps and a 32 volt discharge. This is a 
factor of 18 greater than the first ionization po- 
tential for mercury. For SERT 11^ this factor is 
about 22. Ramsey's magnetoelectrostatic thrusters 5 
typically operate with factors of 16 for both argon 
and xenon as did Sovey's argon line-cusp thruster.® 
Multipole thrusters of Isaacson,' using xenon, and 
Longhurst,® using mercury, operated with discharge 
power losses around 18 and 26 times the respective 
Ionization potentials. 

Dugan and Sovie 5 have shown that the energy 
required to produce an argon ion in plasmas similar 
to that commonly found in ion thrusters should be 
on the order of 50 eV. That this energy is greater 
(by approximately a factor of 3) than the ionization 
potential is due primarily to the fact that the ion- 
ization process is accompanied by some excitation of 
the neutral atoms. Inherent in this analysis is the 
assumption that electrons lose energy through in- 
elastic collisions with the neutral atoms exclus- 
ively, i.e. that electron energy losses to the 
walls are negligible. In the case of ion thrusters, 
this condition is approximately satisfied when the 
primary electrons are prevented from reaching the 
anode by a magnetic field and when only low energy 
electrons are removed from the plasma at the anode. 
Thus, it should be expected that a properly designed 
Ion thruster (i.e., one which prevents primary elec- 
trons from having direct access to the anode) would 
produce argon ions in the plasma at an approximate 
energy cost of 50 eV. 

This indicates that the high energy cost per 
beam ion of typical thrusters is not the result of 
inefficient plasma ion production, but rather it is 
the result of the fact that the fraction of these 
ions being extracted into the beam is substantially 
less than unity. If every ion produced in the 
plasma was extracted into the beam, the thruster 
would be operating at a discharge power level of 
50 eV/beam ion (50 watts/beam ampere). This is 
probably the maximum energy efficiency that could 
be expected for an argon thruster. 

To approach this maximum efficiency one must 
understand the behavior of, and be able to control, 
the ion currents inside the discharge chamber. The 
flexible magnetic field thruster was developed to 
study the behavior of these ion currents. The 
unique features of this thruster are that it allows 
drastic changes in its magnetic field configuration 
to be accomplished quickly and easily, and that it 
is designed so ion currents can be measured to its 
various interior surfaces. Changes in magnetic 
field configuration are made possible through the 
use of a copper wire anode through which a direct 

* For operation at the "knee" of the discharge 
power-propellent utilization curve. 
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current on the order of 100 amps is passed. This 
current sets up a circumferential magnetic field 
around the anode. The anode wire may then be bent 
into any shape inside the discharge chamber to pro- 
duce the magnetic field configuration of interest. 
Axial, mildly divergent, strongly divergent, cusped 
and multipole magnetic field configurations have 
all been created with this thruster. 

This paper provides a detailed description of 
the flexible magnetic field thruster design and 
some of the results obtained with it. These results 
include the measurement of ion currents to all 
thruster cathode potential surfaces. From these 
data the total ion production rate, energy cost per 
plasma ion and extracted ion fraction were calcu- 
lated. In this study these performance parameters 
were monitored as magnetic field configuration and 
strength, discharge current, discharge voltage and 
neutral flow rate were varied. 


Apparatus and Procedure 


The primary feature of the flexible magnetic 
field thruster is the anode. It is made from a 
3.2 mm diameter copper tube through which a rela- 
tively large direct current is passed. In de- 
signing the thruster the magnitude of the currents 
needed to produce fields that will prevent primary 
electrons from reaching the anode surface were 
first computed to insure they would be reasonable. 
This was done by recognizing that the magnetic flux 
density (B) at a radius (r) from the centerline of 
a wire carrying a current (Jf) is given by 


where u 0 is the permeability of free space. It has 
been assumed in writing this equation that the 
radius (r) where the field is sought is greater 
than the wire radius (r Q ). In order to prevent 
collection of electrons having an energy E e (In eV) 
Isaacson found the following condition on the in- 
tegral of the magnetic flux density had to be sat- 
isfied in a similar geometrical situation 10 

r i K 

/ Bdr > 6.74 x 10 . (2) 


In this expression r! is either the radius at which 
the magnetic flux density drops to a negligible 
value or the radius at which electrons are injected 
into the magnetic field, whichever is less. Com- 
bining Eqs. 1 and 2 and performing the indicated 
integration one can solve for the current required 
through the wire. After substituting for the per- 
meability one obtains 


J 


f 


33.7 vf“ 
e 



(3) 


For 40 eV primary electrons, a 1 .6 inn radius wire 
and a radius of electron injection in the radius 
range of 1 to 2 cm the required magnetic field cur- 
rent given by Eq. 3 is about 100 A. This is large, 
but not unacceptably so. For a 1.6 im radius anode 
of the order of a meter in length further analysis 
suggests that copper can carry the required current 
with a small voltage drop over the one meter tube 


length 0.1 V) as long as the wire is maintained 
at a sufficiently low temperature (a. 80°C). Compu- 
tations suggest this temperature can be maintained 
by using commercially available 3.2 mm dia. copper 
tube and circulating water through it at a flow 
rate of a few tens of ml/mm. 


The basic stainless steel, sheet metal chamber 
used to house the anodes in this study is shown 
schematically without an anode in Figure Id. For 
all of the data presented here this chamber was 
17 cm in diameter and 12 cm long. The other 
sketches in Fig. 1 show perspective cutaway views 
of anode configurations examined in this study. 

The magnetic field shapes produced by these windings 
are also shown. The refractory wire cathode used 




b) CUSPED CONFIGURATION 



RESULTING FIELD SHAPE 




Fig. 1 Flexible Magnetic Field Thruster 
Configurations. 
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in the tests was positioned 2.5 cm upstream of the 
grids in the manner suggested in Fig. Id for all 
anode configurations. The cathode employed 0.25 mm 
dia. tungsten wire about 30 cm in length. This wire 
was heated using alternating currents in the range 
12 to 15 A from a center-tapped power supply. 

Argon propellant was used because of the ease with 
which flow rates could be set and maintained. Tests 
were conducted at flow rates of 740, 1270 and 
1580 mA equivalent. Thruster performance was mea- 
sured at discharge voltages through the range of 
30 V to 50 V above the cathode center tap potential. 
The discharge current was adjusted through the 
range 2 A to 6 A by controlling the alternating 
current through the refractory cathode. 

The grids used in the tests had a 67% open area 
screen and a 30% open area accelerator. They were 
operated at a total accelerating voltage of 1500 V 
and a net-to-total accelerating voltage ratio of 
0.67 with a cold grid spacing 0.75 mm. Beam neu- 
tralization was accomplished using a refractory 
wire cathode. All tests were conducted in the 1.2 m 
dia x 4.6 m long vacuum test facility at Colorado 
State University. 

To facilitate the measurement of ion currents 
to various thruster surfaces the following features 
were incorporated into the design of the discharge 
chamber. The thruster back surface, cylindrical 
side wall, screen grid and cathode were all elec- 
trically isolated from each other. The positive 
high voltage was connected directly to the cathode 
center tap. The back surface, side wall and screen 
grid were connected to this point through one ohm 
resistors and through a power supply that was used 
to bias these surfaces negative of cathode potential. 
During thruster operation, these surfaces could be 
biased to -70 V (relative to cathode potential). At 
this bias, electron collection on these surfaces 
was eliminated and the ion current to each surface 
could be determined by measuring the voltage drop 
through each resistor. It is noted that ion cur- 
rents going to the anode surface and recombining 
there are not counted when measurements are made in 
this way. Due to the difficulty of measuring the 
ion current to the anode surface it was decided 
to treat the ion-electron recombination there as 
if it were a volume recombination process. That is 
the energy lost as a result of ion neutralization 
at the anode surface was included in the overall 
average energy cost to produce a plasma ion. 

Further, the anode area exposed to the plasma was 
relatively small compared to the area of the ca- 
thode potential surfaces and was held constant for 
each of the configurations of Fig. 1. Thus the ion 
current to the anode could be safely neglected when 
measurements of the total ion production rate were 
made. 

Applying a large negative bias to the various 
discharge chamber surfaces did not significantly 
alter the discharge chamber operation. When the 
bias was applied, the current from. each surface 
rose from zero (at floating potential) to the value 
corresponding to the ion current being collected. 

As the negative bias on a surface was increased the 
discharge and beam currents would also increase. The 
discharge current increase was caused by the col lec- 
tion, at the anode, of additional electrons which had 
previously neutralized the ions at the walls by 
coming directly from the plasma. When the negative 
bias was applied these electrons could no longer 
neutralize the ions at the walls directly from the 


plasma, but rather they had to be collected by the 
anode and then be passed through the discharge power 
supply so they could reach the walls. The beam cur- 
rent increase associated with biasing discharge 
chamber surfaces is a result of the increase in dis- 
charge current brought on by the phenomenon j'ust 
described. In conducting the tests it was found 
that when the discharge current was restored to its 
original value the beam current would also return 
to its original value. This indicates that the 
thruster discharge chamber efficiency was the same, 
regardless of whether the back, side and screen 
surfaces were allowed to float, or were biased 
strongly negative of cathode potential — provided 
the discharge current was held constant. At first 
glance, this is a rather remarkable result. In the 
floating case the ions reaching the walls are neu- 
tralized there by high energy electrons from the 
plasma. In the other case the ions are neutralized 
by electrons from the discharge supply which are 
primarily low energy electrons collected by the 
anode and then pumped up to the discharge chamber 
surface potential through the anode power supply. 
That these two conditions should result in identi- 
cal discharge chamber efficiencies requires an 
explanation. 

With the walls floating the discharge current, 
Jp, is the sum of the cathode emission current, Jg, 
and the beam current, Jg, so that. 


J D = J E + J B * (4) 

In this case the ions are neutralized at thfe walls 
by high energy electrons emitted by the cathode. 
When the walls are biased sufficiently negative of 
the cathode potential so that no electrons from the 
plasma can reach them, the discharge current 
becomes, 


where is the magnitude of the ion current to the 
walls — a current that remains unchanged as the 
walls are biased negatively. Now, if the discharge 
currents of Eqs. (4) and (5) are to be made equal 
and if, at this condition the beam currents are 
equal, then the cathode emission current (Je) in 
Eq. (5) must be decreased by an amount equal to Jy. 
That is the rate at which high energy electrons 
are injected Into the plasma for the case where the 
walls are highly negative must be reduced by an 
amount exactly equal to the rate at which high en- 
ergy electrons were lost to the walls in the first 
case. This results in the density of high energy 
electrons being the same in both cases, and since 
everything else was held constant, the efficiency 
must be unchanged. 

By measuring the ion current to the thruster 
body and screen grid in the manner suggested above 
one is provided with a powerful tool for evaluating 
thruster performance. It is really an extension 
of the Langmuir probe technique proposed by Sovey 6 
and is essentially the same approach used by 
Siegfried in studying hollow cathodes. ^ In this 
case however the back surface, side wall and screen 
grid act as probes. The value of the measurement 
lies in the fact that the sum of the body and screen 
grid ion currents represents the current of ions 
produced in the chamber but not extracted through 
the grids. The sum of this current and the beam 
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cu-rent is the total Ion production rate (expressed 
as a current). Knowing this total ion production 
rate one can compute the fraction of ions produced 
that are extracted (beam current- to-total production 
rate ratio) and the energy cost of a plasma ion 
(discharge power-to-total production rate ratio). 

The first of these parameters Indicates the effic- 
iency of the discharge chamber in directing ions 
through the grids; the second the efficiency of dis- 
charge chamber electrons in producing ions. 


Results and Discussion 
Existing Thruster Designs 

In the first part of this study three thruster 
magnetic field configurations were examined; the 
divergent field, the single cusped field and the 
multipole configuration. Each of these configura- 
tions was easily created with the flexible field 
thruster. Figure 1 shows these configurations and 
how the anode wire must bend in order to create 
each particular magnetic field shape. In each of 
these three configurations the loop separation was 
approximately 2 cm and the total anode length ex- 
posed to the plasma was 280 cm. In addition, the 
volume and surface area of the region bounded by 
the anode was the same for each configuration. The 
effect of increasing the magnetic field strength 
for the divergent field thruster configuration is 
illustrated in Fig. 2. In this figure the fractions 
of total ion current produced (J p ) which go to; the 
back surface (Jback/J p )i the side wall (Jside/Jp). 
the screen grid ( Jscreen/Jp) . and the beam (Jg/Op) 
are plotted as functions of normalized magnetic 
field current through the anode wire. These data 
were taken at a discharge current of 2.0 A, a dis- 
charge voltage of 50 V and a neutral flow rate of 
1580 mAeq. However, it was found that the ion cur- 
rent distribution in the discharge chamber was In- 
dependent of discharge current, discharge voltage, 
and neutral flow rate over the ranges investigated. 
Because of this independence from discharge condi- 
tions, the curves shown in Fig. 2 can be considered 
generally applicable to the divergent magnetic field 
configuration. Examination of Fig. 2 indicates that 
the fraction of the total ion current extracted in 
the beam increases with increasing field current. 

In addition, the fractions of the total ion current 
going to the back surface and screen grid increase 
slightly while that to the side wall decreases. 

This is clear evidence that the magnetic field 
strength can influence the distribution of ion cur- 
rents inside the discharge chamber. The decrease 
in the fraction of ion current to the side wall is 
believed to be primarily the result of a radial 
plasma density distribution in the discharge chamber 
that becomes increasingly non-uniform as the field 
current is Increased. Most of the ion current going 
to the back surface is assumed to be going there 
through the "hole" in the magnetic field on the 
centerline of the thruster. 

The distribution of ion currents for the cusped 
field thruster (Fig. lb) configuration is shown in 
Fig. 3. Again the fractions of ion currents to each 
surface are plotted against the normalized magnetic 
field current. Here the same basic trends exhibited 
by the divergent field configuration are seen. 
However, in this case the redistribution of ion 
currents is less pronounced than in the previous 
case. This appears to indicate that the plasma 



Fig. 2 Ion Current Fractions to Various Thruster 
Surfaces vs. Field Current for the 
Divergent Field Configuration. 



NORMALIZED FIELD CURRENT f J C /J F 1 
L r r max J 

Fig. 3 Ion Current Fractions to Various Thruster 
Surfaces vs. Field Current for the Cusped 
Field Configuration. 

density remains more uniform at high field 
strengths in the cusped configuration than in the 
divergent field geometry. This plasma uniformity 
is also reflected in the respective ion beam pro- 
files which are more uniform in the cusped field 
configuration than they are for the divergent one.’ 
It is believed that an additional factor which 
causes the side wall current fraction to be higher 
for the cusped configuration than for the divergent 
one is the plasma leakage that occurs at the cusp. 

Finally the distribution of currents for the 
multipole configuration (Fig. lc) are shown in 
Fig. 4. Here, interestingly enough, the magnetic 
field strength has no effect on the gross distribu- 
tion of ion currents inside the discharge chamber. 
Further, the fraction of current to each surface 
was found to be exactly equal to the area of that 
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Fig. 4 Ion Current Fractions to Various Thruster 
Surfaces vs. Field Current for the Multi- 
pole Configuration. 


Fig. 5 Normalized Ion Current Directed Toward 
Grids vs. Field Current. 


surface divided by the total area of the primary 
electron region. This was true only at zero field 
current for the first two configurations. This 
appears to Indicate that the Ions move toward the 
walls with nearly equal probably in all directions in a 
multipole thruster. It should be noted here that 
this result applies only to low field strength 
multipole thrusters. It was not possible to create 
the 0.3 Tesla magnetic fields now found at the pole 
pieces of some multipole designs. 12 The improved 
performance of these designs suggests that some re- 
direction of ion currents is indeed occurring. In 
any case, the uniform ion current distribution of 
low field strength multipole thrusters is clearly 
demonstrated here. The fact that ions are lost 
uniformly in all directions from the plasma implies 
that the extracted ion fraction should always be 
less than 50% for multipole thrusters of low field 
strength design. 


Substituting Eq. (7) into Eq. (6) yields, 
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[ J D-( J B + J M )] V D 
J B + J W 


This can be rewritten as, 
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Recognizing that (Jp-Jg) Vp/Jg is the discharge 
power Pp (in eV/ beam ion), and letting f be the ex- 
tracted ion fraction defined by, 


f = 


J B +J W 


(10) 


For comparison purposes, the fraction of ions 
directed toward the grids, (Jp + Jscreen)/Jp» ls 
plotted against the magnetic field current in Fig. 5 
for the three thruster configurations. The advan- 
tage of the divergent field geometry in terms of the 
fraction of ions directed toward the grids is clear. 

Plasma Ion Production 


The energy required to create an ion in the 
plasma has, as one would expect, a strong influence 
on the overall thruster efficiency. The relation- 
ship between this energy cost per plasma ion and 
that per beam ion can be related using quantities 
measured in the tests described in this paper. The 
energy cost per plasma ion (Pp-in electron volts) 
is defined by 

. < J D - V V D 

P P = Jn 


(6) 


where Vp is the discharge voltage. The ion produc- 
tion rate (jp) can be written as the sum of the 
beam current and ion current to the discharge 
chamber walls ( Jy) , i.e.. 


+ di.r 


(7) 


then, Eq. (9) can be written in the form. 
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This is the desired relation between the energy cost 
per plasma ion and the energy cost per beam ion. The 
first term on the right hand side of this equation 
reflects the energy lost because ions are produced 
that recombine on the walls and are not extracted 
into the beam. The second term represents the 
energy used to accelerate the plasma ions that go 
to the walls of the discharge chamber into these 
surfaces. 

The plasma ion energy cost has been determined 
for each thruster configuration and at each opera- 
ting point considered in this study using Eq. (6). 

A sample of the results obtained with the cusped 
field thruster are shown in Fig. 6. These data 
indicate that the energy cost per plasma ion de- 
creases with increasing magnetic field strength 
(due to increased anode shielding from primary 
electrons). In addition, it indicates that the 
cost increases with discharge current and decreases 
with neutral flow rate. However, it is relatively 
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Fig. 6 Plasma Ion Energy Cost in the Cusped Field 
Configuration. 

independent of discharge voltage. While these data 
were obtained with the cusped field configuration 
results obtained with the other thrusters showed 
similar trends. These and other data indicate that 
the energy cost per plasma ion can be made as low 
as 50 to 60 eV for the argon propellant used. 

In the course of measuring the total ion pro- 
duction rate it was noticed that this rate was 
frequently much greater than the neutral flow rate 
in equivalent mill lamperes. An example of this is 
shown in Fig. 7 for data collected on the cusped 
field thruster. The ion production rate being 
greater than the neutral flow rate is an indication 
that the walls of the discharge chamber behave as a 
virtual source of propellent atoms. Also shown on 
this figure is the beam current. This figure indi- 
cates rather dramatically how the ion production 
rate can be significantly greater than the neutral 
flow rate and yet the beam current is still signifi- 
cantly less than the flow rate. This does not mean, 
however, that if every ion produced could be ex- 
tracted into the beam that one would get a beam cur- 
rent greater than the flow rate. Instead, the 
neutral density would decrease causing a correspond- 
ing decrease in the ion production rate such that 
the beam current (assuming only singly charged ions) 
would always be less than the neutral flow rate. 

Other Considerations 

In addition to the tests described above, 
several other configurations of the flexible field 
thruster were tested in an effort to reduce ion 
losses to the walls still further. The number and 
type of configurations that can be tested is limited 
only by the imagination of the researcher. To date, 
no configuration of the flexible field thruster has 
directed a greater fraction of ions toward the grids 
than the divergent field configuration. However, 



Fig. 7 Comparison of ion Production Rate and Beam 
Current to the Neutral Flow Rate for the 
Cusped Field Configuration. 

in the course of this effort a very interesting 
phenomenon was observed. Tests were being run on the 
flexible field thruster in an axial magnetic field 
configuration as shown in Fig. 8. In addition to 
the solenoidal anode the thruster was equipped with 
an external solenoid. This allowed the magnetic 
field in the discharge chamber to be created by 1) 
passing a current through the anode wire only, 2) 
passing a current through the external solenoid 
only, or 3) passing currents in opposing or adding 
directions through both wires simultaneously to 
produce a wide range of magnetic field conditions. 
With this thruster configuration it was observed 
that at any given flow rate, discharge current and 
discharge voltage, the ion current to the side wall 
was smaller when the magnetic field was created 
with the external solenoid alone than it was when a 
magnetic field of identical strength was created 
with the internal solenoid alone. This indicates 
that the plasma is being contained better when the 
magnetic field is created by the external solenoid 
than it is when the internal anode wire is used. 

Why this should occur was not immediately apparent. 

The reason the external coil gives better per- 
formance than the internal one can be understood, 
however, by considering Fig. 8 further. Figure 8a 
shows the magnetic field lines for the case where 
the field is produced entirely by passing a current 
through the anode wire. Figure 8b shows the field 
lines when the field is produced entirely by the 
external solenoid. The anode wire in this second 
case serves only as the anode. In Fig. 8b we 
notice that certain field lines intercept the anode 
surface directly. As we progress outward from the 
centerline of this thruster along a radius, the sur- 
face of revolution of the first field line that we 
encounter which intercepts the anode surface is 
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Fig. 8 Axial Magnetic Field Thruster Configurations. 

called the virtual anode surface. It is so called 
because once an electron becomes attached to a field 
line on this surface it has a high probably of being 
collected by the anode. Consequently, few electrons 
would be expected to be found at radial locations 
greater than that of the virtual anode surface. 

This in turn implies that the plasma density should 
be reduced substantially in the region outside of 
this surface. It appears that the reduced Ion flux 
to the side wall is a result of the reduced plasma 
density in this region. For the case of Fig. 8a 
(magnetic field created by the internal solenoid), 
however, it is seen that no field lines intercept 
the anode directly. Instead they tend to form con- 
centric circles around the anode wire. Consequently 
there is no virtual anode surface. Further, the 
field strength between adjacent turns of the anode 
wire goes to zero. These two factors apparently 
enable the plasma to leak out of the volume bounded 
by the anode wire and the result is an ion current 
to the side wall greater than that observed with the 
external solenoid. 

When the virtual anode surface is present, as 
mentioned earlier, ion losses to the side wall can 
be greatly reduced. It was postulated that this 
reduced ion loss to the wall was the result of a 
greatly reduced plasma density adjacent to it. 
Further, this reduction in plasma density was be- 


lieved to occur across the virtual anode surface 
because of the removal of electrons from the plasma 
by the anode. If this Is indeed the case, that Is, 
if the plasma density is reduced dramatically 
across the virtual anode surface, then the ions 
must cross this surface with an average velocity 
substantially less than the Bohml3 velocity. This 
must be so for if the ions crossed this surface 
from a region of high density to a region of low 
density with the Bohm velocity, then an excess of 
positive charge would very quickly be built up in 
the low density region. This excess charge would 
serve to reduce the ion velocity across the virtual 
anode surface to the point where the flux of posi- 
tive charges across the surface was consistant 
with the condition of quasi -charge neutrality in 
both the high and low density regions. 

To see if this was indeed the case the follow- 
ing experiment was performed. A thruster having 
the configuration shown schematically in Fig. 9 was 
constructed. The anode consisted of a washer-shaped 
piece of 0.025 cm thick stainless steel. An axial 
magnetic field was created in the discharge chamber 
through the use of a six turn solenoid wrapped 
around the cylindrical thruster body. The cathode 
consisted of a tungsten wire bent into a circle 
10.2 cm in diameter. With this configuration the 
entire side wall was protected by the virtual anode 
surface but the physical anode area exposed to the 
plasma was extremely small. This allowed the ion 
current to the side wall to be measured very accu- 
rately since the ion current to the physical anode 
was negligible. The actual ion currents to each 
thruster surface are shown in Fig. 10 for an argon 
flow rate (m) of 1270 mA, a discharge current (Jp) 
of 2.0 A and a discharge voltage (Vp) of 50 V. It 
is noteworthy that this discharge chamber gave a 
higher extracted ion fraction than the divergent 
field design (Fig. la and Fig. 5). The remarkable 
aspect of the results of Fig. 10 is however that as 
the ion current toward the grids (Jp + Jscreen) ln_ 
creases by approximately a factor of 5, due to an 
increase in the magnetic field current from zero to 
100 A, the ion current to the side wall actually de- 
creases. This shows that the virtual anode surface 
does a remarkable job of containing the plasma. It 
is believed that this occurs because ions cross this 
surface with velocities approaching their random 
thermal velocity rather than the Bohm velocity. The 
following rough calculation supports this hypothesis. 

Assuming that the ions reach the accelerator 



Fig. 9 Axial Magnetic Field Thruster with Washer- 
Shaped Anode. 
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Fig. 10 Ion Currents to Various Thruster Surfaces 
for the Axial Field Thruster with Washer- 
Shaped Anode. 

system with the Bohm velocity based on an estimated 
electron temperature of 5 eV, the plasma density 
just inside of the virtual anode surface was esti- 
mated from an ion beam current density profile to be 
2.4 x 10- 16 m -3 . If the Ions crossed the virtual 
anode surface with the Bohm velocity, then the cur- 
rent to the side wall would be given by the formula. 


^side 
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where n is the plasma density(just cited), e is the 
charge for singly charged ions in coulombs, A^ is 
the area of the virtual anode surface and vg is the 
Bohm velocity corresponding to an electron tempera- 
ture (T e ) of 5 eV. If, on the other hand the ion 
current depends on the ion random thermal velocity, 
then 
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their thermal velocity rather than their Bohm 
velocity. 

This phenomenon can also be used to explain 
the low ion flux across the virtual anode surface 
observed in the SERT II thruster. 14 In this case 
the average current density across this surface was 
measured and found to be approximately 15% of the 
average current density toward the grids. Rough 
calculations suggest that if the ions have the Bohm 
velocity at this surface, then the average current 
density would be on the order of 250% of the average 
current toward the grids. If they are assumed to 
have the random thermal velocity however the calcu- 
lated current density is about 10% of the current 
toward the grids. It appears that the ions do not 
cross the virtual anode surface with the Bohm vel- 
ocity. It should be pointed out that these results 
are still preliminary and final judgement must await 
detailed probing of the plasma around the virtual 
anode surface. 

In light of this ion velocity phenomenon it 
should be noted that the flexible magnetic thruster 
does not produce any virtual anode surfaces regard- 
less of its configuration (for the reasons discussed 
with regard to Fig. 8a). This implies that the 
flexible field thruster does not imitate exactly 
those thruster designs which do create virtual 
anode surfaces. This includes virtually all cur- 
rent designs. However, based on results obtained 
with the thruster configurations shown in Figs. 

8a and 8b, it is believed that the ion flux dis- 
tributions obtained with the flexible field thrus- 
ter are still qualitatively correct, and meaningful 
comparisons can be made between different magnetic 
field configurations using this technique. 

Cone! uslons 


The ability to make drastic changes in magnetic 
field configurations, together with the ability to 
measure ion currents to different thruster surfaces 
makes the flexible magnetic field thruster a power- 
ful tool for the study of magnetic field configura- 
tion effects on thruster operation. With this 
thruster it was observed that the distribution of 
ion currents inside the discharge chamber was 
strongly dependent on the shape and strength of the 
magnetic field. This distribution of ion currents, 
however, was found to be independent of the dis- 
charge current, discharge voltage and the neutral 
flow rate. With no magnetic field the ions move to 
the thruster walls with equal probability in all 
directions, and the fraction of current to each sur- 
face bounding the plasma is equal to the area of 
that surface divided by the total surface area of 
the primary electron region. 


and T} is the ion temperature (assumed to be equal 
to 600°K) , k is Boltzmann's constant and m^ is the 
ion mass. Using the density cited previously, 

Eq. (12) gives a ion current to the side wall of 
1010 mA, whereas Eq. (13) gives a current of 40 mA. 
The measured current at the maximum magnetic field 
was 113 mA. If the one-fourth factor in Eq. (13) 
is eliminated, then the equation would correspond 
to the current due to a directed thermal velocity 
and the result is a current of 160 mA. Admittedly 
these calculations are somewhat crude but they are 
believed to be accurate to better than a factor of 
two. On the basis of these results it is suggested 
that the ions cross the virtual anode surface with 


The fraction of ion current to each thruster 
surface was found to vary as a function of magnetic 
field strength for both the divergent and cusped 
field configurations. Generally the fraction of 
Ion current to the side wall would decrease and the 
fraction of ion current toward the grids would in- 
crease as the field strength was increased for these 
configurations. Ho such redistribution of ion cur- 
rents was observed for the multipole configuration. 
In this case, the ion current distribution remained 
the same as it was for the zero magnetic field 
condition. 

Measurements of the energy cost per plasma ion 
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indicated that this cost decreased with increasing 
magnetic field strengths. This effect was attrib- 
uted to increased shielding of the anode from the 
primary electrons. Energy costs per argon plasma 
ion as low as 50 eV were measured. 

The energy cost per beam ion was found to be 
a function of the energy cost per plasma ion, the 
extracted ion fraction, and the discharge voltage. 
Part of the energy cost per beam ion goes into 
creating many ions in the plasma when only a frac- 
tion of them is extracted into the beam. The rest 
of the energy goes into accelerating the remaining 
plasma ions into the walls of the discharge chamber. 

The flexible field thruster does not exactly 
imitate the discharge chambers of existing thruster 
designs due to its inability to produce virtual 
anode surfaces. However, the Ion current distribu- 
tions measured with it are still believed to be 
qualitatively correct. The virtual anode surface 
provides a boundary between high and low density 
plasma regions. Ions from the high density region 
cross this boundary with velocities approaching 
their random thermal velocity. This indicates that 
ion losses from typical ion thruster discharge 
chambers can be reduced by bounding the plasma with 
virtual anode surfaces. 
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Abstract 


An experimental investigation into the ion ex- 
traction capabilities of accelerator systems with 
small screen hole diameters (<_ 2.0 mm) at net- 
accelerating voltages of 100, 300, and 500 V is 
described. The apparatus developed for this study 
is shown to be well suited for measuring the lm- 
pingement-1 lmited perveance, electron backstreaming, 
and electrical breakdown characteristics of two-grid 
ion accelerator systems. Results suggest that the 
impingement-limited perveance is not dramatically 
affected by reductions in screen hole diameter to 
1.0 mm, but impingement-limited performance was 
found to be dependent on the grid separation dis- 
tance, the discharge-to-total accelerating voltage 
ratio, and the net-to-total accelerating voltage 
ratio. Results obtained with small hole diameters 
and closely spaced grids suggest a new mode of grid 
operation where high current density operation can 
be achieved with a specified net acceleration volt- 
age. In this mode the grids are operated at a high 
net-to-total acceleration voltage ratio rather than 
the low ones that have heretofore been considered to 
yield higher current densities. Under these con- 
ditions it appears to be the onset of electron 
backstreaming that ultimately determines the ex- 
tractable ion current density. Beam current dens- 
ities as high as 25 mA/cm 2 were obtained using grids 
with 1.0 mm diameter holes operating at a net ac- 
celerating voltage of 500 V. 

Introduction 

In order to reduce the number of electrostatic 
ion thrusters required for a space propulsion mis- 
sion and hence the cost and mass of the thruster 
subsystem, it is generally desirable to increase 
the thrust densities of these devices above those 
attainable with currently available ion thrusters. 

To keep the thrusters operating at their optimum 
specific impulse this should be accomplished by in- 
creasing the ion beam current density extracted 
from the ion optics system without increasing the 
exiting ion velocity. This can be done under the 
constraint of the space-charge effect by reducing 
the grid spacings and hole diameters below the 
values used on present-day accelerator systems. The 
need for this type of accelerator system becomes 
even more apparent when one considers the current 
trend away from heavy propellants such as mercury 
to lighter propellants such as argon and xenon. To 
achieve the optimum specific Impulse for a typical 
mission, these lighter propellants require a lower- 
ing of the net accelerating voltages below those 
+ Work supported under NASA Grant NGR-06-002-11 2 

* Research Assistant, Department of Mechanical 
Engineering 

** Professor, Department of Mechanical Engineering, 
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being used for current accelerator system designs. 

In terms of the extracted current densities this 
means that the trend to lighter propellants neces- 
sitates the use of accelerator systems with smaller 
grid holes and closer grid separations in order to 
maintain current densities at their present levels. 

Although existing experimental performance 
data for accelerator systems with relatively large 
screen hole diameters (>_ 2 nm) is quite exten- 
sive,^ 2 * 3 this same information for smaller screen 
hole diameters is somewhat limited. Previous ex- 
perimental results have suggested that there might 
be a reduction in the ion extraction capabilities 
of thrusters utilizing grids with small screen hole 
diameters (< 2.0 mm).'» 2 » 3 Because of this possible 
degradation due to hole size, it was felt that be- 
fore any future high thrust density designs could 
be attempted an improved understanding of the ion 
extraction capabilities of grids utilizing these 
small diameter holes (< 2.0 mm) was needed. It is 
for this purpose that the present detailed investi- 
gation was undertaken. 

Background and Theory 

The essential components of a basic two grid 
ion accelerator system are shown in Fig. la. Al- 
though the side view of a single aperture is de- 
picted, it should be noted that a multi -aperture 
system is implied. The accelerator system is 
placed adjacent to a discharge chamber containing a 
low density ionized plasma at a potential of a few 
tens of volts above that of the screen grid. Around 
each screen hole a plasma sheath is formed which de- 
fines the boundary between the discharge plasma and 
the acceleration region. Potentials applied to the 
grids produce an electric field that accelerates 
the ions and repels the electrons coming into the 
sheath from the discharge plasma. Figure lb Illus- 
trates the variation of these electrical potentials 
associated with the grid geometry of Fig. la. The 
solid line of Fig. lb shows how the potential drops 
slightly from the plasma to the relatively high 
positive potential in the grid. It remains con- 
stant through the screen grid area then drops almost 
linearly in the space between the grids. It is 
constant through the accel grid and then rises to 
the potential far downstream of the accelerator 
grid. In this downstream region electrons are 
available to maintain charge neutrality In the 
plasma associated with the Ion beam. The dotted 
line of Fig. lb shows how for a normally operating 
grid the potential varies with position through the 
center of the grid holes rather than through the 
grids themselves. The difference between the solid 
and dotted lines is induced by the space charge 
effects of the ions passing through the grids. 

Table 1 defines and names the symbols shown in 
Figs, la and 1 b. 
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Fig. 1 Two-Grid Acceleration System 


Table 1 

Accelerator System Nomenclature 
t = Accelerator grid thickness 

a 

t = Screen grid thickness 
d a = Accelerator hole diameter 

a 

d g = Screen hole diameter 

Jig = Grid separation distance 

Ji e = Effective acceleration length 

V f j = Net accelerating voltage 

V-|. = Total accelerating voltage 

R = Net-to-total accelerating voltage ratio 


The current density extraction capability of 
an accelerator system aperture like the one shown 
in Fig. la is generally limited by space charge 
effects to the value given by Child's Law. 
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In this equation, representative of a one-dimension- 
al model of the Ion acceleration process, e 0 is the 
permittivity of free space, q is the ion charge, m-j 
is its mass and l represents the length through 


which the Ions are accelerated by the applied total 
accelerating voltage (Vj). In the case of the 
parallel plate vacuum diode, for which this equation 
was originally derived, i refers to the physical 
distance between the parallel plates. However, in 
the case of the extraction aperture the emissive 
surface is not defined by a plate whose position 
is fixed, but rather by a plasma sheath whose po- 
sition is dependent on the Ion current being ex- 
tracted. The current density of ions arriving at 
the plasma sheath, from the upstream plasma is de- 
termined by the Bohm condition for a stable sheath. 



In this equation the discharge plasma ion density 
(ni) is equal to the associated electron density 
(n e ) and the electron temperature (T e ) of the 
plasma is in electron volts. The quantity "e" is 
the electron charge. 

Experimental results obtained by Aston showed 
that increasing the Bohm current density, by in- 
creasing the plasma density caused the sheath to 
move closer to the screen hole.^ These results 
lead one to the conclusion that the acceleration 
length l adjusts itself so that Child's Law is al- 
ways satisfied. In normal operation the shape of 
the plasma sheath is similar to that shown in 
Fig. la. As the plasma density and as a result the 
ion current density Increase, the sheath moves 
closer to the screen hole and begins to flatten out. 
Eventually, an operating point is reached where the 
sheath can no longer focus the Ions through the ac- 
celerator grid aperture and direct ion Impingement 
upon the accelerator grid begins to occur. At this 
point the acceleration length t required by Child's 
Law has reached a value less than a physical accel- 
eration length required to produce a well focused 
ion beam. This physical acceleration length is 
dictated by the geometry of the accelerator system. 
In past theoretical and experimental investigations, 
the physical acceleration length used to correlate 
the performance of various grid geometries is the 
effective acceleration length i e where 

, d 2 . 

2 e = l>g + r F * (3) 

Replacing i in Eo. 1 with Ji e , and setting the beam 
current density (j) equal tothebeamlet current (Jg) 
divided by the area of a screen grid hole, Eq. 1 can 
be written 
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The left side of this expression is known as the 
normalized perveance per hole and it will be used 
in this paper to correlate the performance of ac- 
celerator grid systems as a function of the various 
grid geometries and the operating conditions in- 
vestigated. Of particular interest will be the 
impingement-limited normalized perveance per hole 
defined as the perveance per hole at the onset of 
excessive accelerator grid ion impingement. It 
should be noted that the right-hand side of Eq. 4 
has a fixed value for a particular propellant. This 
value is indicative of the perveance where one 
would expect to encounter the impingement limit. 


no 



In the case of argon, the propellant used in this 
study, this limiting value Is computed to be 
6.8 x 10 -9 A/V 3 / 2 . Assuming the grids can be held 
parallel at any separation distance there are two 
additional phenomena that limit the operating range 
of an ion acceleration system. These are electrical 
breakdown through the region between the grids and 
electron backstreaming into the discharge plasma 
from the plasma downstream of the grids. In prac- 
tice the occurrence of these phenomena constrains 
the absolute lower limits of grid separation 
distances that can be realized and thus limits the 
ion current which can be extracted from the ac- 
celerator system. Electrical breakdown occurs when 
the electrical field strength between the grids 
(total accelerating voltage minus the discharge 
voltage all divided by the grid separation distance) 
exceeds the breakdown limit. When this occurs ex- 
cessive electron currents begin to flow between the 
grids. This limit is determined by the geometry of 
the grids as well as such factors as the grid sur- 
face finish and the density of charged and neutral 
particles between the grids. 

Electron backstreaming occurs when electrons 
in the downstream ion beam plasma find a path along 
which they can be accelerated through grid apertures 
directly into the discharge chamber. This condition 
occurs when the potential variation through the 
apertures is similar to that shown by the centerline 
in Fig. lb. When this type of potential profile 
exists, the electrons in the beam plasma are no 
longer repulsed by the negative potential on the 
accelerator grid and they are therefore able to 
backstream into the discharge chamber. 

Many factors influence the onset of backstream- 
ing, and traditionally those factors pertaining to 
the accelerator system have been lumped together to 
form a dimensionless parameter known as the elec- 
tron backstreaming parameter. 5 The parameter used 
herein is 

B = (1 - W (df)( exp [d“] )* 

a a 

This parameter is slightly different from the one 
used in previous theoretical investigations where 
instead of (tg/d a ) the factor (z c ,/d a ) was used. 5 
This change was made because the factor (l„/ d.) 
gave a better correlation of the electron Back- 
streaming data collected in this study. 

Apparatus 

All testing has been conducted on a 8 cm dia- 
meter by 10 cm long ion source with a mildly di- 
vergent magnetic field. The magnetic field was 
derived from a long solenoidal winding wrapped 
around the outside of the thruster body. A cylind- 
rical copper anode was used and apart from this, 
nonmagnetic stainless steel construction was 
employed throughout the source. Refractory tungs- 
ten filaments were used for both the neutralizer 
and cathode emitters. The emission levels for both 
of these was controlled by adjusting the current 
flowing through the filaments from their ac power 
supplies. This source design was selected because 
it provided the high current density capability 
needed to assure a proper test of the grids, and it 
provided a stable discharge. To minimize the dis- 
charge chamber noise further, the voltage drop 
across the refractory cathode was reduced by using 
two short filaments in parallel. It is believed 


that this reduces the noise because it reduces the 
range of primary electron energies coming off of the 
cathode. 

The ion current extraction capability of var- 
ious nineteen aperture grid sets investigated in 
this study were each characterized by a different 
screen aperture diameter. The grid aperture pattern 
used was a hexagonal array with a center-to-center 
hole spacing of 2.5 mm for all aperture diameters. 
The grids were made from thin sheets of commercially 
available graphite, and to help ensure proper grid 
hole alignment, the screen and accelerator grids 
were simultaneously match drilled. As a result of 
this the screen and accelerator grids both had the 
same diameter apertures. The thicknesses of the 
screen and accelerator grids were both the same in 
each test and as hole size was varied for the 
different tests the grid thicknesses were also 
varied to hold the ratio of hole diameter to grid 
thickness constant at 0.25. 

Because the grid separation distance was a 
very important parameter in this study, the test 
apparatus was designed so it could be varied con- 
tinuously while the ion source was being operated. 
This is accompl ished, using the apparatus shown 
schematically in Fig. 2, by moving the fork-spaced 
wedge in the manner suggested by the large arrow. 

As this wedge is moved toward the grids It forces 
the accelerator and screen grid support plates 
apart and hence the screen and accelerator grids 
attached to these plates are also separated. All 
materials used in the construction of this appa- 
ratus are capable of withstanding high temperatures. 
Graphite was used for both the grids and grid sup- 
port plates to minimize the differential thermal 



Fig. 2 Variable Grid Spacing Apparatus 
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expansion between these components that could cause 
grid warpage. The distance over which the grids 
were supported was also minimized to minimize the 
effects of grid warpage. Initial alignment of the 
grids was accomplished manually by positioning and 
clamping each grid and then checking the alignment 
visually beneath a large illuminated magnifying 
glass. This alignment is maintained during opera- 
tion by the boron nitride guide posts shown in 
Fig. 2. Electrical isolation is assured by the 
sheet of Isomica between the screen grid support 
plate and the wedge. The wedge is mechanically 
connected to a micrometer used to measure and adjust 
Its translation from the outside of the vacuum 
system. 

Immediately after each period of data collec- 
tion, while the source was still hot, the grid 
separation distance was reduced until the measured 
electrical resistance between the grids went essen- 
tially to zero. This reference point on the micro- 
meter (assumed to correspond to zero separation) 
along with the geometry of the wedge could then be 
used to calculate the grid spacing from the micro- 
meter readings recorded during the testing. This 
procedure for defining the reference point was 
usually repeated two or three times for each test 
setup to determine if the reference point had 
changed for some unknown reason. Based on repeated 
observations. It is believed that the grid separa- 
tion distance was known to within f 0.01 mm. 

Test Procedure 


Because the objective of this study was to de- 
fine the limitations of the accelerator system and 
not those limitations that might be source related, 
every effort was made to operate the ion source so 
the highest possible impingement-limited perveance 
would be realized with each grid configuration. 

This was accomplished by optimizing the propellant 
flow rate, discharge voltage (Vg) and magentlc 
field strength for each grid configuration tested. 

In general, it was found that discharge voltage and 
magnetic field strength could be maintained at a 
single optimum value for a particular net accelera- 
ting voltage. However, the propellant flow rate 
was more critical and had to be optimized for every 
grid separation investigated. This optimum flow 
rate was almost always equal to the minimum flow 
rate at which the source operated stably over the 
full range of beam currents Investigated. The con- 
duct of the tests involved the establishment of 
this optimum flow with a particular grid set, grid 
spacing and net accelerating voltage and then in- 
creasing the cathode emission current in steps to 
Induce increases in the beam and impingement cur- 
rents. These currents were recorded manually from 
high accuracy digital gauges at each step. This 
process continued until the impingement current 
reached a high value (=10% of the beam current). At 
this time, the cathode emission was reduced and the 
procedure for collecting electron backstreaming 
data was executed. 

The collection of backstreaming data was ac- 
complished by reducing the accelerator grid poten- 
tial (Vt-Vn) while watching for the rather sudden 
increase in beam current that is indicative of elec- 
tron backstreaming into the discharge chamber. This 
procedure yielded a minimum accelerator grid voltage 
which could be used along with the net accelerating 
voltage { Vjg) to calculate the upper limit on the 
net-to-total accelerating ratio (Rmax) which could 


not be exceeded without inducing electron back- 
streaming. Since the grid separation and grid 
geometry were known, the value of the electron back- 
streaming parameter could be calculated from these 
data. It will be shown later that the electron 
backstreaming parameter is somewhat dependent on the 
beam current or normalized perveance per hole. The 
value of this parameter at the impingement-limited 
perveance Is probably of the most interest, however, 
so most of the backstreaming data was collected at 
this high perveance condition. After the impinge- 
ment-limited perveance and backstreaming data had 
been collected, the grid separation distance was 
reduced and the preceeding procedures were repeated. 

Because of the potentially destructive nature 
of the electrical breakdown phenomenon, the elec- 
trical breakdown characteristic of each of the 
graphite grid pairs were measured after the per- 
veance/impingement characteristics and electron 
backstreaming data had been obtained. The back- 
ground pressure was held constant for electrical 
breakdown testing at a pressure of 2.0 x 10* 5 To rr 
because this represented an average value of the 
bell jar pressures encountered during all testing. 
The data was collected by establishing a given grid 
spacing with the ion source operating at a moderate 
beam current. The screen and accelerator voltages 
were then increased in magnitude until electrical 
breakdown occurred. These tests indicated that 
either the screen or accelerator voltage could be 
increased to induce breakdown, i.e. the sum of the 
two voltages was the critical variable in determin- 
ing when breakdown occurred. 

Results 


Figure 3 displays a set of typical curves 
showing how impingement current (normalized using 
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Fig. 3 Effect of Grid Separation on Impingement 
Curves - Nineteen Hole, Two-Grid Systems 
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beam current) was observed to vary with beam cur- 
rent and grid separation distance. The figure shows 
that closer grid separations facilitate operation at 
higher beam currents before direct ion impingement 
levels become excessive. The data of Fig. 3 also 
show that the sharpness of the "knee" or rate of 
increase to high impingement becomes less dramatic 
as the grid spacing is reduced. These curves per- 
tain specifically to the impingement data collected 
for the screen and accel grid hole diameters of 
1.5 mm at a net accelerating voltage of 500 V and a 
net-to-total accelerating voltage ratio (R) of 0.7. 
Data similar to those in Fig. 3 were collected at 
each of the screen hole diameters, net accelerating 
voltages and net-to-total accelerating voltages in- 
vestigated. Instead of presenting all of the per- 
formance data in the form shown in Fig. 3, it is 
more convenient and more meaningful for comparison 
purposes to define a maximum beam current for each 
grid separation, i.e., the beam current at which 
the direct accelerator grid impingement current be- 
comes excessive. Since the total accelerating 
voltage, the grid separation distance, and the hole 
size are known, this impingement-limited beam cur- 
rent can be divided by the number of apertures and 
then used to calculate the impingement-limited 
normalized perveance per hole using Eqs. 3 and 4. 

The normalized perveance based on this definition 
facilitates comparison of the ion extraction cap- 
abilities of various grid geometries over a range 
of operating conditions. The impingement condition 
selected somewhat arbitrarily to define the lmpinge- 
ment-limited beam current corresponds to an lmpmge- 
ment-to-beam current ratio of one percent. It 
should be noted here that a normalized impingement 
current of 1% is only slightly above the baseline 
(charge exchange induced) impingement level shown 
in Fig. 3. This limit on the beam current was 
based on such a modest change in impingement current 
because it yielded consistent results close to the 
operating condition where an ion thruster would be 
expected to operate. A careful examination of repre- 
sentative data showed that this criterion yielded 
impingement-limited beam currents that on the 
average were 10% t 3% higher than those beam cur- 
rents that might be associated with negligible 
direct ion impingement. 

Effect of Hoi e Size 

Using the impingement-limited normalized 
perveance per hole, obtained from plots like those 
of Fig. 3, as the performance yardstick, the results 
of Figs. 4, 5, and 6 were obtained. These figures 
compare the performance of the various screen hole 
sizes investigated as a function of the normalized 
grid separation distance. Each figure presents data 
for a different net accelerating voltage. For each 
voltage, the curves corresponding to the different 
hole sizes all exhibit similar behavior. As the 
normalized grid separation is reduced the perveance 
remains constant or increases slightly and then It 
falls off as the grid separation drops below the 
value where focusing becomes poor and accelerator 
impingement begins to limit performance. The poor 
focusing that causes this reduction in performance 
is believed to be a result of intra-grid equlpoten- 
tial lines which become progressively less planar 
as the grid separation is reduced. 

Figures 4, 5, and 6 do not appear to show any 
dramatic trends that would indicate a decrease in 
performance with decreasing hole size. This obser- 
vation differs from what was expected based on pre- 
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Fig. 4 Effect of Aperture Diameter on Two-Grid 
System Ion Extraction Performance for 
a Net Acceleration Voltage of 500 V 
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Fig. 5 Effect of Aperture Diameter on Two-Grid 
System Ion Extraction Performance for 
a Net Accelerating Voltage of 300 V 
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Fig. 6 Effect of Aperture Diameter on Two-Grid 
System Ion Extraction Performance for 
a Net Accelerating Voltage of 100 V 
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vious experimental investigations that reported a 
decrease in performance for screen hole diameters 
below 2.0 mm. 1 * 2 * 3 It is believed that results 
obtained in these tests differed from those measured 
in previous ones because of the experimental pro- 
cedures used. During the present study it was 
found for example that it was important that the 
ion source be capable of producing Bohm current 
densities toward the grids significantly greater 
than the current densities being demanded by the 
grids when they were operating near the impingement 
limit. This ensured that It was the grids rather 
than the ion source that was limiting ion extrac- 
tion. To accomplish this the source must be opera- 
ted at the progressively higher plasma densities 
and hence smaller Debye lengths that are demanded 
as hole size and grid separation are reduced. It 
is noted in addition that some adjustment in flow 
rates was required to achieve the discharge chamber 
operating condition that would yield the highest 
impingement-limited current densities. It was also 
determined that the stability of the discharge 
chamber influenced the ion extraction capabilities 
of the grids. During a segment of the tests the 
discharge chamber was changed to a design that 
operated less stably (arc discharge noise being 
being apparent on an oscilloscope). The onset of 
high impingement currents occurred at lower beam 
current levels with this design than it did with 
the more stable one. 


To help visualize the dramatic effect that 
decreasing hole size has on current density the 
data of Fig. 4 has been replotted in Fig. 7 with 
current density rather than perveance as the de- 
pendent variable. The current densities shown here 
were calculated by dividing the impingement-limited 
beam current per hole by the area of a screen hole. 
It is noteworthy that the current densities con- 
tinue to rise as grid separation is reduced below 
the value where perveance began to drop off in 
Fig. 4 (z g /d s = 0.6 to 0.7) and that current den- 
sities Increase markedly as hole size is reduced. 
The current density of n, 25 mA/cm 2 measured with 
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Fig. 7 Effect of Aperture Diameter on Extracted 
Ion Current Density 


the 1.0 mm diameter holes is observed to be the 
maximum value realized at a net acceleration volt- 
age of 500 V for the case considered here where the 
screen and accel grid hole diameters were equal and 
the impingement was limited to 1% of the beam cur- 
rent. When impingement levels of several percent 
were allowed this maximum current density increased 
to about 35 mA/cm 2 . This substantial increase in 
current density can be realized when Impingement is 
allowed to increase because impingement does not 
rise percipitously with beam current when the grids 
are spaced very close together (see Fig. 3 for 
example) . 

Effect of Net Accelerating Voltage 

Although there is no apparent reduction in the 
impingement-limited perveance capabilities of grids 
as hole size is reduced. Figs. 4, 5, and 6 do show 
that there is a definite trend toward lower lmpinge- 
ment-limlted perveance per hole values as net accel- 
erating voltages are reduced. Aston 3 has also 
reported that the net acceleration voltage level can 
affect the impingement-limited perveance capabili- 
ties of grid systems. He suggests that the ratio 
of discharge voltage-to-total accelerating voltage 
is a suitable correlating parameter for data of this 
sort. This parameter has been used in Fig. 8 to 
show how changes in total accelerating voltage as- 
sociated with the data in Figs. 4, 5, and 6 induce 
changes in the impingement-limited perveance per 
hole. The data symbols and error bars in Fiq. 8 
pertain to the data from the horizontal portions of 
the curves in these figures. Figure 8 shows an 
almost linear decrease in performance as the ratio 
of discharge-to-total accelerating voltage 
increases. 
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DISCHARGE TO TOTAL ACCELERATION VOLTAGE RATIO 

Effect of Discharge-to-Total Acceleration 
Voltage Ratio on Ion Extraction Performance 


Discharge and total accelerating voltages are 
known to alter the screen hole sheath shapes and 
positions and this could be the cause of the effect 
shown in Fig. 8. One might also speculate that 
this reduction results from a plasma related block- 
age of the screen grid hole. Such a situation 
could occur at low net accelerating voltages where 
low current densities with their attendant low 
plasma densities and larger Debye lengths are ob- 
served. It is also considered possible that 
discharge voltage influences the off-axis velo- 
city component of ions entering the sheath while 
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total accelerating voltage determines their on-axis 
component. The ratio of these two voltages would 
then be expected to determine the relative trajec- 
tory of ions toward the accelerator grid webbing. 

In any case, this effect is believed to be a real 
effect and to avoid substantial decreases in per- 
formance at very low total accelerating voltages, 
it appears that it will be necessary to use thrus- 
ters which exhibit stable operation at high current 
densities with low discharge voltages. 


Effect of Net-to-Total Accelerating Voltage Ratio 


Figure 9 shows the effect of variations in the 
ratio of the net-to-total accelerating voltage on 
the impingement-limited perveance per hole charac- 
teristics of grids as a function of normalized grid 
separation for a typical case (d s = 1.5 mm, 

Vfj = 300 V). It shows that increases in the net-to- 
total acceleration voltage ratio (R) cause the 
characteristic impingement-limited perveance vs. 
grid separation curves to shift to the left. This 
means that operation is possible at a given net 
acceleration voltage at closer grid separation and 
hence higher current densities at the higher net-to- 
total acceleration ratios (R). This appears to be 
a result of the improved focusing associated with 
the higher values of R. Although Aston^ reported 
large variations in beam divergence with R, his 
data did not indicate such a dramatic effect of R on 
the impingement-limited perveance per hole as that 
suggested by the results of Fig. 9. The reason for 
this difference is uncertain at this time but it may 
be related to a difference in the screen-to-accel 
grid hole diameter ratio which was different in the 
two tests. 
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Fig. 9 Effect of Net-to-Total Acceleration Voltage 
Ratio on Ion Extraction Performance for 
1.5 mm Diameter Grid Apertures 

The significance of the results shown in Fig. 9 
can be seen more readily if the impingement-limited 
current density (current per hole/hole area) is used 
as the dependent variable in place of the lmpinge- 
ment-limited normalized perveance per hole. This 
has been done in Fig. 10. Even though the poten- 
tial difference between the grids decreases as the 
net-to-total acceleration voltage ratio is increased, 
Fig. 10 indicates the current density can be in- 
creased by doing this. This can be accomplished 
because increasing the net-to-total acceleration 
voltage ratio facilitates better focusing and the 


opportunity to reduce the grid separation to pro- 
gressively lower values before accelerator grid 
impingement becomes excessive. Unfortunately, the 
process of increasing the net-to-total accelerating 
voltage cannot continue indefinitely, because the 
operation of the grids is eventually limited by the 
onset of backstreaming. It should be noted that the 
results of Fig. 10 were obtained with an accelera- 
tor grid hole diameter equal to the screen hole 
diameter and the observations just made could be 
altered somewhat if the accelerator grid hole dia- 
meter were reduced. 
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Fig. 10 Effect of Net-to-Total Acceleration Voltage 
Ratio on Current Density Capability of 
1.5 mm Diameter Grid Apertures 


Electron Backstreaming 

Figure 11 shows the effect of variation in the 
screen/accelerator hole size on the electron back- 
streaming parameter as a function of normalized 
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Fig. 11 Effect of Screen/Accelerator Grid Aperture 
Diameter on the Electron Backstreaming 
Parameter 


115 



grid separation length. Qualitatively, all the 
curves are observed to be similar. As the norma- 
lized grid separation Is reduced, the electron back- 
streaming parameter remains fairly constant, and 
then falls off slightly at the lower grid spacings. 
The reason for this drop off is uncertain at this 
time, but it should be noted that at the lower grid 
separations under some conditions it was difficult 
to distinguish between Increases In beam current 
due to Improved focusing and increases due to back- 
streaming. There Is also a trend In Fig. 11 toward 
lower electron backstreaming parameters with de- 
creasing hole size. It Is also possible that this 
lowering of the backstreaming parameter might be 
related to an Increase of screen webbing area be- 
tween the holes which accompanied reductions in 
screen hole diameter In the present experiments. 

This occurred because the hole-to-hole spacing was 
held fixed while the hole diameter was varied. 
Additional experiments are needed to determine which 
of these parameters induce the observed effect. The 
less significant effect of net accelerating voltage 
on the electron backstreaming parameter, for the 
typical case where the screen hole diameter is 
1.5 mri. Is shown in Fig. 12. 
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Fig. 12 Effect of Net Acceleration Voltage on the 
Electron Backstreaming Parameter for 1.5 urn 
Diameter Grid Apertures 

Figure 13 shows the typical results indicative 
of the effect of the perveance at which a grid Is 
operating on the electron backstreaming parameter. 
The rate at which the electron backstreaming para- 
meter decreases appears to increase with decreasing 
perveance levels. Similar qualitative results were 
reported by Kaufman on the basis of numerical so- 
lutions for the electron backstreaming problem. 0 
In summary, Figs. 11 to 13 Indicate that there are 
secondary effects due to hole size (or open area 
fraction), net accelerating voltage, and perveance 
on the electron backstreaming parameter. Disregard- 
ing these minor effects the figures suggest that 
there is an upper limit on the electron backstream- 
ing parameter of n, 0.10. 

Electrical Breakdown 

The apparatus used in the conduct of the ex- 
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Fig. 13 Effect of Perveance on the Electron Back- 
streaming Parameter for 2.0 mm Diameter 
Grid Apertures 

periments described herein were particularly well 
suited to the conduct of tests on the electrical 
breakdown characteristics of grid pairs for ion 
sources. The reasons for this were that 1) the 
grid separation could be varied to determine when 
breakdown would occur and 2) the grids were small 
so the effects of electrostatically and thermally 
induced deformations could be minimized. Figure 14 
shows that the electrical breakdown characteristics 
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14 Electrical Breakdown Characteristics of 
Graphite Grid Pairs 
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of each of the different graphite grid pairs re- 
mained fairly constant over the range of grid sepa- 
rations investigated in the range of 3 to 5.5 kV/mm. 
All of the data lies above an electrical breakdown 
limit of 3 kV/mm. This is slightly higher than 
generally accepted limit of 2 kV/mm for ion thruster 
accelerator systems. ^ No special effort was made to 
ensure the quality of the finishes on the grid sur- 
faces and it is believed that this caused the scat- 
ter in the electrical breakdown limits shown in 
Fig. 14. There is no effect of grid aperture dia- 
meter or grid spacing on the breakdown electric 
field apparent from the data of Fig. 14. Additional 
tests also suggested that neither the net accelera- 
tion voltage nor the perveance at which the grids 
were operating influenced this breakdown electric 
field measurably. 


Conclusions 


Apparatus for testing two-grid accelerator 
systems has been developed which is well suited to 
evaluating the impingement-limited current density 
capabilities and the electron backstreaming and 
electrical breakdown characteristics of such 
systems. The basic relationships defining the cur- 
rent density capabilities of ion optic systems 
appear to be valid for screen grid hole diameters as 
small as 1 mm and it is expected that these relation- 
ships will be valid for even smaller screen grid 
hole diameters. Tests show, however, that the 
impingement-limited perveance per hole at which the 
grids can be operated degrades as the ratio of dis- 
charge-to-net accelerating voltage increases. The 
results also show that operation at the highest 
current densities are realized at close grid sepa- 
rations when the net-to-total accelerating voltage 
ratio is high. In general, at close grid separa- 
tions there is an optimum range of net-to-total 
accelerating voltages. To operate in this optimum 
range, the ratio of net-to-total accelerating volt- 
age must be high enough to avoid poor focusing yet 
low enough to prevent electron backstreaming. Elec- 
tric fields greater than 3 kV/mm were obtained be- 
fore electrical breakdown occurred between the pairs 
of graphite grid tested. 
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Abstract 

An experimental investigation Is described in 
which a comparison of the ion beam neutralization 
characteristics of a local neutralizer (within 
5 cm of the beam edge) with those associated with 
a distant one (n- 1 meter away from thruster) were 
compared. The Influence of magnetic fields in the 
vicinity of the neutralizer cathode orifice which 
were either parallel or normal to the neutralizer 
axis were studied. Plasma property profiles which 
reflected the Influence of the magnetic fields were 
measured. The results indicate that magnetic fields 
at the region of a neutralizer cathode orifice in- 
fluence its ability to couple to the ion beam. 

They show that there is a potential jump from the 
neutralizer cathode orifice to the plasma which 
exists close to the orifice. This potential drop 
increases as the axial component of magnetic flux 
density increases. A magnetic field perpendicular 
to the neutralizer axis induces a potential rise a 
few centimeters further downstream of the neutral- 
izer cathode. Test results Indicate the path 
length from the neutralizer to the beam along mag- 
netic field lines also has an effect on the beam-to- 
neutralizer plasma potential difference but this 
effect was typically smaller than the other two 
effects. The results suggest that in order to min- 
imize the adverse effects of stray magnetic fields, 
a region as free of magnetic fields as possible 
should be selected for placement of the neutralizer. 


Introduction 

Recent tests conducted in space on the SERT II 
Spacecraft included an Investigation of the process 
of ion beam neutralization from an electron source 
located far from the Ion beam. This was accom- 
plished by using the neutralizer or unaccelerated 
thruster plasma of a thruster located a meter or so 
away from the thruster operating at high voltage.’ ^ 
One of the results observed In this test was that 
the potential difference between the neutralizer 
and beam plasma potential could be lower with the 
distant neutralizer than It was when the local one 
was used. In order to determine If local magnetic 
fields could be the cause of this and to develop 
an understanding of neutralization phenomena that 
would aid in the design of future neutralizers, 
the study described herein was conducted. The ex- 
periments undertaken included a comparison of the 
ion beam neutralization characteristics of a local 
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neutralizer with those associated with a distant 
one. A study of the Influence of magnetic fields 
In the vicinity of a neutralizer cathode orifice 
on ion beam neutralization was also undertaken. 
The plasma property profiles which reflected the 
effects of changes In magnetic fields near the 
neutralizers were measured. 


Apparatus and Procedure 

The experiments were conducted In the 1 .2 m 
diameter by 4.6 m long stainless steel tank shown 
schematically In Fig. 1. The operating pressure 



Fig. 1 Neutralization Test Apparatus 


in the vacuum system was in the range of 5 v 
7 x 10" 6 torr during the tests. At this pressure 
the electron collision frequency Is small compared 
to a typical cyclotron frequency so the conditions 
of this test should model the colllslonless en- 
vironment of the space tests adequately. Further 
the potential of the plasma In the vacuum tank was 
typically several volts above ground, so most neu- 
tralizer electrons should have been reflected from 
the tank walls; a condition that also facilitates 
modeling of the space plasma. The SIT-8 mercury 
ion thruster 3 identified In Fig. 1 Is located at 
one end of the test facility In such a way that its 
Ion beam is directed along the tank axis. Neutral- 
ization of this beam is accomplished by either the 
standard SIT-8 hollow cathode mercury neutralizer 3 
located adjacent to the thruster or by the distant 
neutralizer shown in Fig. 1. Both neutralizers 
utilize conventional 0.32 cm dla. mercury hollow 
cathodes. Each of the neutralizers Is equipped 
with a rolled tantalum foil insert, and a toroidal 
keeper electrode located a. 0.08 cm downstream of 
the orifice plate. The distant neutralizer cathode 
has a 0.04 cm dla. orifice while the local one Is 
0.025 cm in diameter. In order to achieve stable 
and reproducible operation of the neutralizers, the 
distant one was operated at a mercury flow rate of 
^ 60 mA while the local one operated at ^ 12 mA. 

The lower flow rate was presumably required for the 
local neutralizer because It had a smaller orifice 
diameter and because it utilized an enclosed keeper 
configuration. 

The distant neutralizer was attached to a rod 
passing through the top of the tank as shown In 
Fig. 1. When this rod was rotated the neutralizer 
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moved through a 180° arc from one side of the tank 
to the other. Because the vacuum tank Is stainless 
steel the influence of the earth's magnetic field 
Is felt within It. This field has a magnitude of 
■v 0.5 gauss and is Inclined at an angle of 65° from 
the direction of the Ion beam In the plane of 
section A-A In Fig. 1. Movement of the neutralizer 
support rod through its 180° travel Is sufficient 
to place the neutralizer at positions where It lies 
on earth magnetic field lines that either Intersect 
the ion beam centerline or miss It by distances up 
to -v 30 cm. Each neutralizer cathode was equipped 
with power supplies that controlled keeper poten- 
tial, cathode tip heater current and neutralizer 
bias potential. Keeper currents for the cathodes 
were maintained at either 0.3 or 0.4 A and the neu- 
tralizer bias was adjusted so It was sufficiently 
negative to effect the neutralizer emission re- 
quired to match the ion beam current. Auxiliary 
magnetic fields of variable magnitude were gen- 
erated In these experiments by using the two elec- 
tromagnets shown in Fig. 1. These magnets were 
located Immediately adjacent to the particular neu- 
tralizer they served. Magnetic fields of the order 
of 15 to 30 gauss could be Induced at the neutra- 
lizer cathodes with these magnets. 

Evaluation of the effects of parametric varia- 
tions introduced in the experiments required the 
measurement of the neutralizer bias voltages, pro- 
files of the plasma potential, and magnetic field 
shapes and magnitudes. The neutralizer bias volt- 
age could be measured directly but the beam plasma 
potential had to be measured by an emissive probe 
as suggested by Fig. 1 . 

In a second series of experiments the distant 
neutralizer was moved to the vicinity of the vacuum 
feedthrough in the configuration shown in Fig. 2. 

In this configuration the neutralizer cathode axis 
was pointed at the ion beam rather than away from 
it as it had been in the configuration shown in 
Fig. 1. It was also installed so emitting and non- 
emitting spherical Langmuir probes could be used to 
measure plasma properties along the neutralizer 
axis through that region ranging from Immediately 
adjacent to the distant neutralizer into the ion 
beam. The emissive probes were used to measure 
plasma potential profiles while the non-emitting 
one gave electron density and temperature informa- 
tion as well. The details of the operation of 
these probes are described in Refs. 4 and 5. In 
this second series of tests two electromagnets were 
installed to facilitate control of the magnetic 
field in the vicinity of the distant neutralizer In 
the manner suggested In Fig. 2. The intensities 
and shapes of magnetic fields in all of these tests 
were determined from gaussmeter and Iron filings 
map measurements. 


Experimental Results 

Neutralization in the Unaltered Environment of 
the Magnetic Fields of the Thruster and Earth 

The initial Investigations into the neutrali- 
zation of the ion beam were conducted with the 
electromagnets shown In Fig. 1, turned off. The 
magnetic field in the region between the distant 
neutralizer and ion beam Is In this case the weak 
one (t. 0.5 gauss) associated with earth only. In 
the region between the local neutralizer and ion 
beam the net magnetic field as induced by the 
thruster and earth was about 18 gauss. At a beam 




Fig. 2 Neutralization Test Apparatus 

current of 100 mA with the distant neutralizer 
positioned as shown in Fig. 1, so electrons emitted 
from It would pass about 30 cm from the beam 
centerline, the plasma potential profiles repre- 
sented by the dotted and dashed curves of Figure 3 
were observed. These profiles were measured 
through the Ion beam at the Indicated distances 
downstream of the Ion source screen grid and during 
data collection for these curves, the local neu- 
tralizer was not operating. When the distant neu- 
tralizer was turned off and the local neutralizer 
was started the higher plasma potentials represent- 
ed by the sol id and centerl Ines were observed at the 
two axial locations Identified. The higher plasma po- 
tentials observed here with the local neutralizer 
are reminiscent of results obtained during some of 
the SERT-II tests 1 * 2 conducted in space. The 
plasma potential profile at the 6 cm location Is 
also seen to be broader when the local neutralizer 
is being used, a condition that was also observed 
in the space tests. Both of these effects are be- 
lieved to be manifestations of an impedence between 
the beam and local neutralizer plasmas that Is 
higher than the one between the distant neutralizer 
and beam plasma. This postulate is also supported 
by the observation that the neutralizer cathode 
bias voltage relative to ground (V c ) required to 
effect neutralization was also considerably lower 
in magnitude when the distant neutralizer was used. 

In a second experiment the local neutralizer 
was secured, the distant neutralizer was moved with 
respect to the ion beam while the thruster was 
operating. The effect of this movement on the Ion 
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Fig. 3 Effect of Neutralizer Location on Ion 
Beam Plasma Potentials 

beam plasma profile was examined. This neutralizer 
movement was accomplished by rotating the arm sup- 
porting the distant neutralizer (Fig. 1) so the 
neutralizer plasma was moved from the earth mag- 
netic field line missing the ion beam centerline 
by n. 30 cm to one passing through the ion beam 
centerline and the beam potential profiles shown 
in Fig. 4 were obtained. As indicated on the fig- 
ure, profiles were measured at axial locations 
downstream from the grids (x) of 6, 29 and 52 cm. 
The results shown in this figure are as expected 
in that the beam plasma potentials dropped when the 
neutralizer was moved to the field line intercept- 
ing the beam axis. The fact that this drop was 
small ('v 3v) could be due to the fact that the neu- 
tralizer plasma was not moved very far away from 
the edge of the beam plasma. In this regard it is 
noted that 1) the beam plasma potentials suggest 
substantial beam spreading with axial position and 
2) the electron gyro radius in the tank plasma and 
earth's magnetic field is about 10 cm. The fact 
that the gyro radius is less than the 30 cm dis- 
tance between the magnetic field line passing 
through the neutralizer plasma and the beam center- 
line is consistent with the observed change in the 
plasma potential profiles. It is also noted that 
the neutralizer bias voltage (Vg) required to 
effect neutralization of the 100 mA beam current 
tracked with the changes in plasma potentials. It 
has been suggested 6 that fluctuations in beam 
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Fig. 4 Effect of Distant Neutralizer Position in 
the Earth's Magnetic Field 

plasma potentials of the order of several volts 
observed during the SERT-II space tests could have 
been caused by spacecraft motion in the earth's 
magnetic field. These fluctuations were originally 
attributed to a spacecraft wake effect but the 
relative motion of the distant neutralizer and beam 
in the earth's magnetic field could have also 
caused this effect in much the same way it was in- 
duced in the test just described. 


The Influence of Magnetic Field at the Distant 
Neutral iz~er 

For this series of tests the electromagnet at 
the distant neutralizer in the Fig. 1 test config- 
uration was positioned so its solenoidal magnet 
field was aligned with the neutralizer axis as sug- 
gested by the inset sketch in Fig. 5. As the elec- 
tromagnet current was varied the shapes of the 
field lines shown in Fig. 5 were not altered but 
the field strengths did change. The beam plasma 
potential profiles at axial locations of 6, 29 and 
52 cm downstream of the grids were measured as this 
change was effected. The relative shapes and mag- 
nitudes of these profiles at these three axial 
positions were however consistent with the trends 
shown in Figure 5. As a result the plasma poten- 
tial at the beam centerline 6 cm downstream of the 
grids was selected as representative of the state 
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Fig. 5 Effect of Magnetic Field at Distant 

Neutralizer on Beam Plasma and Neutralizer 
Potential s 

of the beam plasma potential. Figure 5 shows that 
the magnitudes of both this representative beam 
plasma potential and the neutralizer bias voltage 
increase as the magnitude of the solenoidal mag- 
netic field at the neutralizer tip is increased. 

The same behavior is observed at beam currents of 
100 mA and 200 mA although the magnitudes of both 
voltages are higher at the higher beam current. 

Figure 6 shows that the coupling voltages be- 
tween the neutralizer and beam plasmas (l.e., total 
voltage drop from beam plasma to neutralizer tip) 
taken from the data of Figure 5 are nearly directly 
proportional to magnetic flux density at both beam 
currents. Further, these voltage difference levels 
increase In direct proportion to the beam current 
level. Visual observation of the neutralizer 
plasma when the solenoidal magnetic field current 
was Increased showed that an increase in the field 
strength caused the luminous region associated with 
the cathode discharge to contract. As the magnetic 
field strength was Increased to n. 15 gauss, the 
hollow cathode discharge appeared as a very tiny 
bright spot at the orifice. 

In order to further clarify the influence of 
the magnetic field on the neutralizer's ability to 
couple to the ion beam, another series of tests 
were undertaken. In these tests the distant neu- 
tralizer was positioned so its axis was perpendic- 
ular to the thruster axis as shown in Fig. 2. The 
magnetic fields were arranged so one was aligned 
with the neutralizer axis and another was normal to 
it as suggested by the sketch of section B-B in 
Fig. 2. Each of the magnets could be controlled 
independently. Again as the current was varied to 
either electromagnet the shapes of the field lines 
were not altered but the field strengths did change, 

The plasma potential profiles along the neu- 
tralizer axis from the cathode to the Ion beam 
center observed when the axial magnetic field was 
on and the normal magnetic field was off are shown 
in Fig. 7a. These curves show that there Is a 
steep potential jump between the neutralizer 
cathode and the plasma close to the cathode orifice. 
Most of the voltage increase between the neutra- 
lizer cathode and the beam plasma Is seen to occur 
in this region. Further the magnitude of this 



Fig. 6 Effect of Magnetic Field at the Distant 
Neutralizer on Voltage Difference Between 
the Neutralizer and Ion Beam 


voltage rise is observed to increase with axial 
magnetic flux density at the neutralizer tip. It 
is noteworthy that the data points at the cathode 
orifice (on the vertical axis) In Fig. 7a corres- 
pond to the neutralizer bias voltage while those 
further downstream were measured using the emis- 
sive probes. After the initial voltage jump, its 
rise is observed to be approximately linear with 
distance and of the order of a few to several volts. 
As the axial magnetic field (the component along 
the direction of the neutralizer axis) was in- 
creased by increasing the current to the axial 
electromagnet (Fig. 2), Fig. 7a shows the magni- 
tudes of the plasma potentials Increased and the 
cathode potential required to Induce a 100 mA emis- 
sion current decreased. Fig. 8 shows this same 
effect In the form of plots of the beam plasma 
(i.e., 60 cm from the neutralizer), the neutralizer 
plasma (i.e., ■v 0.3 cm from the neutralizer) and 
the neutralizer cathode potentials as a function of 
the magnetic flux density at the neutralizer ca- 
thode orifice. It is noted the plasma density pro- 
file measured along the neutralizer axis and shown 
in Fig. 9 as function of axial magnetic field does 
not seem to be affected significantly by this field. 
These results suggest that the axial magnetic field 
affects the capacity of a neutralizer to couple to 
the ion beam principally because it induces a po- 
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a) PLASMA POTENTIAL PROFILES ON NEUTRALIZER AXIS 



b) AXIAL MAGNETIC FLUX DENSITY PROFILES 
ON NEUTRALIZER AXIS 


Fig. 7 Effect of Axial Magnetic Field 


100mA BEAM CURRENT 




Fig. 9 Electron Density Profile for Axial Magnetic 
Field Case 

tential drop close to the cathode orifice. Finally 
it is noted that the coupling voltage between the 
neutralizer and beam plasma is nearly directly pro- 
portional to the magnetic flux density. 

Plasma potentials were also measured for the 
case where the axial magnetic field was turned off 
and the normal magnetic field was turned on. 

Figure 10 shows how this magnetic field component 
affects the plasma potential profile. In contrast 
to Fig. 8, where the axial field was observed to 
cause an increase in the potential difference very 
close to the cathode (< 0.3 cm), the normal magnetic 
field induces a potential difference over a much 
longer distance (a few centimeters). The Increase 
in potential difference through this region is pre- 
sumably caused by the plasma impedance induced by 
the magnetic field aligned perpendicular to the di- 
rection of current flow. Figure 11 shows that the 
variation in plasma density along the neutralizer 
axis when the normal magnetic field is applied. 
Comparison of this Figure with Fig. 9 reveals a 
much more dramatic drop in plasma density along 
the neutralizer axis when the normal field is 
applied in place of the axial one. In Figure 12 
the potentials of the beam plasma, neutralizer 
plasma and neutralizer cathode are plotted as a 
function of the normal magnetic field flux density. 
This Figure again shows that the potential rise at 
the cathode does not increase in the immediate 
vicinity of the cathode as it did in the case of 
the axial magnetic field. As suggested by Fig. 10 
most of the potential rise occurring between the 
beam and neutralizer plasmas is observed within a 
few centimeters of the neutralizer cathode. 


The Influence of the Magnetic Field at the Local 
Neutralizer 


Fig. 8 Effect of Axial Magnetic Field on The test results obtained with the local SIT-8 
Neutralizer and Beam Potentials neutralizer showing the influence of stray magnetic 

fields in its vicinity on neutralization perfor- 
mance indicators are shown in Fig. 13. In this 
case the neutralizer is located in the magnetic 
field associated with the thruster Itself; this 
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b) NORMAL MAGNETIC FIELD FUJX DENSITY 
PROFILES ON NEUT AXIS 


Fig. 10 Effect of Normal Magnetic Field 

field is oblique to the local neutralizer axis. The 
test involved altering this field with an external 
magnet and measuring the resulting changes in beam 
and neutralizer bias potentials. The central iron 
filings map in Figure 13 shows the unaltered mag- 
netic field associated with the thruster as well as 
the location of the neutralizer in this field. In 
the filings maps to the left of center the external 
magnet acts to cancel the thruster field thereby 
reducing the magnitude of magnetic flux density at 
the neutralizer and bending the field lines passing 
through the neutralizer toward the ion beam. In 
the filings maps on the right, the external magnet 
is used to increase the magnetic field magnitude 
and draw the field lines passing through the neu- 
tralizer still further away from the beam. The 
arrows associated with each filings map indicate 
the magnitude of magnetic flux density at the neu- 
tralizer corresponding to that filings map. 

Figure 13 shows that when the magnetic flux density 
magnitude is Increased and the field lines passing 
through the neutralizer are drawn further away from 
the beam, the neutralizer bias and beam plasma po- 
tentials are Increased. Conversely, reducing the 
magnetic flux densities and directing the field 
lines toward the beam reduces the magnitudes of 
these potentials somewhat. As with the distant 
neutralizer (Fig. 5), operation at high beam cur- 
rents requires higher potential differences between 
the beam and neutralizer In order to effect neu- 
tralization. The results of Fig. 13 are observed 
to be qualitatively consistent with results obtain- 
ed with the electromagnets used in tests associated 
with the distant neutralizer. 



Fig. 11 Electron Density Profile along Neutralizer 
Axis - Normal Magnetic Field 


MAG FIELD 100mA BEAM CURRENT 

DIRECTION 



Fig. 12 Effect of Normal Magnetic Field on 
Neutralizer and Beam Potentials 
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Fig. 13 Effect of Magnetic Field at Local 

Neutralizer on Beam and Neutralizer 
Potentials 


A Phenomenological 
Discussion of Test Results 


The Influence of Axial Magnetic Field on 
Neutralization 


The experiments have suggested that the plasma 
between the neutralizer cathode and the ion beam 
can be described as three separable regions. The 
first will be defined as the thin layer (sheath) 
immediately adjacent to the cathode orifice. The 
next (intermediate) region is assumed to extend 
from this cathode region through the first few cen- 
timeters immediately downstream of it. The final 
(distant) region is assumed to extend from the edge 
of the second region to the ion beam. Based on the 
experimental results obtained the following model 
describing the influence of low intensity stray 
magnetic fields in these three regions is proposed: 

1) The sheath region is affected by the axial 
magnetic field, 

2) The intermediate region is affected by the 
normal magnetic field, and 

3) The distant region is not affected signif- 
icantly by either component of magnetic 
field. 

At the hollow cathode orifice emitted electrons 
and ions are both present so a double sheath is 
formed. When no magnetic field was applied in the 
experiments the bulk of the neutral! zer-to-beam po- 
tential rise occured across this sheath. As the 
axial magnetic field strength was increased this 
sheath potential rise increased in magnitude. It 
is postulated that with no magnetic field present 
the sheath is curved because the electrons exiting 
the orifice can move both radially and axially 
toward the keeper and/or the beam. As the axial 
magnetic field strength is increased it is argued 
that the electrons are confined to move along field 
lines and the sheath moves from a curved to a flat 
configuration. In order to model the effect of the 
axial magnetic field on the sheath potential rise 


consider the sheath with no field present to be 
hemi-spherically convex. 


For the purposes of this analysis the poten- 
tial through the sheath will be measured relative 
to the downstream sheath plasma boundary where the 
potential is taken to be zero. Position will be 
measured from the neutralizer cathode orifice exit 
surface where r = r 0 as shown in Fig. 14. Assuming 
the thickness of sheath to be one Debye length (Xd) 
analysis shows in the test under consideration that 
Xd < X m is satisfied where x m is the electron mean 
free path for elastic collisions. Therefore the 
particles are assumed to move without collision or 
ion production within the sheath. The sheath 
plasma is assumed to be composed of three particle 
groups, namely electron drawn from the cathode in- 
terior, ions from the downstream plasma, and elec- 
trons from this same plasma. The electrons from 
the cathode interior are assumed to enter the 
sheath with a negligible velocity and to have a cur- 
rent density j e and a mass m e . Ions entering from 
the other side have a current density ji, a mass mi 
and a velocity determined by the Bohm criterion. 

The potential equivalent of their velocity is given 


V 


o 


e 

2e 


(1) 


where T e is the electron temperature in the down- 
stream plasma. 



The plasma electrons entering the sheath from 
the downstream plasma will have a Boltzmann dis- 
tribution so that their density n e will decay ex- 
ponentially into the sheath. Their charge density 
can be expressed as 

P p = en = en exp [eV/kT ] (2) 

e l e e 0 e 
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where n eQ is the density of electrons at sheath 
edge where V = 0. These three plasma constituents 
combine to determine the potential variation 
through the plasma sheath. In this case we assume 
spherical symmetry and a sheath dimension that is 
small compared to the orifice radius so the effect 
of curvature can be neglected.' For this case 
Poisson's equation can be used to describe the 
physical situation in one dimension. 


d 2 V 

dr' 1 


\W \f 

"o\ m e o\ 


2 e(V 0 -V) 


+ ]: (en e )exp(eV/kT e ) 


(3) 


where c 0 is the permittivity of free space, 
the definition 
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Using 


and applying the space charge limiting boundary 
conditions 

= 0 
r = r o 


dV 

3F 


d V 
df 



where x s is the sheath thickness Eq. (3) integrates 
to 


d V = 
dr 
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The sheath thickness can be expressed as 
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where . . m h , 

T ( 2§) (-V c )^ 2n eo kT e 

a = 


— {— )* [(V -V )** - V **] 
A v 2e' LV o V o J 


(5) 

( 6 ) 


In this expression J is the sum of the keeper cur- 
rent and the ion beam neutralization current and A 
is the surface area of sheath through which the 
electron current flows. This area is taken to be 
the surface area of the sector of a sphere of radius 
R intersecting the cathode orifice circle. If one 


assumes the potential of electrons entering the 
sheath from the cathode Interior is approximately 
equal to the potential of the cathode and that the 
potential at the downstream sheath surface is equal 
to the potential which was measured ^ 3 mm down- 
stream of the cathode the sheath thickness in the 
absence of a magnetic field is found to be 
^ 0.024 mm from Eq. (5). 

In the presence of a magnetic field, as men- 
tioned earlier the surface should tend to become 
more planar (i.e., R increases). In the limit of 
a high magnetic field R -► ~ and the sheath surface 
is equal to the area of the cathode orifice. If we 
assume the sheath thickness is independent of the 
axial magnetic field strength then the difference 
in the potential rise from the cathode to the neu- 
tralizer plasma for these two extremes (R = ro and 
R = ») can be computed. This change In potential 
for the conditions pertaining to the present prob- 
lem is found to be -v 12v, (assuming T e = 0.8 eV 8 , 
and using measured values for the other quantities 
in Eqs. (5) and (6)). 

The measured changes in the neutralizer plasma- 
to-cathode potential difference ranged in value up 
to 7 or 8 volts at axial magnetic flux densities 
near 25 gauss. It is believed that magnetic flux 
densities near 50 gauss would be sufficient to in- 
duce the computed change in potential difference 
of 12 volts. At any rate the relative agreement 
between the measured and computed values is con- 
sidered significant. 

For this case of an axial magnetic field the 
experiments suggest that the potential drop through 
the intermediate and distant regions is only a few 
volts and that it varies linearly with distance. 
Axial magnetic field strength seemed to have a 
minimal effect on the voltage rise through these 
regions. It is postulated that this condition 
exists because the magnetic and electric fields are 
approximately parallel. Under this condition the 
conduction current dominates so the plasma conduc- 
tivity is not Influenced substantially by magnetic 
field strength. In addition the collision frequency 
is relatively small and this is also essential to 
the small rise observed through these regions. The 
small increase in potential rise that is observed 
through these regions with changes in axial mag- 
netic field is considered to be due primarily to 
small changes in the path length from the neutra- 
lizer to the ion beam. 


The Influence of a Normal Magnetic Field on 
Neutral i zation 


The experiments indicate that a normal mag- 
netic field effects a voltage rise in the inter- 
mediate region rather than the sheath region. It 
is argued in this case that electrons in the inter- 
mediate region are affected by both the electric 
and magnetic field that exist there and that these 
fields tend to be mutually perpendicular. It is 
also argued that in this intermediate region the 
neutral density is sufficiently low so the ratio 
of gyro frequency to collision frequency is 
moderate and as a result the Hall Current is sig- 
nificant. As the normal magnetic field is increased 
under these circumstances there should be a sub- 
stantial increase in the path length of electrons 
passing between the sheath and distant regions. It 
is believed that this increase in path length 
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through a region of moderate to low collision 
frequency causes the observed increase In potential 
difference. 

For this case where a normal magnetic field is 
applied the voltage rise through the sheath Is not 
affected by the field strength. This is expected 
because the sheath thickness is small compared to 
the gyro radius. In the distant region the magnet 
induces a small change in the magnetic field. The 
magnetic field there is quite weak and generally 
parallel with the electric field. Hence the effect 
of field strength on the potential rise through the 
distant region is relatively small as it was In the 
case where an axial magnetic field was applied. 


Combined Axial and Normal Magnetic Field Effects 

When the local neutralizer is used both normal 
and axial magnetic field components are present. 
Under this condition the axial magnetic field com- 
ponent should cause the potential rise at the 
sheath to increase while the normal magnetic field 
component should cause a potential rise through 
the intermediate region. Comparison of the results 
of Figs. 8 and 12 with those of Fig. 13 suggest 
that the potentials observed in this configuration 
are however higher than one would predict consider- 
ing both of the separate voltage rise components 
due to separate axial and normal magnetic fields. 

It is considered likely in this case that the con- 
duction path length In the distant Is Increasing 
substantially as magnetic field Is varied so that 
the potential rise in the distant region also be- 
comes significant. 

Conclusions 


and a neutralizer located on the order of 1 meter 

away. 
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The magnetic field at a neutralizer cathode 
orifice influences the neutralizer's ability to 
couple to the ion beam, A magnetic field which is 
aligned with the neutralizer axis influences the 
potential rise at the neutralizer cathode sheath 
which exists between the neutralizer cathode and 
the plasma immediately adjacent to it. This po- 
tential difference between the neutralizer and 
adjacent plasma increases linearly with axial mag- 
netic flux density. When a magnetic field that is 
perpendicular to the neutralizer axis Is applied 
the region over which the potential rise occurs is 
extended over several centimeters. 


Stray magnetic fields associated with Ion 
thruster themselves are sufficiently large to cause 
substantial increases in the potential difference 
between the neutralizer and the beam plasma. In 
order to minimize this effect a region as free of 
magnetic fields as possible should be selected for 
placement of the neutralizer. The neutralizer does 
not seem to have any difficulty coupling to the ion 
beam over distance of the order of 1 meter as long 
as the magnetic field levels are low (of order 0.5 
gauss or less). If a neutralizer must be located 
in the region of a stray thruster magnetic field, 
an external magnet which is properly arranged to 
alter the magnetic field so it is directed toward 
the ion beam can improve coupling. Application of 
such a field does not seem to alter thruster per- 
formance measurably. Magnetic fields of the order 
of those associated with the earth (^ 0.5 gauss) can 
effect small beam plasma potential changes (of order 
lv) when they are interposed between an ion beam 
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THE ELECTROTHERMAL RAMJET* 
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ABSTRACT Because of the high acceleration levels at 

which these devices could operate, the associated 
launch tracks could be relatively short and the 
An electrothermal ramjet conf irguration is terminal launch velocity could be achieved while 

examined as a possible alternative to rail guns and the P»yload was close to the earth's surface, 

mass drivers for high acceleration launch missions. Achieving terminal velocity in a short distance is 

For a specific mission (earth escape) the idealized considered desirable because the launch could be 

performance of the electrothermal ramjet, the elec- effected in a tube of reasonable size which would 

trothermal rocket and the electromagnetic accelera- be sealed during the launch sequence if it became 

tion system are compared. This comparison indi— necessary to abort the launch. This could facili— 

cates that the gross performance of the ramjet com- tate containment if toxic or radioactive materials 

pares favorably with that of the ideal electromag- were being launched, 

netic acceleration system. A specific configura- 
tion for the ramjet is chosen and models for the Another propulsion concept which may be useful 

dynamics, thermodynamics and fluid mechanics are 1° these high acceleration |evel missions is the 

presented. Results of calculations for a typical electrothermal rocket concept. In these devices 

supersonic launch cycle suggest that pressure, tear- on-board propellant is heated electrically and then 

perature and power demand profiles associated with discharged through a nozzle to produce thrust, 

ramjet operation should be reasonable. A light gas Because the heating is accomplished electrically, 

gun is proposed to accelerate the vehicle to the hydrogen propellant with its attendant high 

critical velocity where efficient ramjet operation specific impulse capability can be selected. There 

can begin. The theoretical performance of the ram- are several types of electrothermal rockets. In 

jet is also shown to be substantially better than the resistojet. for example, the propellant is 

that of the light gas gun at high velocities. heated by passing it over current-carrying, resis- 

tive elements and exhaust velocities as high as 
8S00 m/sec have been achieved. Higher exhaust 
velocities can be realized if the material limita- 
INTRODUCTION tions imposed by the heating elements are elim- 

inated and joule heating is accomplished by passing 
current directly through the propellant in an 
The relatively massive expenditure of chemical arcjet thruster. Exhaust velocities as high as 

fuel and oxidizer associated with the launch phase 15000 m/sec have been demonstrated in these dev- 

of space missions adds substantially to the cost of ices, 
such missions. Not only are the propellants 

costly, but the additional mass of tankage and A basic electrothermal rocket vehicle which 

structure required to hold the propellant mass might be used in this application is shown concep- 

detracts from the payload mass and therefore adds tually m Figure la. This rocket would carry on- 

substantially to the cost of the delivered payload. board propellant which would be metered into the 

It is also desirable to simplify these launch sys- heating chamber where its temperature could be 

terns and improve their reliability. For some mis- raised by energy either beamed there electromagnet- 

sions where high acceleration levels can be ically or carried there in the form of electrical 

tolerated such alternative launch schemes as light current from rails located in the launch tubo. 

gas guns 1 , rail guns , and mass drivers' 3 have been 

proposed as substitutes for conventional rockets. For the mission under consideration here. 

These missions could include ones in which various where the launch is accomplished in a relatively 

materials would be directed into apace either for short distance close to the earth's surface, it is 

disposal or application there. While the capital also practical to distribute the propellant in the 

investment in such systems is expected to be large, launch tube. When this is done a ramjet engine can 

they hold the promise of lower cost per launch be used to collect the propellant, heat it and pro- 

because of the large number of launches that should duce the desired thrust as it moves down the tube, 

be realized over their design lifetimes. Physically the ramjet would operate In the manner 

— — — suggested by Figure lb by simply drawing propellant 
• Work performed under NASA GRANT NGR-06-002-112 in fron the tnlle through an intake diffuser, heat- 

•* Professor, Dept, of Mech. Eng’g. ln * 11 electrically and then expelling it through a 

+ Research Asst., Dept, of Mech. Eng'g. nozzle. Energy could again be transported into the 

heating chamber either electromagnetically or in 
This paper Is declared a work of the US the forn of electrical current. No compressor or 

Government and therefore Is In the public domain other rotating equipment would be needed although 
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an initial velocity would have to be imparted to 
the vehicle. This might be accomplished by a 
booster rocket, a light gas gun or a low terminal 
velocity electromagnetic driver. 



Fig. 1. Electrothermal Thruster Concepts 


The purposes of this study are 1) to determine 
the potential of the electrothermal ramjet concept 
by comparing its gross performance capabilities 
with those of rocket and electromagnetic accelera- 
tor systems, 2) to select a ramjet configuration 
that appears to be promising and 3) to examine the 
fluid mechanics and thermodynamics of operation of 
this configuration. The execution of this work 
should point out any basic physical limitations 
that make the concept impractical and illustrate 
its potential advantages over other launch vehicle 
concepts. 


GROSS PERFORMANCE COMPARISON 


In order to compare the performance of various 
propulsion systems a high acceleration earth escape 
mission was selected. In this mission sufficient 
velocity would be imparted to the payload at the 
earth's surface so it would be at earth escape 
velocity when it reached an altitude where atmos- 
pheric drag would be negligible. The thrusters to 
be compared are characterized by their exhaust 
velocities with the following values being selected 
as typical of the capabilities of the devices indi- 
cated: 

Exhaust 

Thruster Velocity (D) 


Chemical Rocket 
Resistojet Rocket 
Arejet Rocket 
Resisto-ramjet 
Arc-ramjet 


4.500 m/sec 

8.500 m/sec 
15,000 m/sec 

8,500 m/sec 
15,000 m/sec 


Each of these thrusters will, for this prelim- 
inary analysis, be assumed to operate at a constant 
exhaust velocity. This is probably a reasonable 
assumption for the rockets, but it is not obvious 
that the ramjet is modeled adequately as a constant 
exhaust velocity device. Vork contained in the 


next section of this paper shows, however, that 
from a theoretical point of view it is quite rea- 
sonable to model the ramjet in this way. In addi- 
tion an ideal electromagnetic launcher will be used 
in the comparison. This device is idealized in the 
sense that all of the electrical energy input is 
assumed to be transferred to the moving payload 
mass and it therefore represents the ultimate 
launcher. 

The velocity required at the surface of the 
earth by a payload intended to escape the earth and 
the maximum acceleration that such a payload can 
endure are also nooessary as inputs to the problem 
under consideration. These quantities can be 
defined if one considers tho expulsion of the pay- 
load to consist of two phases, namely the launch 
phase illustrated in Figure 2 and the subsequent 
phase of vehicle passage through the atmosphere. A 
vertical launch trajectory is specified and ini- 
tially the vehicle is assumed to be at rest and to 
have a mass m^. When the thruster is started it 
generates a constant thrust F which produces an 
initial acceleration o . The launch phase is 
assumed to continue until time t^ when a terminal 
launch velocity v^ is reached. At this point the 
vehicle mass is m , the acceleration is o. and the 
axial location measured from the initial location 
is x.. Thrusting stops and the phase of vehicle 
passage through the atmosphere begins at this 
point. During this second phase the vehicle velo- 
city decreases as a result of atmospheric drag and 



Fig. 2. Launch Phase Nomenclature 


130 




The maximum acceleration occurs at the end 
of the launch phase. In order to minimize track 
length this acceleration should be as large as the 
vehicle can withstand, a value that is ultimately 
determined by the mechanical strength of the vehi- 
cle. If one assumes the basic vehicle is cylindri- 
cal and has a diameter d, the compressive stress S 
at the base of the vehicle is equal to the thrust 
force F divided by its cross sectional area. 


4F = 4m f g f g (1) 

nd 2 nd 2 



Fig. 3. Velocity Loss Through Atmospheric Drag 


In this equation o^, the final (maximum) 
acceleration, is in G's and g is the acceleration 
due to gravity. The ballistic coefficient of the 
vehicle B which characterizes vehicle drag during 
its movement through the atmosphere after launch, 
is given by 


B = 


c n « d 

4 m. 


(2) 


In this equation C is the vehicle drag coefficient 
and it should have a value nea£ 0.1 for the veloci- 
ties under consideration here. Combining Eqs. 1 
and 2 and solving for the maximum acceleration one 
obtains 



(3) 


Figure 3, which was obtained from data supplied by 
the NASA Lewis Research Center", shows the ratio of 
the velocity at the start of passage through the 
atmosphere to that after this passage as a function 
of the ballistic coefficient of the vehicle. For 
the case under consideration here the velocity 
after atmospheric passage is the earth escape velo- 
city as the vertical axis title suggests. From 
this figure it appears that a low but realistic 
ballistic coefficient might be 2x10 m /kg. At 
this value the ratio of atmospheric entrance 
velocity to the earth escape velocity (v “ 
11,200 m/sec) is just under 1.2 and the terminal 
launch velocity (v_) would therefore have to be 
about 13,000 m/sec. If one specifies an allowable 
stress in the vehicle of 1.4 x 10 nt/m ( 200,000 
psi) then the maximum allowable acceleration (Eq. 
3) becomes 30,000 G's. This acceleration, while 
high, is not considered unattainable. Light gas 
guns 2 have in fact been operated at projectile 
acceleration levels an order of magnitude higher 
than this. For our case, however 30,000 G's should 
not be exceeded if compressive stresses are to be 
held below 1.4 X 10 9 nt/m 2 during launch to a velo- 
city of 13 km/sec. This velocity will be suffi- 
cient to enable escape from the earth when the 
vehicle diameter is given by Eq. 2. For a 10 kg 
mass having the ballistic and drag coefficients 
mentioned this would correspond to a vehicle 
several centimeters in diameter. 


Results 


A simple mathematical model based on a dynami- 
cal force balance and an energy conservation 
analysis can be applied to each of the thruster 
cases considered. In each of these cases a con- 
stant thrust and constant exhaust velocity in the 
frame of reference of the rocket or ramjet have 
been assumed. The basic equations used are given 
in Reference 7. The results were obtained for a 
particular case where a 10 kg final mass was 
launched vertically on the earth escape trajectory. 
The 10 kg final mass was selected rather arbi- 
trarily but the basic trends and comparisons made 
here are not affected by this selection. 

CONSTANT THRUST 
CONSTANT EXHAUST VELOCITY 
<r, » 30000 G's 



Figure 4 shows the velocity profiles for the 
six thruster options being compared. The effect of 
increasing rocket exhaust velocity from 4,500 to 
15,000 m/sec is observed to be a three-to-fourfold 
reduction in the launch time and a two-to-threefold 
reduction in the required launch track length (x^) . 
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This points oat the great benefit of using the high 
exhaust velocity electrothermal rocket rather than 
the conventional chemical rocket motor which is 

characterized by the lower exhaust velocity. By 
using a ramjet rather than a rocket further reduc- 
tions in launch times and track lengths are real- 
ized at a given exhaust velocity. This occurs 
because there is no on-board propellant being 
accelerated in the case of the ramjet. At an 

exhaust velocity of 15,000 m/sec the ramjet per- 
forms essentially the same as the electromagnetic 
accelerator in that it requires the same track 

length (285 m) and launch time (44 msec). The 
8,500 m/sec ramjet is seen to exhibit performance 
comparable to the 15,000 m/sec rocket. 

Figure 5 shows the time profile of the vehicle 
acceleration for the various thrusters. The 

acceleration is limited to 30,000 G's and this 
implies considerably lower initial accelerations 
for the conventional rockets in which the vehicle 
mass must decrease with time as on-board propellant 
is discharged. The shorter launch times and track 
lengths given in Figure 5 occur of course for the 
cases where the acceleration remains closer to the 
limiting value throughout the launch. 

CONSTANT THRUST 

CONSTANT EXHAUST VELOCITY 

cr, ■ 30000 G's 



A discussion of the variation in acceleration 
for the 8.5 km/sec ramjet will facilitate an under- 
standing of some of the physical characteristics of 
this thruster. Before a ramjet can be operated it 
is necessary to accelerate it, using an alternative 
force, to that velocity where it is processing ade- 
quate propellant mass. For the ramjet cases in 
Fig. 5 this is arbitrarily presumed to be accom- 
plished using an auxiliary rocket motor. During 
the first ten milliseconds or so the acceleration 
of the 8.5 km/sec ramjet increases because this 
auxiliary rocket is operating and on-board mass is 
being expelled to accelerate the vehicle to the 
point where the ram pressure is sufficient to 
effect ramjet operation. This decrease in mass 


causes the acceleration level to increase with 
time. During the interval from 10 to 50 mil- 
liseconds, ramjet operation occurs and the vehicle 
acceleration is constant. Near the 50 millisecond 
point the vehicle velocity gets sufficiently close 
to the thruster exhaust velocity so that the ramjet 
can no longer produce adequate thrust.* This means 
on-board propellant must again be used and 
acceleration again increases with time. This use 
of on-board propellant continues until the desired 
final velocity (13 km/sec) is realized. 

For the 15 km/sec ramjet the final phase 
wherein on-board propellant is used is not required 
because the desired final velocity is significantly 
less than the exhaust velocity. As this discussion 
suggests it is desirable to have a ramjet exhaust 
velocity that is substantially greater than the 
desired vehicle velocity at the end of the launch 
phase. It is interesting to note at this point 
that constant exhaust velocity, constant thrust 
operation of a ramjet implies constant mass flow 
rate operation. This would seem at first to be 
difficult to achieve in this case where the projec- 
tile velocity is continuing to increase. A con- 
stant mass flow rate can be realized however by 
either varying the ramjet intake area or the dis- 
tribution of mass in the launch tube. 

Fignre 6 shows the variation in vehicle mass 
with time for each of the thrusters. These curves 
simply reflect the variations in acceleration 
observed in Fignre 5. It is noteworthy however 
that the 4,500 m/sec conventional rocket requires 
an initial mass of 180 kg to deliver a final mass 
of 10 kg at 13 km/sec while the 15,000 m/sec ramjet 
can do it with 11.5 kg and the ideal electromag- 
netic thruster requires only the desired 10 kg as 
an initial mass. It is also important to note that 
propellant tankage, structure and thruster masses 
are neglected in this analysis. If these masses 
were included in the analysis the differences 
between required initial masses would be even 
greater. 

CONSTANT THRUST 

CONSTANT EXHAUST VELOCITY 



* Application of the momentum balance equation 
demonstrates that the thrust (F = m [U-v]) 
goes to zero as the vehicle velocity (v) ap- 
proaches the exhaust velocity (U) . 
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The power profiles for each operating condi- 
tion are compared in Figure 7. The total energy 
(E) required to accelerate the payload to the 
specified final velocity for each option is 
obtained by integrating under the power curves and 
is also shown in this figure. The variation in 
total energy required for the launch is seen to be 
relatively modest between the various modes, but 
the ramjet and electromagnetic designs are seen to 
require substantially greater peak power levels. 
Figure 7 shows that the peak power levels of the 
high exhaust velocity ramjet and the electromag- 
netic gun are comparable. The electromagnetic gun 
is observed to require the least energy of any 
option, but it is important to note that this is a 
highly idealized model. The energy stored in the 
magnetic fields of this device, when the vehicle 
leaves the launch track, as well as joule heating 
losses will be substantial and they have not been 
included here. Their inclusion could well make the 
energy required for the electromagnetic gun more 
comparable to values obtained with the rockets and 
ramjets. 



The modeling that has been used in this study 
is highly idealized. For example losses due to 
aerodynamic and frictional drag and residual ther- 
mal energy in the propellants have been neglected 
for the launch phase, as have propellant pumping 
powers and losses associated with energy transport 
through conductors or windows. These could change 
absolute numbers significantly although the com- 
parative results are considered correct, and the 
suggest the hydrogen ramjet thruster could exhibit 
performance that is comparable to that associated 
with idealized electromagnetic thruster. It is 


therefore appropriate to examine specific ramjet 
configurations and to determine the mechanical com- 
plexity that might be required and the extent of 
the thermodynamic losses associated with their 
operation. 

RAMJET DESIGN FEATURES 

The conventional ramjet design shown in Fig. 
lb could be used in the present application. In 
this configuration it would fly down the launch 
tube processing a portion of the hydrogen and 
bypassing the remainder through the region between 
its exterior surface and the launch tube wall. 
Such a design has several drawbacks typified by the 
following examples of major design problems: 

• Transmission of electrical 
or electromagnetic energy 
through the ramjet sidewall 
into the heat addition core. 

• Control of the ramjet 
geometry to effect shock 
swallowing and efficient 
operation over the large 
vehicle velocity range that 
will be encountered. 

• Control of the radial posi- 
tion of the vehicle within 
the tube. 


A more suitable design would be one in which 
the heat addition zone was located external to the 
ramjet and the tube wall could serve to guide it 
down the tube. Figure 8 shows a schematic diagram 
of such a design which is designated the annular 
flow ramjet configuration. This axisymmetric 

design utilizes a launch tube of circular cross 
section. All of the propellant flow is forced to 
pass through the ramjet. The converging-diverging 
diffuser and nozzle passages are observed to be 
formed by the undulating configuration of the pro- 
jectile coupled with the essentially straight 
launch tube wall. The conceptual advantages of 

this configuration over the conventional design 
are 

• Beat addition can be 

effected outside the ramjet 
body directly from the tube 
wall. This eliminates the 

need for windows in the ram- 
jet for the case of beamed 
energy heating and sliding 
current contacts for arc 
discharge heating schemes. 

• Shock losses will probably 
be lower with the annular 
flow design because a coni- 
cal shock attached to the 
projectile would be expected 
while a bow shock would be 
driven ahead of the conven- 
tional ramjet. 

• Direct mechanical or fluid 
mechanical forces induced 
between the tube wall and 
projectile should tend to 
stabilize and center it. 
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• The structural design of the 
projectile appears to be 
relatively simple. 

• Variations in flnid flow 
areas needed to achieve 
proper flow conditions as 
the projectile velocity 
increases can be introduced 
by varying the cross- 
sectional area of the launch 
tube as a function of axial 
position. 



Fig. 8. Annular Flow Ramjet Configuration 

It is noteworthy that by eliminating the need to 
transfer electrical energy to the projectile itself 
and introducing the self-centering capability just 
alluded to. one mitigates two of the problems asso- 
ciated with such electromagnetic devices as the 
rail gun. These problems relate to the design of 
sliding electrical contacts and the maintenance of 
a tight tolerance on launch track straightness. 
Relaxation of the tolerance should for example 
permit longer track lengths and lower acceleration 
levels than those required for railguns. 

FLUID MECHANICS MODEL 


Since fluid mechanical forces propel the vehi- 
cle it is appropriate to begin an analysis of the 
annular flow ramjet by determining the magnitudes 
of thermodynamic and fluid mechanical variables 
required to produce the desired thrust level. This 
analysis must be considered preliminary and as a 
result very much simplified. The first major 
assumption in this analysis is that the process can 
be considered quasi steady. It is recognized that 
at the acceleration levels being considered this 
cannot in fact be accurate. However, including 
unsteadiness introduces major complications in the 
analysis. Estimates of the errors involved have 
indicated that though they are fairly substantial, 
they are not so large as to render these prelim- 


inary calculations valueless. 

Starting the projectile through the applica- 
tion of an external force will be addressed later. 
For purposes of this analysis we assume the projec- 
tile has been accelerated to a substantial super- 
sonic velocity. Though it is likely that there will 
be a conical shock attached to the nose of the pro- 
jectile, the strength of this shock can be con- 
trolled to some extent by the contour and half 
angle of the projectile nose. Consequently, for 
this preliminary analysis the diffusion process is 
modeled as an isentropic one. After diffusion the 
hydrogen working fluid undergoes heat addition at a 
subsonic velocity in a constant area duct and this 
is followed by an isentropic expansion to a super- 
sonic exhaust velocity. For this analysis propel- 
lant velocities were computed in the frame of 
reference of the moving projectile. The fact that 
the tube wall which is the outer boundary to this 
flow, is moving in this frame of reference was 
assumed to have no effect on the behavior of the 
flow. The basic equations of compressible fluid 
mechanics (conservation of mass and energy coupled 
with a thermodynamic path equation) were applied 
separately to the diffusion, heat addition and 
expansion processes. 

Thrust was computed by applying the momentum 
equation in the frame of reference of the control 
volume moving with the projectile. It should be 
noted that the boundary conditions associated with 
the annular flow ramjet are different than those 
associated with conventional ramjet operation. In 
the present case the inlet and outlet flow areas 
are equal and the flow area in the heat addition 
zone must be less than this inlet/outlet area. In 
the conventional ramjet case the static inlet and 
outlet pressures are set equal while in the annular 
flow ramjet the static pressure of the outlet 
stream will typically be greater than that at the 
inlet. 

The thermodynamic properties of the hydrogen 
working fluid will vary over the range of tempera- 
tures and pressures encountered during ramjet 
operation primarily because of the effect of disso- 
ciation. For this preliminary analysis however it 
was desirable to avoid the complications of tem- 
perature and pressure dependent specific heats and 
gas constants. Figures 9 and 10 show how the 
values of the ratio of specific heats and gas con- 
stant for hydrogen vary with temperature and pres- 
sure. In general, during ramjet operation, 
increases in temperature are accompanied by 
corresponding increases in pressure. This being 
the case Figures 9 and 10 suggest that it is not 
too unrealistic to treat these two parameters as 
constants. For the results contained in this paper 
ideal gas behavior was assumed and a ratio of 
specific heatS o (y) of 1.4 and a gas constant (R) of 
4160 joules/kg°K were used. 

Extensive analysis of the early phase of ram- 
jet operation when a bow shock was propagated ahead 
of the projectile was also undertaken. Modeling of 
this phase of operation was successful and showed 
that with proper launch tube design the bow shock 
would become attached as an conical shock without 
difficulty. However it also showed that positive 
ramjet thrust could not be produced at the subsonic 
diffusion inlet condition that accompanies opera- 
tion with a bow shock. It was therefore assumed 
that acceleration to supersonic speeds and initia— 
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tion of the isentropic diffusion process would be 
accomplished bp some alternative acceleration 
scheme . 



Fig. 9. Ratio of Specific Heats for 
Equilibrium Hydrogen 



TYPICAL SUPERSONIC LAUNCH CYCLE 

If one assumes that the projectile will be 
acted upon by either the ramjet force or some 
externally applied force needed to produce the 
desired acceleration, the projectile trajectory 
will be known. With the projectile velocity 
profile known the propellant flow velocity relative 
to the ramjet is known at each position along the 
launch track. One can now prescribe the diffuser 
throat area required at each position along the 
tube so that the flow will remain sonic at this 
throat. One can also prescribe the heat supplied 


in the subsonic heat addition region so that the 
specified thrust will either be produced 'or the 
mach number at the exit of the heating region (in 
the projectile frame of reference) reaches its max- 
imum allowable value of unity. After heat has been 
added, the flow is accelerated isentropically in 
the nozzle to the maximum supersonic mach number 
achievable under the constraint of the area change 
available. The static pressure at the nozzle exit 
is determined by this expansion and it will vary 
with projectile position. Proper control of this 
pressure at the base of the projectile would have 
to be achieved in an actual device by proper design 
of the launch tube. Such a tube would consist of 
two sections, a cavity section into which hydrogen 
processed by the ramjet would be collected and the 
launch section through which the projectile would 
pass. It has been assumed in this analysis that 
the cavity section could be designed to have the 
length and cross sectional area needed to accommo- 
date the propellant and facilitate maintenance of 
the necessary projectile base pressure. While it 
seems possible to do this intuitively, the examina- 
tion of the unsteady fluid mechanics describing the 
process has not yet been attempted. 

Figure 11 shows the results of such an 
analysis for the case where a constant thrust of 3 
x 10 Newtons is maintained on a 10 kg mass and the 
hydrogen in the launch tube is initially at 300°E 
(T ) and 30 atmospheres (P^) . The projectile 
selected for this example has a fixed 16 cm diame- 
ter at the diffuser throat and an 8 cm diameter in 
the heat addition region. These diameters are con- 
sidered adequate to assure the structural integrity 
of the projectile during launch. As Figure 11a 
shows the projectile velocity (in the launch tube 
frame of reference) increases linearly with time. 
Since the propellant is initially at rest, relative 
to the tube wall the projectile velocity is equal 
to the propellant inlet velocity in the projectile 
frame of reference. The exhaust velocity (in the 
projectile frame of reference) is observed to be 
equal to the inlet velocity during the first ~8 
msec of operation so during this time interval 
essentially no net thrust is produced through the 
mechanism of momentum change of the propellant. 
After about 16 msec of operation, however, the 
exhaust velocity has risen to about 15 km/sec and 
begun to level off. During the time between -8 to 
-16 msec all of the heat that can be added without 
choking the flow (achieving unity mach number at 
the end of the heat addition region) is being 
added. 

As Fig. lib shows the net ramjet thrust is 
inadequate to meet the 3 x 10° neuton thrust 
requirement at times less than 16 msec. Thus, an 
alternative accelerating scheme such as the con- 
stant pressure gun (discussed in the next section) 
would have to be used during this time interval. 
As the projectile velocity increases, more heat can 
be added and this causes the increase in exhaust 
velocity shown in Fig. 11a. After -16 msec both 
the exhaust velocity and the mass flow rate h«.ve 
risen to the point where the desired thrust level 
can be produced and the mach number after heat 
addition begins to drop below unity. 

Figure 11c, which indicates projectile posi- 
tion as a function of time for this case of con- 
stant acceleration, shows a launch tube about 400 m 
long would be required to accelerate the 10 kg mass 
to the terminal velocity of 15,000 m/sec in —50 
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msec. Figure lid indicates the variation in tube 
and projectile diameters that would have to be 
effected to satisfy fluid mechanical constraints 
imposed on the ramjet. As the sketch suggests, the 
projectile diameter at the diffuser throat and in 
the heat addition region have been treated as fixed 
at 16 cm and 8 cm respectively. The reduction in 
tube diameter required to yield unity mach number 
at the diffuser throat as the projectile velocity 
increases is shown by the solid line. Data from 
this line taken together with the position profile 
of Fig. 11c specify how the tube diameter would 
have to vary along its length. Figure lid also 
shows that the projectile diameter at the nozzle 
throat would have to increase as the projectile 
velocity increases and the input power is varied to 
induce the desired thrust level. The fact that 
this diameter must change during the launch cycle 
represents a significant mechanical complication in 
the system. It is noted however that the amount of 
change in this diameter becomes very small after 
about 20 m sec when the projectile velocity is ~5 
km/sec. It is anticipated that a pressure or tem- 
perature sensitive scheme might be designed to 
facilitate this change in projectile diameter at 
the nozzle throat. Alternatively one could hold 
the projectile diameter in the nozzle region fixed, 
reduce the tube diameter slightly more rapidly that 
Fig lid suggests and abrade the diameter of the 
projectile at the diffuser throat using protrusions 
from the tube wall. 


a, = 30000 G's 
m = I Okg 
T, = 300°K 




b. Thrust Profile 







Fig. 11. Typical Electrothermal, Hydrogen Ramjet 
Launch Cycle 

Figures lie and Ilf indicate how the static 
pressure and temperature at the nozzle exit should 
vary to conform to the problem specifications. The 
peaks in temperature and pressure occur at the 
point where the power input is a maximum (Fig. 
llg). This in turn occurs at the time when the 
thrust has just reached its design value and the 
mach number at the end of the heat addition region 
is still at unity. At this condition the exhaust 
velocity is still relatively low and the bulk of 
the thrust is being produced by the static pressure 
difference across the projectile. If one did not 
start using the ramjet until after ~20 msec into 
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the launch, we see that the exhaust pressure and 
temperature would lie in a very moderate range 
throughout the period of ramjet operation. Note 
that it is the pressure profile of Fig. lie that 
would have to be matched through proper design of 
the cavity section of the launch tube. This might 
be achieved by selecting not only the length and 
cross sectional area of this section but also by 
the inclusion of accumulator chambers that could be 
connected into the tube at the proper time. 

The temperatures of Figure Ilf are observed to 
be relatively low. So low in fact that it would 
probably not be possible to add energy to the 
hydrogen using electromagnetic radiation to couple 
to the plasma using the mechanism of inverse 
bremmstrahlung . It is noted however that this is 
not the case in the heat addition region where tem- 
peratures in the 10,000°K range needed to effect 
this energy addition mechanism do prevail. 

Figure llg shows the thermal power input pro- 
file required to effect operation as well as the 
actual profile of mechanical power input to the 
projectile. After the ~20 msec period over which 
the auxiliary starting force would be applied it is 
observed that the ratio of projectile power to 
thermal input power rises above ~50%. (An approxi- 
mate mathematical relationship for this ratio is 
given in the next section.) The ramjet is observed 
to become increasingly efficient as the launch 
cycle proceeds. This suggests that the ramjet will 
be more efficient at the higher velocities where 
the efficiencies of electromagnetic launchers drop 
off. The difference between the projectile and 
thermal input powers represents the power being 
deposited as kinetic and thermal energy m this 
hydrogen propellant. It should be noted here that 
frozen flow losses in the nozzle and electrical- 
to-thermal power conversion losses as well as fric- 
tional losses have been neglected in this prelim- 
inary analysis. 

Analyses similar to the one leading to the 
results of Fig. 11 have been carried out at other 
values of parameters such as the initial tube pres- 
sure and the diameter of the projectile at the 
diffuser throat. They suggest that inadequate mass 
flow rates are brought on by low projectile diame- 
ters and low initial tube pressures. This in turn 
leads to excessively high exhaust pressures and 
temperatures as well as extremely high thermal 
power demands. It has also been determined that 
the device will work if the nozzle is used as a 
diffuser and the propellant is exhausted at a sub- 
sonic mach number and a high static pressure. 


CONSTANT PRESSURE GDN - RAMJET COMBINATION 


energy would be added to a fixed quantity of gas at 
the base (in back of) the projectile and the pro- 
jectile would completely occupy the tube cross sec- 
tion. If acceleration is to be constant and the 
tube cross sectional area and projectile mass are 
constant, then the pressure in the tube must be 
maintained at a constant value as well during the 
acceleration. It is anticipated that the transi- 
tion from the gun to the ramjet phases of operation 
would be accomplished when the projectile passed 
through a diaphragm containing the hydrogen in the 
section of the tube to be used for ramjet opera- 
tion. The section of tube between the initial pro- 
jectile location and the diaphragm would be at low 
pressure initially so no significant amount of gas 
would be compressed ahead of the projectile during 
the constant pressure phase. 

A simple relationship between the required 
rate of energy deposition to the gas Q and the pro- 
jectile velocity v is obtained easily for the con- 
stant pressure (gun) phase of acceleration from the 
principle of energy conservation if frictional 
effects are ignored. This relationship is. 


Q = (m f o f g v) 


Since the rate of kinetic energy increase of the 
projectile is m^ g v the efficiency of the pro- 
cess (ratio of rateof kinetic energy increase to 
rate of energy supplied to driven gas) must be 
(y -l)/r. If hydrogen is taken to be the driver 
gas, y “ 1.4, the efficiency is about 28.6%. This 
is constant throughout the acceleration process. 
The energy not appearing as kinetic energy of the 
projectile must reside in the thermal energy of the 
driver gas. The temperature of the gas driving the 
projectile must therefore increase substantially 
during the acceleration. The magnitude of the rise 
is, of course, dependent on the mass of driver gas 
used. 

Once the velocity necessary for effective ram- 
jet operation is reached, operation of the projec- 
tile as a ramjet in the mode described in the 
preceeding section is initiated. It is of interest 
to develop an expression for the efficiency of this 
process and to compare it with the efficiency of 
the constant pressure gun process. Application of 
conservation of energy to the ramjet results in the 
following expression: 


• • u “V • 

Q = m — + m (h^-hp (5) 


n 


As discussed above efficient ramjet operation 
does not occur until some critical projectile velo- 
city is attained. Below this speed an auxiliary 
thrust device is required in order to maintain the 
constant acceleration level assumed. In the exam- 
ple whose characteristic parameters are shown in 
Figure 11 this velocity was reached at approxi- 
mately 20 milliseconds. At this point the projec- 
tile velocity was about 5900 m/s and the distance 
down the tube about 59 m. One way of attaining 
this necessary velocity would be to accelerate the 
projectile as in a gun and delay ramjet operation 
until this velocity was reached. In this process 


where the subscripts e and i denote exit and inlet 
conditions for the device, and h and m represent 
respectively enthalpy and mass flow rate. Combina- 
tion of this relationship with the thrust equation 
F “ m (D — v) + (P -P.)A, where A is the tube cross 
sectional area, Ind P is the pressure gives after 
rearrangement! 


q = FjU+v)( 


(P e- P i )A . " C p <V T i> \ 

F FOJ+v) / 


( 6 ) 
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Since F=m, 
written 


Q = m f a f gv 


g v, the 



expression above can 


(P e _P i ) . 


A + - 


F(U+v ) 
2 




be 

( 7 ) 


Comparing this to the rate of kinetic energy gain 
of the projectile, F = m^o^gv allows the definition 
of the ramjet efficiency. 


where v = (v. + v,)/2 is the average projectile 
velocity between t and t . Since m^o,g(x,-Xj) is 
equal to the change of kinetic energy of the pro- 
jectile, the overall efficiency for the ramjet 
acceleration process becomes 



(13) 
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F(U+v ) ) 
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For the example of Figure 11, = 77.31. For the 

combined constant pressure acceleration-ramjet 
acceleration process the efficiency is approxi- 
(8) mately 681. 


CONCLUSIONS 


The efficiency clearly varies with the projectile 
velocity v. At the beginning of ramjet operation, 
for the case given in Figure 11, at t = 20 mil- 
liseconds. 


1 

n R " (1.25 + 0.5) (1-.062 + .212) 


0.497 


( 9 ) 


and at t *= 50 milliseconds (end of acceleration) 


n R (.624 + .5) (1-.007 + .023) ” 0,876 * 


( 10 ) 

I 


These calculations indicate that the efficiency of 
the ramjet is always substantially higher than that 
of the constant pressure acceleration (0.286) and 
increases with increasing projectile velocity. 

The overall efficiency of the ramjet accelera- 
tion process can be estimated easily if the term 


(P - p )A 
1 ■ + 


C (T -T ) 
p e 1 

F(U+v) 

2 


is taken to be unity and if U is taken to be con- 
stant. The first assumption introduces a fairly 
small error of about 151 at the start of ramjet 
operation and a very small error of about 1.61 at 
the conclusion of the acceleration. The constancy 
of U depends on the ramjet parameters chosen. For 
the example of Figure 11 it is not a bad approxima- 
tion. With these two assumptions the expression 
for the rate of energy consumption becomes 


Q + h ) m f °f gv • (11) 

Integrating the above expression over any period of 
time, say between t^ and t^ 


Q = m, a. g(x 7 - x ) — + h , (12) 

|_2v 


Much analysis remains to be done to investi- 
gate completely the effects of all of the parame- 
ters of concern. In particular the effects of high 
acceleration and unsteady gas dynamic effects need 
to be considered and solutions that yield higher 
system efficiencies and/or more readily accommo- 
dated geometry changes need to be sought. Results 
presented in Figure 11 suggest however that the 
electrothermal ramjet designed to produce thrust 
after being given an initial velocity cannot be 
rejected on the basis that it violates a basic 
principle of fluid mechanics or thermodynamics. 
Moreover, on the basis of the same results it 
appears that the efficiency of the ramjet energy 
conversion process has the potential of being 
fairly high. 


REFERENCES 


1. Siegel, Arnold E., "Theory of High-Muzzle- 
Velocity Guns," appears in Interior Ballis- 
tics of Guns, H. Erier and M. Summerf leld 
eds., vol. 66, Progress in Astronautics and 
Aeronautics Series, American Institute of 
Aeronautics and Astronautics, 1979, pp. 135- 
175. 

2. Eolm, H. et al, "Electromagnetic Propulsion 
Alternatives," AIAA Paper 79-1400, Fourth 
Princeton/AIAA Conference on Space 
Manufacturing Facilities, May 1979. 

3. Snow, W.R. and G. E. O'Neill, "Construction 
and Testing of the 2.5 m Mass Driver," AIAA 
Paper 79-2095, Prince ton/AIAA/DGLR 14th Elec- 
tric Propulsion Conf., Oct. 1979. 

4. Jahn, R.G. , Physics of Electric Propulsion . 
Chapter 6, 1968, McGraw Bill. 

5. Hildebrand, R.B. "Aerodynamic Fundamen- 
tals", appears in Handbook of Astronautlcal 
Engineering . H. H. Eoelle, ed., McGraw Hill, 
1961, p. 5-24. 

6. Eerslake, W.R. , Private Communication of 
unpublished results, NASA Lewis Research 
Center, Oct. 1980. 

7. Wilbur, P.J., ''Ion and Advanced Electric 
Thruster Research" NASA CR-165253, Dec. 1980, 

p. 118. 

8. Eemp, N.H. od R.G. Root, " Analytical Study 
of laser supported Combustion waves in Hydro- 
gen," NASA CR-135349, Aug. 1977. 


138 



c 


DISTRIBUTION LIST 


Copies 

National Aeronautics and Space Administration 


Washington, DC 20546 

Attn: RS/Mr. Dell Williams, III 1 

RTS-6/ Mr. Wayne Hudson 1 

RTS-6/ Mr. Jerome Mull in 1 

MT/ Mr. Ivan Bekey 1 


National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 

Attn: Research Support Procurement Section 

Mr. G. Gol inski, MS 500-306 1 

Technology Utilization Office, MS 3-19 1 

Report Control Office, MS 5-5 1 

Library, MS 60-3 2 

Mr. N. Musial, MS 500-113 1 

Dr. M. Goldstein, Chief Scientist, MS 5-3 1 

Mr. T. Cochran, M.S. 501-7 1 

Dr. F. Terdan, MS 501-7 1 

Mr. R. Finke, MS 77-4 1 

Mr. B. Banks, MS 11-4 1 

Mr. D. Byers, MS 501-7 1 

Mr. W. Kerslake, MS 501-7 20 

National Aeronautics and Space Administration 

Lyndon B. Johnson Space Center 

Houston, TX 77058 

Attn: Mr. Hu Davis 1 

National Aeronautics and Space Administration 

Marshall Space Flight Center 

Huntsville, AL 35812 

Attn: Mr. Robert Bechtel 1 

Research and Technology Division 

Wright-Patterson AFB, OH 45433 

Attn: (ADTN) Mr. Everett Bailey 1 


NASA Scientific and Technical 
Information Facility 
P. 0. Box 8757 
Baltimore, MD 21240 
Attn: Accessioning Dept. 


1 



-140- 


Copies 


Case Western Reserve University 
10900 Euclid Avenue 
Cleveland, OH 44106 

Attn: Dr. Eli Reshotko 1 

DST 1 

Ministry of Defence 
Metropole Building 
Northumberland Avenue 
London, WC2 N5BL ENGLAND 

Attn: Dr. D. G. Fearn 1 

United Kingdom Atomic Energy Authority 
Culham Laboratory 
Abingdon, Berkshire 
ENGLAND 


Attn: Dr. P. J. Harbour 1 

Dr. M. F. A. Harison 1 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, MD 20771 

Attn: Mr. W. Isley, Code 734 1 

Mr. A. A. Vetman 1 

Dr. David H. Suddeth 1 

COMSAT Laboratories 
P. 0. Box 115 
Clarksburg, MD 20734 

Attn: Mr. B. Free 1 

Mr. 0. Revesz 1 

Comsat Corporation 
950 L' Enfant Plaza, SW 
Washington, DC 20024 

Attn: Mr. Sidney 0. Metzger 1 

Intelsat 

490 L' Enfant Plaza, S.W. 

Washington, DC 20024 

Attn: Mr. Rolland Schreib 1 

Rocket Propulsion Laboratory 

Edwards AFB, CA 93523 

Attn: LKDH/Mr. Walt Detjen 1 

LKDH/Dr. Robert Vondra 1 


DFVLR - Institut fur Plasmadynamik 

Technische Uni vers i tat Stuttgart 

7 Stuttgart-Vaihingen 

Allmandstr 124 

West Germany 

Attn: Dr. G. Krulle 


1 



-141- 


Copies 

Giessen University 
1st Institute of Physics 
Giessen, West Germany 

Attn: Professor H. W. Loeb 1 

Jet Propulsion Laboratory 
4800 Oak Grove Drive 
Pasadena, California 91102 

Attn: Technical Library 1 

Mr. Eugene Pawlik 1 

Mr. James Graf 1 

Mr. Dennis Fitzgerald 1 

Dr. Graeme Aston 1 

Electro-Optical Systems, Inc. 

300 North Halstead 
Pasadena, California 91107 

Attn: Mr. E. James 1 

Mr. W. Ramsey 1 

TRW Inc. 

TRW Systems 
One Space Park 

Redondo Beach, California 90278 

Attn: Mr. Sid Zafran 1 

National Aeronautics and Space Administration 

Ames Research Center 

Moffett Field, California 94035 

Attn: Technical Library 1 

National Aeronautics and Space Administration 
Langley Research Center 
Langley Field Station 
Hampton, Virginia 23365 

Attn: Technical Library 1 

Hughes Research Laboratories 

3011 Malibu Canyon Road 
Malibu, California 90265 

Attn: Mr. J. H. Mol i tor 1 

Dr. R. L. Poeschel 1 

Dr. Jay Hyman 1 

Dr. J. R. Beattie 1 

Dr. W. S. Williamson 1 

Dr. H. J. King 1 

Princeton University 
Princeton, NJ 08540 

Attn: Dean R. G. Jahn 1 

Dr. Arnold Kelly 1 



-142- 


Boeing Aerospace Co. 

P. 0. Box 3999 

Seattle, Washington 98124 

Attn: Mr. Donald Grim, M.S. 8K31 

Lockheed Missiles and Space Co. 
Sunnyvale, California 94088 
Attn: Dr. William L. Owens 

Dept. 57-24 

Electrotechnical Laboratory 

1-1-4, Umezono, Sakura-Mura, 

Niihari-Gun 

Ibaraki, Japan 

Attn: Dr. Katsuya Nakayama 

Sandia Laboratories 
Mail Code 4537 
Albuquerque, NM 87115 
Attn: Mr. Ralph R. Peters 

Ion Tech Inc. 

P. 0. Box 1388 
1807 E. Mulberry 
Fort Collins, Colorado 80522 
Attn: Dr. Gerald C. Isaacson 

EG & G Idaho 

P. 0. Box 1625 

Idaho Falls, Idaho 83401 

Attn: Dr. G. R. Longhurst, TSA-104 

Michigan State University 
East Lansing, MI 48824 
Attn: Dr. J. Asmussen 

Dr. M. C. Hawley 

The Aerospace Corporation 
Space Sciences Lab. 

P. 0. Box 92957 

Los Angeles, California 90009 
Attn: Dr. Y. T. Chiu 

The Takagi Research Laboratory 
Department of Electronics 
Kyoto University 

Yoshidahonmachi Sakyo-ku Kyoto 606 
JAPAN 

Attn: Dr. Toshinori Takagi 



1 


1 


1 


1 

1 


1 


1 



-143- 


Department of Aeronautics 
Faculty of Engineering 
University of Tokyo 
7-3-1, Hongo, Bunkyo-ku 
Tokyo, JAPAN 

Attn: Prof. Itsuro Kimura 

Prof. Tom Maul 

P. 0. Box 98182 

Tsim Sha Tsui Post Office 

Kowloon, Hong Kong 

British Crown Colony 

Mr. Susumu Masaki 
Department of Electronics 
Tokyo National Technical College 
No. 1220-2 

Kunugida-cha, Hachioji 193 
Tokyo, JAPAN 

Dr. Pradosh Ray 
Tuskegee Institute 
School of Engineering 
Tuskegee Institute, Ala. 36088 

Dr. John Barber 

International Applied Physics Inc. 
7546 McEwen Road 
Dayton, Ohio 45459 

Dr. Ron Hawke 

Lawrence Livermore Laboratory 
Livermore, California 94550 

Dr. Richard A. Marshall 
University of Texas 
Taylor Hall 167 
Austin, Texas 78712 

Dr. V. V. Zhurin 
Computing Center of the USSR 
Academy of Sciences 
Vavilova 40 
117333 Moscow, B-333 
USSR 

Dr. M. Krishnan 
Dept, of Applied Physics 
P. 0. Box 2159 
Yale Station 

New Haven, Connecticut 06520 



1 


1 


1 


1 


1 


1 


1 



-144- 


Dr. Ed Goldman 
Physics International Co. 

2700 Merced Street 
San Leandro, CA 94577 

Mr. An Shih-Ming 

Vice Director of Electric Propulsion Laboratory 
Space Science and Technology Center 
Chinese Academy of Sciences 
Beijing, China 

Mr. Lee W. Parker 

252 Lexington Road 

Concord, Massachusetts 01742 

Dr. Chris Olson 
Dept, of Physics 
University of Huntsville 
Huntsville, AL 35899 

Dr. Kevin Rudolph 
M.S. M0482 

Martin Marietta Aerospace 

P. 0. Box 179 

Denver, Colorado 80201 

Dr. Rolf Buhler 

Institut fur Raumfahrtantriebe 
Universitat Stuttgart 
7 Stuttgart 80 
Pfaftenwaldring 31 
West Germany 

Internal Distribution: 

Dr. H. R. Kaufman 
Dr. R. S. Robinson 


Copies 


) 


1 


1 


1 


1 


1 


1 

3 

3 



End of Document 



